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ABSTRACT. In this paper it is shown that, under certain assumptions about the 
metric of the space-time, energy momentum tensor of the total field, when Einstein’s 
field equations are satisfied, vanishes identically. 

A new version of the unified field theory is suggested, and it is indicated briefly that 
a genuine energy momentum tensor exists which is conserved and has a structure similar 
to that of general relativity. 

From the form of the energy momentum tensor and from the field equations of the 
proposed theory it is concluded that the law of force for a charged mass particle moving 
in a gravitational and electromagnetic field would have the expected form. 


no interaction of a moving charged mass particle with the electromagnetic 

field. It seems that Infeld’s assumptions about the identification of 
symmetric and antisymmetric parts @,, and 4, of the field variables g,, with the 
gravitational and electromagnetic fields, respectively, cannot be true even at 
sufficiently large distances from the centre of the particle. In a unified field 
theory matter is accounted for by an electromagnetic field, therefore it is necessary 
first to define a metric that represents space-time structure. Once a metric is 
obtained, one finds the matter tensor that corresponds to this metric which 
follows from the Hamiltonian derivatives of the Lagrangian with respect to the 
metric. This is the usual procedure adopted in all classical field theories 
(Maxwell’s classical electrodynamics, Born-Infeld electrodynamics (1934), 
general relativity). 

Thus, if m,, is the metric of the space-time, then the matter tensor is given by 

ee ee ee, (1) 
Omg Ox" Omyp, y 

In Einstein’s theory there is no obvious way for the choice of a metric; at 
least in its present form it does not render itself to such a treatment. But one is 
inclined, in comparison with the matter tensor and the metric of general relativity, 
to choose 6” =q%/4/(— Det 9%), to represent the required metric. If this is 
granted, then, by taking the non-invariant form of the Lagrangian (Schrodinger 
1948) of the total field, ie. B = 9° B,, where B,,=V3,0),— PP nas it can 
easily be shown that 


he his recent letter Infeld (1950) states that Einstein’s theory (1950 a) gives 


$B = — Ba + Biyoa",, wees (2) 
where #B?,=T',3— 531, 
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From (2) it follows that 


Te = es = 7 Bits = = hy/(—5)(B%2” + BP” — be) Rog + $0 OE Rag 
Qo Qo,y 


og LNs (3 
where b,,,b%” =8; and Det|d,,|=5 =Det|9#|. os 
On raising the indices of R,, by b”’, (3) may be written as 
Te = = V/(—B(RE-WR yy) +POOT Rage eee (4) 
Now, when the field equations R,,=0 are used in (3) we obtain 
TlSaOAF HAN We 9.7) Dee eee (5) 


From (4) it is obvious that what vanishes identically is the matter tensor of the 
field. If we were to look for the conservation laws (Schrodinger 1948) of the 
total field the result comes out to be 

3 (Boy — ae acl rere UO eee ame: (6) 

Thus the conservation of energy-momentum for the field is expressed only 
in terms of a pseudo quantity, and it does not contain a tensor term. On the 
other hand, it is an experimental fact that energy and momentum are not pseudo 
quantities, and (6) is not an acceptable representation of matter. Therefore we 
may safely state that the actual energy momentum tensor is a missing one in 
Einstein’s theory. A similar result as in (3) follows, when one assumes a,, to 
be the metric of the space-time as was done by Infeld. 

The above results indicate in a simple manner that a genuine energy-momen- 
tum tensor does not exist corresponding to the metrics 6,, and a,,, i.e. there is 
no source to create these ‘potentials’. ‘Therefore it is not surprising if one does 
not find the expected interaction of the moving particle with the electromagnetic 
field. An attempt to find the equation of motion in a unified field theory with the 
similar method of Einstein, Infeld and Hoffman (1938), can only be made after a 
metric and the corresponding energy tensor are obtained. It is not sufficient to 
associate various field quantities in an a priori manner with different fields, if 
one is not sure of the existence of the sources that create these fields. 

Now, we may draw two conclusions from the above considerations: 

(i) If Einstein’s theory is to be maintained, then the above assumptions 
which seem most natural in deriving an energy-momentum tensor for the total 
field have to be abandoned. 

(ii) That (1) is true finds confirmation also from Infeld’s result, since he 
assumes @,,, as in our second alternative, to represent the gravitational field. 

A unified field theory free at least from the deficiencies indicated briefly 
above is in the process of publication. For the sake of comparison we give the 
field equations of the new theory: 


Rip= Pasi Bea). oy ound Ge (7) 
Rast Riya? hep = ay (8) 
9F5 =00...° onli ie ad eee fe eee (9) 
rl=0, “(ge s=0)* 9) oe (10) 


where L,5,= yp.) +sy,.+%y,8 (Current density), and p has the dimension 
(length)'. The conservation laws of the total field are given by 


v!,=0, " |. (11) 
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where 
— ©) = (92 Ry, — $5;0"R,,) + (GP Ray — 3550 R,,») + 40°, Be, — SFB)... (12) 


‘The expression (12) in the absence of the antisymmetric se reduces to the usual 
mixed energy tensor (pseudo + tensor density) of general cca: When the 
field equations (7) are used (12) leads to 


Ba PCa Bish) + BH a) VBI) 
=F (GYR, — 3559R,,,) 5 (9° 7B oe — 7B), oan + i of (13) 
where we used the expression for 9 given as 


gt = a%a/(—g)— oo" 
‘The expression in square brackets in (13), when the cubes of the ¢’s are neglected, 
becomes }8°4""¢,,,—¢°",,(=T% say). The place of the matter tensor of the 
general relativity is taken in this theory by an electromagnetic energy tensor. 
If Einstein’s field equations are used in (12), then (11) reduces to (6). The 
right-hand side of the equations (7), with the same approximation for the ¢’s. 
reduce to T,,. The fact that we have a Maxwell’s energy tensor on the right- 
hand side of our field equations is sufficiently convincing for us to expect a correct 
law of force (Lorentz force) for a charged mass particle moving in a gravitational 
and electromagnetic field. ‘The same was already proved by Infeld and Wallace 
(1940) in general relativity. 

The details for the derivation of our field equations and conservation laws. 
from the Bianchi—Einstein identities for non-symmetric field (1950 b) as well as 
from a variational principle with a Lagrangian £=q”’R,, + 2p"[ \/(— 6) — /(—g)] 
will appear elsewhere. 

REFERENCES 
Born, M., and INFELD, L., 1934, Proc. Roy. Soc. A; 144, 425-451. 
EINSTEIN, A., 1950 a, The Meaning of Relativity (London: Methuen) ; 1950b, Canad. F- 
Math., 2, 120. 
EINSTEIN, A., INFELD, L., and HorrMan, B., 1938, Ann. Math., 39, 65. 
INFELD, L., 1950, Nature, Lond., 166, 1075. 
INFELD, L., and WALLACE, P. R., 1940, Phys. Rev., 57, 797. 
ScHRODINGER, E., 1948, Proc. Roy. Irish Acad., 52 A.1., 7. 


The Two-Body Scattering Problem with Non-Central Forces 
I—Non-Relativistic 


By G. J. KYNCH 
Department of Mathematical Physics, The University, Birmingham 


MS. received 6th July 1951 


ABSTRACT. A perturbation theory is developed for the two-body problem so that 
central and non-central forces can be treated in a consistent manner. A phase matrix S(r) 
is introduced which satisfies a first-order Riccati equation. ‘The asymptotic value of S 
occurs in the formula for the cross section, and for a special choice of functions is the 
Heisenberg matrix of the problm. The formalism leads to a simple comparison of 
neutron-proton and proton-proton scattering, an analysis of the degree of separation of the 
simultaneous equations which occur, and of the most suitable choice of S(r) for calculation. 


VelIN RO DUCTION 
HE object of this paper is to expound a consistent method of treating 
two-body scattering problems which can be used with tensor forces as 
well as central forces and leads directly to the values of the phases required 
in the calculation of cross sections. The method is a new type of perturbation 
6-2 


84 G. F. Kynch 


theory for degenerate states which is a generalization of a well-known equation. 
In the present paper we apply it to neutron—proton and proton—proton scattering 
with and without tensor forces; in the following paper (Kynch 1952) it is applied 
to the relativistic one-body and two-body problems. 

The radial wave function of a particle in an S state, acted on by a central force, 
can be written u=A sin (kr+8), where A and 6 are functions of 7, which tend 
to constant values as 7 tends to oo, provided the correct condition on A and 6 
is applied. The function can also be written u=(u,+%,S)a, where u=sin kr, 
u=cos kr and S=tan &(r). Our choice of condition on A and 6 is such that 
dS/dr =V(r)(u,+u,S)?, where V(r) is proportional to the potential. This 
procedure can also be carried through for other states. 

To determine the scattering cross section we need the asymptotic value of 
5 =8(00) asr+oo. This is found by integrating the equation with the condition 
S=0 at r=0. The two-body problem with central forces can be dealt with in 
the same way. 

When there are spin-orbit or tensor forces the two-body radial functions 
for a given total angular momentum J, M satisfy a set of three simultaneous 
equations. To obtain scattering cross sections it is necessary to determine the 
asymptotic forms of the three independent solutions of these equations. This 
we achieve by expressing the radial wave functions in the same form as before, 
where S(r) is now a square matrix whose elements satisfy a set of equations 
similar in form to that of the single S with central forces. These equations are 
independent of the amplitudes A(r). 

As a means of calculating phases these have the advantage of being first- 
order equations. They readily lend themselves to approximations, and to 
expansions in powers of the energy. Such an expansion leads in fact to 
linear equations for all but the lowest order term. ‘They also enable us to link 
together a number of odd results previously obtained. ‘Thus a relation between 
n—p and p-p scattering with the same nuclear potential was obtained by Chew 
and Goldberger (1949) for S states in the approximation where the coulomb 
term is small near the origin. It is now seen to apply to all states, even when 
tensor forces are present. 

Our equation can be modified so that d(r) rather than S=tan d(r) is 
calculated. For attractive central forces it is a monotonic increasing function 
of r. This modification is very desirable when there are bound states, as 8(r) 
then passes through the value 37 and sometimes higher angles where S diverges. 

Owing to the general method of approach we are able to attempt to find the 
best method of formulating the calculation with tensor forces in order to obtain 
as directly as possible the matrix needed in the scattering problem. It would 
seem that there is no especially favoured method when angular distributions 
are required. 

Heisenberg (1943) has introduced into quantum mechanics a scattering 
matrix S;=e". By a suitable choice of the free wave functions u, and u, our 
method gives a differential equation for a phase-matrix S(r) with the property 
that S;,=lim S(r). ‘The meaning of S(r) for other values of 7 is that it is the 

rT—> 0 
Heisenberg S-matrix if the potential were cut off at this distance. In fact S(r) 
is more general than Sy since it determines the potential as well as being 
determined by it. 

To find the energies of bound states it is possible to consider the analytical 
continuation of S}y=S(0o) in k-space, where k? is the energy. A more direct 
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method based on our equations is to replace u, and uw, by solutions for negative 

energies — i,” and choose ky so that the wave function has the correct asymptotic 

behaviour. If uj~exp—kor and u,~exphkyr then the necessary condition 

is lim S(r)exp2kyr=0. This is stronger than the condition (co) =0, which is 
T> co 

known to be unsatisfactory. For potentials whose solutions are known analytically 

the method is not very good, but this might be expected from the method of 


splitting up the wave function. 


[2 PERTURBATION: THEORY 2FOR*THESPHASE, FUNCTION 
Given two independent solutions u,(r) and u,(r) of the differential equation 


d du 
rataay at 5 sat El an a aed WG oS (1) 
we wish to find the asymptotic solution of the equation 
day rau 
5 {PG tqu= Vu bio sec (2) 


which is regular at the origin, i.e. w=0 when r=0. In the applications it is 
also true that lim rV(r)=0. 


T> @ 


The required solution of eqn. (2) can be written 
emia Gee ACT Reet apa seer (3) 


where the phase function S(r) and the amplitude a(r) both tend to constant 
values as r-> 00. 
We impose the condition 


da d 
ae + Uy 7 (Sa) =0 stravale).oca. (4) 
du (du, | du, 
so that Pi (FZ + 5) CORN ape gos ties (3) 
and substitute into eqn. (2). Using eqn. (1) and the Wronskian relation deduced 
from it, that 
: dus du 
p (uF — Us ms) a a eee Pena (6) 
where « is a constant, we can eliminate the derivatives of u,, vu, and a, leaving 
OL so Gag Se) Via UO )E. Sek seas (7) 


This equation is independent of a(r), so that an equation is obtained for S(r) :* 
a(dS/dr) =(uy+Sug)V(uy+ugS). == we anne (8) 
This generalized Riccati-type equation for the function S(r) can be solved 
numerically, or otherwise, and forms the basis of our subsequent discussion. 
A slightly more general form for u, 


TARY aS Sie: Oe > en BBA (9) 
dS dS A 
leads to the equation 0 (s, a — S, =) i ape rant AICS are ABE (10) 


where a=S,A and S,= SS}. 

* Some trivial details are included in the derivation, and the order of the symbols has been 
maintained, so that the same equations can be used formally to solve a set of simultaneous equations 
in § 5, where the symbols represent matrices. 
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§3. NEUTRON-PROTON SCATTERING WITH CENTRAL FORCES 
The scattering is determined by the asymptotic values of the phases 6 in 
the regular solution u~sin (kr—4L7+5) of the Schrodinger equation 


a ae (we ae) a —v(r)u. Sie ‘eo 0, (11) 


Following the method just described, we use for uw, and uz two solutions of 
this equation with V(r) =0. If we choose the solutions u, =) (kr) ~ sin (kr — 3 In), 
Uy =n, (kr) ~ cos (kr —4L 7), then S(r)—tanSasr—+o. The equation for S is 


RdS/dr=0(r)\Gn+Sape 20 0s ie Tee (12) 
Alternatively, by writing S=tan 8(r) (or S,=cos 8(r), S,=sin 6(r) in (10)), 
we find * : laa 
k d8/dr =v(j,, cos 6+, sin 8), (13) 
1.€: ké(r) = re o(j,cosStn,sind)dr. | 


The integral form uses the result that both S(r) and 8(r) are zero at the origin, 
and tend to zero so rapidly as to overcome the singularity in the irregular 
solution n,(kr)« r~” near the origin. 

Both these equations show that d(r) and S(r) are monotonic increasing 
functions of 7. The second is more convenient for numerical work when the phase 
is large. For S states in particular, where there is a bound state, 6 exceeds $7, so 
that S(r) must diverge for a finite value of r. Born approximation corresponds 
to neglecting Sn, in comparison with /,. 

An expansion in powers of the energy A=! tan d(r) =A) +A, R2+A,R*+ ...- 
leads to the equation for Ay, dAy/dr =v(r)(1+7A9)?. The other terms dj, Ag,...- 
satisfy linear first-order equations. Such an expansion is not possible if 8(r) 
exceeds $7 and a subterfuge is necessary, which is to write A=tan 7(7) and 
expand 7(7) in powers of R?, 


$4. PROTON -ERO TONS CAT LE RING 


For proton—proton scattering the wave equation with the coulomb term is 

— (e- SE) - F) n= Vu ee (14) 
where a =(2h/ Me”) is the coulomb length. In addition to the coulomb scattering, 
there is an anomalous scattering, due to the nuclear forces V, which is expressed 
in terms of an extra phase factor 8. 

In applying eqn. (8) we shall treat the nuclear potential as the perturbation, 
and use coulomb wave functions for u, and u. If we choose for u, and u, the 
functions defined by Yost, Wheeler and Breit (1936)} the asymptotic form of wis 
given by 
a(co ) 


u=(G,+SF,)a~ S. sin{kr —4$L7—y log 2kr +2 arg [8+ T(in + L+1)]}, 


where ka=(1/n) and S=cot8.  ——— (15) 


The equation for S is similar in form to that for n-p scattering, and only 
involves the nuclear potential, but the functions u4,=G, and u,=F, are 
completely different. 


* This particular equation was derived by Dr. Krook four years ago. It has since been 
published by Bergmann (1950). 


t In this section we use the notation of these authors. 
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A comparison of n—p and p-p scattering on the basis of our equation would 
be simpler if we could choose functions which behaved near the origin, where V 
is large, in much the same way as j,(kr) and ,(kr). Now there are solutions of 
the coulomb equation which go over to j,(kr) and n,(kr) in the absence of the 
coulomb force. From the equations of Yost, Wheeler and Breit (1937) we find 
Fr =(2L+1)CPd(kr); G,=%+¢{q—plog2n+plogr/a}, where p and q are 
functions of k®, and C,? is defined by them. The functions ¢ and % can be 
expanded in power series ¢ =(kr)”*! + higher powers, ys =(kr)-” + higher powers. 
The combinations of these which we are seeking are j,=4/{1.3...(2L +1)}, 
ny=1.3....(2L—1){(%+pdlog(r/a)}. Expressing F,, and G, in terms of these 
we obtain w=(j%,+n%S’)a’, where 


1/S’ = 17.3°...(2L—1)(2L4+1){(2L+1)C;? cot 5+q—p log (27)}, ...... (16) 
_ so that RaSsyar SV} iia Paks © e* Lae ee (17) 


We shall not discuss this equation in the detail in which we discussed the 
equations of the previous section. Everything said in the previous section 
applies in this section without modification. 

We continue our comparison of n—p and p-p scattering by assuming that 
the nuclear forces are the same, as in the charge-independent hypothesis. 
Eqns. (17) and (12) would then differ only in the use of 7{, and n{, in the former 
and j, and 7, in the latter. Now these two pairs of functions have similar series 
expansions and behave in approximately the same way near the origin where the 
potential is large, so that the function S’( 00) defined here* is almost identical with 
tan 6(00) of then—p problem. This result was obtained by Chew and Goldberger 
(1949) for S states. This independent derivation shows that it applies to all 
states, and it is shown later to be true even when tensor forces are present. 

A condition for this approximation can very simply be found, although it tends 
to underestimate the degree of equality. ‘The wave equation (15) can be expressed 
as an equation with no coulomb term but an effective angular momentum v where 
w(vt+1)=L(L4+1)+2r/a, ic. v—L=2r/(v+L+1)<2r/(2L+1). If we can 
neglect dv/dr then we can replace j{(kr) by 7,(kr) and n{(kr) by n,(kr) in the 
equation for S’. The equality of S’ for p—p interaction and tan 6 for the n—p 
interaction then requires that v—JL is small within the range of the nuclear 
force, i.e. that the coulomb term can be neglected inside this range. If the 
range is 7) then we require 27,/{a(2L+1)}<1. For a meson potential 


2r, 2h Me. M (=) 1 


ah ila PLOTS (Tie 


so that this condition is quite well satisfied. 


§5. PERTURBATION OF A SET OF SIMULTANEOUS EQUATIONS 
Consider a set of m simultaneous equations of the second order for n 


independent functions w,, u,....U,, of the type 

d [{ . du, 

5 (0%) Ca a (ies. I) cee eee (18) 
such that u,=0 ((=1, 2,....n) whenr=0. Any solution can be expressed as a 


linear combination of m independent solutions. 


* The function S’ was first introduced by Landau and Smorodinski (1944). 
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Let (w,); and (uy), be distinct solutions of the 7th equation, when all V,;=0. 

Then we shall show that any solution u(r) of eqn. (18) can be written 
Us={(Uy);+ Ei(U2)pSu(}al7), ste (19) 

where the functions S,,(r) are independent of the particular solution chosen. 

It is convenient to use a matrix notation. The solution u(r) and a,r) are 
written as (n x1) column matrices; the (u;),, (us); and also the q, in (18) as the 
diagonal elements of a square (7 xm) matrix; and V,, and S,;, as square matrices. 
Equations (18) and (19) then become elements of the matrix equation 


< (2%) +qusVu. o\ Pes fy) Sexes (20) 
and u=(u,+u,S)a. These are similar in form to eqns. (2) and (3), and to prove 
our statement (19) we carry through formally the process described in §2. All 
the equations (3)-(9) still apply, provided the quantities in them are regarded 
as matrix quantities. As already noted, the order of the terms in these equations 
has been maintained throughout. Eqn. (9) is 


aa dS/dr=(uy+ Sttg)V(uy+ugS)a.. sa wee (21) 


We now assert that this equation is independent of the functions a(r), so that S is 
determined by the matrix equation 


a dS/dr=(a,+Sus)V(uytunS), = aaanes (22) 


which is identical in form to eqn. (10). 

The proof of this is simple. We assume real solutions. Since the system 
of equations (18) is self-adjoint, it follows that S is symmetric, 

S = ),. 0 a ee re (23) 
where S’,;=S,, is the transpose of S. Thus S contains $n(m+1) independent 
functions. Assuming that w satisfies eqns. (20), it is possible to impose $n(n + 1) 
conditions upon these functions, provided that they are self-consistent and 
consistent with eqns. (20). Now it is clear that eqns. (21) impose just this 
number of conditions upon S(r) and that these are consistent with the original 
equations on account of the method of derivation. Hence the transition from (9) 
to (23) is justified. 

This result is particularly relevant to the scattering problem where the 
asymptotic values of a(r) are given by the problem. 

Before we apply the analysis to nuclear problems we shall discuss some 
alternatives. ‘The first is that a more convenient form can sometimes be obtained 
by the substitutions a= S,A and S,=SS, so that u=(u,S, +S )a by choosing 
S, suitably. ‘The second is that u, and uw, need not be real, and they can be chosen 
so that S is unitary. This point is discussed in §7. 

Amongst the remaining variations the most promising seems to be one which 
tries to diagonalize V(r), as this may diagonalize S. 'To do this we suggest a new 
trial wave function u = U(r)(u,+uS)a, where U(r) is an orthogonal matrix, and 
a supplementary condition so that du/dr=U(r)(du,/dr+Sdu,/dr)a. The 
resulting equation for S is 


dS 
a = (uy + Sua)[U-“(V —q)U + glu +usS)+pK, — «.....(24) 


where K = (uy' + Sutg')(U~! dU/dr)(uy + UgS) — (uy, + Suy)(U-1 dU /dr) (uy + uy'S). 
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The matrix q is already diagonal; so our aim of a diagonal S is completely 
achieved only if U-\(V—q)U and K are diagonal, and this is not possible. 
Moreover, to achieve any simplification we need, K=0 or U=constant. Nowa 
constant matrix U other than unity does not commute with q unless q has a special 
form, so that U can at best be chosen to reduce the off-diagonal elements of V. 
This can be completely done if all the elements of V depend onr in the same way. 


§6. SCATTERING WITH TENSOR FORCES 
To show the relevance of §5 to scattering problems we apply it first to the 


wave equation and then show that lim S(r) is in fact needed to calculate the 
'—> 00 


scattering. As these equations have been discussed many times in the literature 
we do not give all the details, but merely give the essential steps in a suitable form. 

Non-central forces affect the triplet, but not the singlet, state of a two-nucleon 
problem. ‘Triplet states of given total angular momentum J, M are derived from 
orbital angular momenta L=J,J +1. With non-central forces the last two are 
coupled. For an attractive potential 


V =0,(r) —2(L. S)v,(r) + {3(04.7)(09.7)/72 — 04.02 }03(7) ; 


the wave functions are of the form 
5 eee 


1 
hou = arr ' eal ujy(kr)f aN sete (25) 


where f¥y=Dp2_iciuY 1, m—mXm are the appropriate spin-orbit functions 


(cf. Rosenfeld 1948). The radial wave functions w” satisfy the equations 


Fag Ye Pe Ce ae 
ae, 


) 6f{S(J 4.1) 2 
fon —1) (0, Ph ae SLE 


dul ad (I A 
oa a (ie ( 2 ) yu = — {v1 +2(v_+ 03) }u7 


agli z (= ies ee) ee 
dr* os 


1/2 
=— {or +2 +2) Ce he i) burt es Oaur it 

In scattering problems we start with waves of given L and M, so that it is 
convenient to treat these three equations together. Comparing them with the 
equations in the previous section we see that S¥/?, is a 3 x3 matrix with only 
off-diagonal elements connecting L=J+1. The equations for S are (22) 
with the appropriate potentials, which we do not write in full. There is a single 
equation for S,, and three simultaneous equations for) Syzygy 9 y24.9o > 


Sy-1341 
See ene a ° ° . Se 
The scattering for given M is described by an incoming plane wave e'**y y 


and a scattered wave ,., which have both to be expanded in terms of the spin-orbit 
functions. If the scattered wave is #,,~kr-te™" ZX, , diy, fi, then the solutions 
for the radial wave functions have the asymptotic form 


Udy ~(i){4r(2L + 1) 2c sin (kr —4L77) + diye™™. 
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Now suppose we choose u, and uz in the equations (22) so that the assumed 
solution and its asymptotic form are* 


U= (uy +U,T)a~{(i)" sin (kr —$Lr)a,+ eT, ay}; 


then clearly the matrices a(co) and T'(0o) are such that 
diy = % Tyyay, a,=chy(tn(2L +1) Pp? sah. 


‘Thus a(o) is known and ¢%,, is determined by T(co). The scattering cross 
sections are given by 
Rdoy=dQ & | x diufin ? =dQ x | py Si) umes LCi uy (an(2L’ + 1) } 2? 
spins J, Imf<l J, 0,27 
and 
Roy=2|djyP= ate tn iy tim =D —3aty(T—T*) 7 pty 
o TIRES TE as ip ig 
since S,=1+42:T is unitary and symmetric. On averaging over M this can be 
reduced further to R8¢=4k?Xoy, =U ,37(2J + 1)Sp{—4(T— T*)}. For central 
forces all these results reduce to the usual forms. 

The expression for the angular distribution can be reduced on averaging 
over M, using relationships between the Wigner coefficients c4%, ; without tensor 
forces T is diagonal and the reduction has been made by Case and Pais (1950) 
using projection operators. Apart from this reduction it is clearly useful to be 
able to express T or Sy as the transform of a diagonal matrix, especially as the 
total cross section depends only on Spur 7, which is invariant under 
unitary transformations. Thus, if the unitary matrix S has eigenvalues 
exp 2187, (A=0 +1), then Spur {—4:(T— T*)}=%, sin? 8,. 

In the previous section we discussed one method of achieving this during 
the calculation of the matrix T(r). The potential was diagonalized so that a 
modified 7 was calculated whose off-diagonal elements were relatively small, 
the bulk of the effect of the non-central forces having been transferred to the 
diagonal elements and the transformation matrix U(r). 

Apart from the term K, the relevant factor in eqn. (24) for our equation is 
[U\(V—q)U+q] with q=L(L+1)/r?. When there is a tensor force this 
connects L=J+1. Let us first consider V alone. If the spin-orbit coupling 
factor v,(r) is a simple multiple of v,(7) (including zero) the matrix V is 
diagonalized by a constant orthogonal matrix 


[ez ( meee 4 beget aB 


—sin@ cos6é 


where 6 depends on the potential constants and J. As U is constant, K=0. 
‘This is very satisfactory, but off-diagonal terms still arise from U-1gU of the 
form (2J+1/r*) sin 26. ‘The precise effect of this term is difficult to say and 
depends on 4, but it is known that centrifugal terms L(L+1)/r? lead to large 
phase shifts in central forces, and the same may be true here. 

The alternative procedure of choosing U so that U-\(V—gq)U is diagonal 
bears some resemblance to that adopted by Christian and Hart (1950). This 
time @ varies with r, but K simplifies in view of the relation 


pad . dé zo \ ee 
—1.J)dr 


dr hay ah 


* We have used T rather than S to avoid confusion later. 
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The remaining factor in K is rather complicated. The whole expression is 
similar to that which would arise from a potential d6/dr and is, therefore, probably 
small. For example, if there is no spin-orbit potential and the tensor force 
U3 =(b«/r) exp —«r, where b has a value less than 2 (the value to give binding for 
a central force of the type), we find 20 = tan-1{2/2b«r/(ber —6 exp «r)}. Then @ is 
always negative, decreases from zero to about —}7 at r=1 and back to zero 
as r increases. ‘Thus these considerations are not inconsistent with the 
calculations of Christian and Hart, who found that the error in the phases when 
these off-diagonal terms were neglected was about 10°%, for states other than 
the SS state (J =1). 

There is need for a brief discussion of two more points. ‘The first is the 
telation between T (or S,) and the phase matrices S,, S,, which would be 
calculated using a trial function of the form u,=(u)S,+uUS,),7A”. If in 
fact we have chosen the free wave functions so that 


yet ~r, e— kr—L) its efkr Tee ~dg e~ ihkr— La) — lo efkr, 


then, after a certain amount of algebra, it is found that 
Sg = 1+ 20T = (uy Sy + wy So)(AvSy + AgS2)™- 

“Thus So — i if Si =| and By As 0, be Ay il 

In general, in calculations we choose S,=1, but in order to make a 
correspondence between the quantities calculated and the phases 6 which 
occur in central forces we have determined two representations where S involved 
three phases, corresponding to the three phases for L=J,/ +1 when there are 
no tensor forces, and a mixing parameter X. ‘These are given in the Appendix 
with the requisite equations. 


§7. THE HEISENBERG S MATRIX 
The S matrix is a unitary matrix which determines the scattering. It can be 
simply defined as a matrix S such that the asymptotic form of any stationary 


‘state is of the form eee he dye (Sf). ere CE (27) 


In a representation where J and L are diagonal, as for central forces, the 
asymptotic form of the radial functions u,(7)~sin (kr+6,— 3L7) shows that 
_ 4S is diagonal in L and has matrix elements Sy, =(—)* exp 275, 8;,,’. When 
there are tensor or other non-central forces S may not be diagonal in L, although 
diagonal in J. In fact the last remarks of the previous section show that 
SYP =(- Er PSo) rr: 

It therefore appears that our equations lead directly to a first-order differential 
equation for the S matrix, and that this equation is obtained by inserting for 
u, and uw, those combinations of j,(kr) and ,(kr) which have the asymptotic 
forms e‘!" and e‘* respectively. If these are h_(kr) and h,(kr) then the 
equation for S is 2ikdS/dr =(h_— Sh,)V(h_—h,S), where S=1 at r=0. Thus 
for the singlet S state with potential v(r) the equation becomes 

2ikdS/dr =v(r) fexp(—ikr)—S exp(tkr)}. ss. (28) 
Since this equation is true for any potential it follows that S(7) is unitary for 
all values of r. 

This unitary matrix is in some ways more desirable than the real matrices 

considered in the previous section, and for this reason it is interesting to notice 
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the following result. The matrix S(7) can by analytic continuation be evaluated 
for theimaginary values k = — kp, and for this value of k the equation (28) becomes 


2k,dS/dr =v(r){exp(—ikyr)—S exp(tkyr)}? sss ees (29) 
with S=1 when r=0, which is a real equation. If an analytic expression for S 


as a function of k, can be obtained then the S matrix for scattering at positive 
energies is found by replacing ky by 7k. 


§8. BOUND STATES AND THE S MATRIX 

The discussion at the end of the previous section is clearly related to the 
existence and the energies of bound states. We can deduce these in two ways, 
and we illustrate these for the singlet S state. 

The first method uses analytic continuation. As usually applied, the 
Heisenberg S matrix is calculated, i.e. S(oo), as a function of k. By analytical 
continuation this is extended to the imaginary values k= —zk,. For these 
imaginary values the energies of the bound states E = —h?k?/M are obtained 
by finding those values of ky such that S(0oo) =0. 

We see how this condition arises by considering the wave function which 
we assume to obtain eqn. (29). It is 


u={exp(—Ror)—'S exp(kyy)ja,, ~ ie © eae (30) 
together with a subsidiary condition on the functions a and S, so that 
du/dr = {—ky exp (—Ror)— SRo exp (Ror) fa. The condition S(co)=0 then 
eliminates the exponentially increasing terms in the brackets. 

This condition is clearly not sufficient. A more stringent condition has to 
be applied: lim,.,,, S exp 2kjy=0. ‘This can be seen from eqn. (30) or the ratio 


1du _ 14+ S exp 2Ror 
udr  - °1—Sexp 2k,r’ 
which eliminates the unknown factor a(r). 

A second method, which eliminates the necessity of knowing S as a function 
of k, is simply to integrate eqn. (29) directly with the conditions S=1 when r=0: 
and lim, ,,, S exp 2kjy=0. ‘The equation is a little difficult as the conditions 
on S and the fact that it continually increases with r show that S(r) diverges 
for some value or values of r and, in fact, the number of divergences equals the 
number of nodes in the wave function, including that at r=0. 

Attempts have been made to find a systematic method of solving the 
equation, not for S, but for y=S exp 2kor, which is 2kgdy/dr—4k,?y =v(1 —y)?, 
but without success. 

§9. CONCLUSION 

It is concluded that this method, which consists in taking seriously a 
well-known method of obtaining a first approximation to a single wave equation 
and then generalizing it, not only solves the equation in a reasonable form but is 
capable of producing useful general results. Its main disadvantage is that it 
requires a representation in which the angular momenta are diagonalized, but 
this may be related to the fact that we are working in ordinary space. It is 
significant that, when the momentum is diagonalized, which is desirable when: 
the incident waves are plane waves, and momentum space is used, the equation. 


for the state function is very similar in form to that given by the integrals of our 
differential equations. 
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APPENDIX 


The following summarizes (i) the wave functions, (ii) the equations needed 
to give the matrices S, and (iii) the matrix S) needed for the scattering cross 
section, with tensor forces. 

The free wave functions used are u,=j,(kr); u,=n,(kr). 

The diagonal element corresponding to L=J is omitted. 

The 6 and A in A are not the same as those in B. 


Method A. uw=(u,+u,S)a: S= ee ca) - (ia oF r ) 


Sess eas 
aa e251 _ Nu (1 + e241) 2ip 
(Nid Qin e212 Qu(1 + 2402) 
pS = (uy + Sua) V(uy +095) 
ANu 


— 2161 2762 
ay Tt Lee. 


Method B. u=(u,S,+u.S,)A: 


ae cos 6; A(cos 5, — sin a) 
#™\X'(cos 6, — sin 5,) cos 5, 


( sin 6, A (cos 6 — sin ay 
Sa= 


X(cos 6, — sin 6,) sin 55 
eS e201 ui (i — e242) Fy 
ae 6 e242 ui(i— e482) 
ds. dS } ; 
k (sy i 23g So’ =) = RQ = (S, uy + So Uz) V(uyS, + US») 
0 os. X'(cos 6, —sin 6,)(cos 6; —sin ~) 
7 Gres 5, —sin 5,)(cos 5, — sin 6,) Dai 


PA _)2(j—e2i91)(f— e202), = aAldr, 84! =d3,/dr. 
1+ pa 
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ABSTRACT. It has been shown that both for one particle scattered by a centre of force 
and in the two-body problem, treated relativistically, equations for the S matrix of the 
Riccati type can be found exactly similar to and as simple in form as those of the non- 
relativistic approximation. The analytic continuation of this matrix, to determine bound 
states, has the same limitations as in the non-relativistic problem. 


§1. INTRODUCTION 

HE perturbation method used in the previous paper (Kynch 1952, to be 

referred to as I) showed that the scattering problem could be reduced to the 

calculation of a phase-matrix S(r) from a set of first-order simultaneous. 
equations. This suggests that a similar technique could be applied to the 
relativistic equations of Dirac. Existing calculations for one particle using these 
equations (e.g. Mott and Massey 1949) have usually reduced the four equations for 
a single particle to a single Schrdédinger equation, by the elimination of all the 
components of the wave function except one. ‘This equation is then solved by 
any of the usual methods. 

We have found that the assumption of a wave function similar to those used in 
I results in a great simplification in the solution. An equation for a phase function 
S(r) is obtained which is not essentially different from that obtained for the non- 
relativistic calculation. ‘The potential appears in the same way, and so do the free- 
particle wave functions. The alteration is due to an extra term arising from the 
relativistic coupling of the states L=J +4. The potentials are only assumed to be 
invariant under rotations and inversion, in order to conserve the total angular 
momentum and parity. Apart from the form of the equation, which makes it easy to 
discuss relativistic corrections, one feature of the method is that the phase function 
SS = tan 6(r) is the same for all components of the Dirac equation. 

In the two-particle problem of given angular momentum //, M, in the absence of 
a potential there are four states which in the non-relativistic limit reduce to the 
singlet L =J, the triplet L =J and the two triplet states L=J +1. Whenthereisa 
potential satisfying the conditions mentioned in the previous paragraph, the last 
two states are mixed. Since the particles are free before and after scattering, we 
aim at the calculation of an S-matrix for these states at most yielding matrix 
elements between the last two, entirely analogous to that obtainedin I. This aim 
has been achieved, and the required equation appears as a matrix product. There 
is some difficulty here because of the occurrence of non-differential relations 
between the components, but this can be overcome. 

One difference between the solution of the Schrodinger and the Dirac equations 
is that the former required the introduction of auxiliary conditions before 
the solution was completed, but the latter does not. The general conclusion 
seems to be that the present method has advantages over former methods for the 
relativistic equations, both at low energies where corrections to the Schrédinger 
equation are required and at high energies where a rigorous calculation is necessary. 
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S27, DIRAC PARTICLE INS TATICNFIELD 
The Dirac equation of an electron in a static potential can be written as 
{—ip,;(o.V)+xpstE+uy=0, creas. (1) 

usingunits 4=c=1. Inanyotherunitsthe potentialis V=vic. Themomentum 
—tV =p has eigenvalues such that E? = x? + p?. 

The solution of this equation in polar coordinates has already been discussed by 
Darwin (1927) and Dirac (1932). ‘To derive the solution we use the relation 
”(o.p) =(o.r)(p.r) —(o.r x L), so that 


| <i. {oe 5 —AlerxL)} +49, +8 | y= —o. ny Tobe (2) 


We now introduce a representation in which J = L+4o and M=/, are diagonal. 
Their eigenfunctions f43;7 are constructed in the usual way from the spherical 
harmonics Y,, ,, and the eigenfunctions «, 8 of o? and Ox For each J, M there are 
two eigenfunctions f+ and f- obtained from L—J = +4 0r —4} respectively. With 
the usual choice of phases in these functions Hine matrix Sianeane of (o.7r) are given 
by (0.7) ft =rf-,(o.r)f-=rf*- 


The matrix elements of (o .r x L) follow easily from the equation 


BEX A vA et ag Et 8 a ol Sere (3) 
They are 
—1i(o.rx L)p 441 =—L(o.r)r .n41 ah Gat X L)7 y(t I) 7) peas 
so that 


i 0 ah 1 Otel: 
[eng se-rxD)| = F-sruls 5) 

1 0 1 Oat 
[log slerxt)| =F @-esna(g + * ) 


Finally we choose p? and p, diagonal, (p=p,, pz, p3) with eigenfunctions 
xe(B=+1asps=+1). Then for given J, M we write nb= LD, ,_.1u8 f” xp- 
The representatives wu’ being written as a column matrix uw, eqn. (2) is a matrix 
equation for su if v is replaced by a suitable matrix. 

The potential V is assumed to commute not only with J and M but with the 
parity. The equation then breaks up into two representations of different parity. 
The two pairs of equations are (cf. Darwin 1927) 


—tkD ,.u, +(E—«)ut = —(vu)t | (4a) 
Sul (her us = S(eu) PP On eee 
kD u=+ (B+) ui = ~(0W)} 

—ikD_ux +(E-«) uz = —(vu)= | Prt a 


(The upper suffix is the value of 2(L —/) and the lower se value of ps.) The 
abbreviations D,, =k 1{0/ér —(L+1)/r}, D. =k \{0/dr + L/r} are used, operating 
one’. 

In the absence of a potential these equations can be solved in terms of Bessel 
functions. Asin I, we use 


Jr lx) = / (Ara)? ST 7 5 (x) ~sin (x —3L7) 

m1 (x) = (g7x)"? J_(14.4) (*) ~Cos (x — 3 Lz) 
with the properties, if «=r, that D,j,= —Jz41,D_J;,=Jr,-1. Clearly we assume 
for a regular solution u~ « j;_,(kr), u* ax j7.,(kr). The equations are consistent 
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if LE? =x? + k? and we find 
wtmp ria (kr)iy we Sy per) py. 0a ears (5a) 
aD=Jyo, (Rr) Ul = SHOGg Rr). etek (55) 
There are other solutions involving the irregular functions 1, instead of j,. The 
factor 2 =(E—x)/(E +x) serves to distinguish small terms in the non-relativistic 
approximation. 
For eqns. (4a) we now put, as in our previous procedure, 


ut =(jyyyt+Sny.,)a=Uta, uy =i(jz_,+ Snz_,)a=Uya. 
This assumption is consistent with the equations, and eliminating da/dr and 
using the relation j,m,,,—"zj,.,=1, we derive the following equation for S, 
without any further assumptions : 


aS =U+*(oU)* + Uz*(OU);. eee (6) 

A similar assumption for eqns. (46) leads to a similar equation for S. 

To appreciate the simplicity of this equation it is best to assume central forces, 
with no p or o dependent factors in the original expression. _ On insertion of the 
expressions for U, eqn. (6) is 

a =(e+K)o(jz_y Snz_3)? + (C—*) OG zagt Si zyyy = 220s (7) 

The non-relativistic approximation is obtained for e—« and 
(e+K)v—>2«v =2mV /h’, and is identical with the equations of I. The corrections 
to this due to the use of Dirac’s equation are, therefore, the inclusion of the second 
term and the replacement of 2x by (e+«). | Both these corrections are of order A?, 
and both increase the magnitude of the phase compared with its non-relativistic 
limit. 

The appearance of the second term is, of course, related to the derived spin- 
dependent potentials which appear when the transition from the Dirac equation to 
the Schroédinger equation is carried through. 

Born approximation is obtained by neglecting terms in S on the right-hand side 
of the equation. ‘The exact equation shows that S is then underestimated, but 
that, in so far as Born approximation gives S' (0) finite in the limit k—-> 0, so also 
does the exact solution. 


§3. RELATIVISTIC TWO-BODY PROBLEM 
The two-body problem where both particles satisfy Dirac’s equation has been 
examined in considerable detail by Kemmer (1937). After eliminating the mass- 
centre coordinates and momenta the wave equation is 


(pot — pio"). p+K(p3 +p3) + 2e + v}p =0 
where 2e is used for the energy to eliminate factors of 2 subsequently, p = 3(p — p®) 


andr =r‘ —r (the neutron and proton mass are assumed equal). Using the sub- 
stitution mentioned in the previous section for r?(o . p) we obtain 


| - Halo.) 5, & sa(orxD)h— hind =@.1) 2 = n(@.rx1)} 


+aps-+ 0420 |=0 ae tes (8) 


where o=0% +0", ¢=o0%—o", and similarly for p and p. Owing to the relation (3) 


the same relation exists between the matrix elements of (o.7) and (o.r x L), or 
(o.r) and (o.r x L) as held in the previous section. 


The Two-Body Scattering Problem: II 97 


We now choose a suitable representation so that the total angalan momentum i, 
is diagonal. The 16 elements in; separate into four groups (p) = + 1, p? = +1, say, 
if these are diagonalized), each containing four rows which can be chosen so that J 
is diagonal. 

Let «,8 be eigenfunctions of o, and a, 6 eigenfunctions of p;. Then any ib 
is expressible as a sum of multiples of such products as aN b’ aN a’. The states 
of neutron and proton can be combined to give eigenfunctions of the total spin : 
Xm(m=0,+1) for the triplet, and y,y for the singlet. Combining these with 
suitable angular functions, we obtain, as in I, normalized eigenfunctions of the total 
angular momentum J, M: 


Triplet: Sia) = 2h Yt,u—mXm —(LE=J =0, +1) 
Singlet : fim =Y 1, Xs 
(Except for normalization, f-, f°, ft are X, B, Zin Kemmer’s notation.) 
In Soe to ee onelize his potentials Kemmer now uses a representation for the 
1, So that p3 ps and p} p; are diagonalized. We prefer to use that in which p; and p? 
are diagonal, as the large and small components are conveniently separated, i.e. we 
use eigenfunctions C, .;, €s defined in peer of a, b in the same way as y,, Xs are 


defined in terms of «, B. The operator p; py then mixes the states ¢,, ;. 
The wave ee for givenJ, M is expanded as* 
r= Zunfl, («,B=0, +1, S) 
where some suffixes.J, M have been dropped. 

As found by Kemmer, there are three non-zero representations, one of parity 
a oke containing six of the uJ, two of parity (—)’ containing four each. The 
functions wv, u)} are identically zero. 

The wave equations are easily written down using the abbreviations D., and D_ 
as before and « = {2(J +1)/(2J + 1)}2k, B = {2 /(2J +1)}2R. 


iBD (ux +u_)—ixD_(uz +ut)+ 2eug = —(vu) 
iaD .(u— —u~)+iBD_(uz —ut)+ 2eus = --(vu)3§ 
I(a) iz D_u} +iBD_u + 2(e—«)u= = —(vu)= 
iBD us —ixD ,u% + 2(e—K)uy = —(vu)> 
—ixD_u? +1BD_u +2(e+K)u* = —(vu)* 
—iBD ue —ixD us + 2(e+x)ul = —(vu)t 
iBD uy —iBD,.ug + 2eu%, = —(vu)e 
taoD_u, —txD_uf +2eu® = —(vu)* 
1(6) izD ,u® +¢BD_u*. +2(e—K)uy = —(vu)y 
—iaD,uS —iBD_uS. +2(e+K)ugy = —(vu)e 
—iaD (uy +us) + 2eu. = —(vu)®, 
11(b) iBD(ug+us)  +2eu®. hey (vu)? 
iBD, u°. —ixD_u®_+2(e+K)uy = —(vu)g 
iBD., u°.—ixD_u® +2(e—K«)us = —(vu°s 
In these equations we have expanded wf = &,, ,(vu)! f*¢,. 


* The following matrix elements are used : 
(c.n)ft=—28rf ; (c.r)\f-=—2orf ; (o.r)f°=—2r(af-+Bf*) 
(c.r)ft=2arfS & (Gon S SS (Co p= AD) 
eP=(J+1)(2I+1); P=J/(27+1). 
PROC. PHYS. SOC. EXYV,. 2—A 7/ 
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Potentials of the type assumed by Kemmer at most connect states u* and uw”, so 
that it is clear that certain components are then multiples of others, the ratios 
depending on the potentials, e.g. in II (b) the two components us and uw; . 

The solution for free particles in terms of the modified Bessel functions are : 


Gee wer ut = —PXa2—PYj_, uf=iadjo, 
: i= 0y ut =a pis, USER. 
(2) (i) La OF = = OBA Ts up = —1BAy » 

uy =J+> ut =$(ar— BP), , — Us = taj - 
I(d) Uo =Jos ug = —X%Fjo, 

w= 1aNf IBN 
II (d) UZ =Jo > ug =NJo , 

u’=—taAj,, u° =—iBd_. 


§4. SOLUTION WITH POTENTIALS 

The aim of our method of solution is to obtain an equation for a single function 
S(r), or a matrix, which only contains S, the free wave functions and the potentials. 
There is some difficulty if we try simple substitution of the form u=(j,+ Snz)a 
because of the components of the wave function which are direct multiples of others. 
The procedure adopted here has the merit of simplicity and is easy to carry through, 
but it is not at all elegant. We shall carry it through for the solution II (d) and 
state the results in other cases. 

In the set of equations II (6) there are four wave functions which we separate into 
two parts ; one gives the asymptotic form of the type previously assumed and 
contains the function S(r); the other is to be proportional to the potentials 

ul, = Ula+ Wi, Us =jo+ Sno, US. = —taA(74.+ Sn), 
Uz =M(jo+ Sn), U9. =—ifNj_+Sn_). 
The arguments of 7, and n, are always (kr), wherek?=e?—x*®. ‘These reduce 
when S=0 to the correct regular solutions with no potential. 

There are six unknown functions, and four equations. ‘Two of the equations 
give relations between the W, the others are equations for Sand a. 'Two further 
conditions can be imposed on the W, which are used to eliminate their derivatives. 

To obtain an equation so that a(r) and its derivatives are eliminated we write 
the functions u(r) as a column matrix, ut as its adjoint (conjugate) row matrix, 
and the potential as a square matrix. 

The derivatives ofthe W are eliminated if Wy + Wg =0; W2 =aK, W°_=6K, 
where K is still to be determined. ‘The remaining terms in K can almost be 
eliminated. Multiplying the equations by U® and adding, the W only remain 
with the potentials aoe ds 


(e+K)? dr” 

From the relations satisfied by the W we find 
4eWy = —4eWy =(vu)s —(vu)s,  2eK =P(wu)® +a(vu)? 
from which it follows that 22W+W=W+tVu. These equations show that in - 


fact the W are proportional to a(r) and that eqn. (9) can be divided by a(r), thus 
leaving an equation for S alone. 


=U+Vuss maging 2 eee (9) 
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This procedure, though cumbersome, is straightforward and easy to carry 
through. The result would be more elegant if the W could be eliminated from the 
equations completely before inserting their values. ‘They influence the scattering 
although they disappear from the asymptotic form of the wave function. 


§5. GENERAL SOLUTION OF TWO-BODY EQUATIONS 

The procedure used in the previous section for the representation IT (5) can be 
used for the others, or, indeed, on the wave equation as a whole. It follows the 
following common pattern. Using obvious abbreviations, eqn. (8) has the form 

{(F0/0r+G)+xpg+2et+v}u=0, © ...e. (10) 
where F and G are operators, and u is a column matrix of the representatives 
wz in eqn. (8). 

When there is no potential there are four regular solutions each of the form 
u,=J,, where the J, are proportional to the functions j,(kr) with an appropriate 
value of L. ‘There are four more irregular solutions with 1,(kr) instead of j,(kr). 
These solutions are arranged in four columns to form matrices J, N. 

We now assume u=(J+NS)a+W, where S and a are (4x4) and (4x1) 
matrices, respectively, which are functions of 7, and the W have still to be chosen. 
The W are chosen so that their derivatives do not appear in the wave equation, 


i.e. FW=0. Now this equation has eight constant solutions, normalized and 
orthogonal : 

CAI MS NTA ES ent alt a SEED Ree 

Chae) ee 2, Upp, Teitmatoe i 

(Q/2eoyo ub aug = 

(a/2khay jeu wea, us = —B 

(/2w5) ig = —ug =1 

(2k, ap” hale 

(2) ug =1. 

(ws) i, ==. 


The last two involve wand u?, which are known to be zero components and can 
therefore be ignored. The most general form for W is, therefore, a sum of 
the remainder multiplied by six arbitrary functions A(r) of r : 


WitcoXss ibd erage eS 2am (11) 


Substitute these expressions in eqn. (10): it becomes 

rs +N — =(Sa)} + (G+xp;+2e)W = —vu. 
‘The equations are sii? oi by the adjoint matrix (J+ NS)*. 
terms in W are zero : 


aJ+  . aN+ 
(J +.NS)+(G+xp,+2e)W = — (= — ++ — a) W 


since F+*W= FW=0. The term in 0a/ér is (J + NS) ei NS) da/ér =0, as the 
multiplier is effectively the expectation value of F in a free state, which is zero. 
Also lei ye BN =e EN —2r(1+)7), where A®=(e—«)/(e+x), so that, 
- finally, if k?=e?—xK? 


— mi ae 


The following 


(eS aT CNS) ta: (13) 
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We now show that a matrix a(r) can be dropped from this equation because the W 
are proportionaltoa. The evaluation of the W starts from eqns. (11) and(12). It 
may easily be shown that GW and p,W are orthogonal to all w,. Consequently, 
since the w are normalized and orthogonal wtw =1, if eqn. ( 12) j is multiplied by 
wt, then 2eA=wtv(J+NS)a+wtcwa. This equation can be solved for the A, 
giving, for example, A=A(J+NS)a. Thus from eqn. (13) we obtain the 
following equation for the matrix S consistent with the wave equations : 
2a bi =e — 
Cemar 
In view of the discussion 3, the phase matrix S and the whole of the equation 
reduce to a single equation each for the S matrix element for the representations 
1 (6) and II (6) and a2 x 2 matrix for S for the representation I (a). 


(J + NS)+o(1 + wA)(J + NS). 


§6. SUMMARY 

The form of the solutions which have been assumed for the wave functions 
shows that the asymptotic values of S atinfinity are related to the phases 
required to calculate the scattering. For example, if S(oo)=tané, then 
Ji t+ Snz,~sin(kr—4L7+8)/cos §. Moreover, if we substitute S(r) —tan 6(r), then 
our equations show that 5(7) increases monotonically withry to its maximum value. 
In I we gave the necessary modifications if S is a matrix. 

It has been established that the relativistic scattering of particles can be deter- 
mined in exactly the same way as for non-relativistic scattering. This means that 
many statements of a qualitative nature apply at both low and very high energies. 
‘Thus it is unnecessary to discuss in detail many statements made in I, although 
there are some points which deserve explicit mention. 

The comparison between n—p and p-p scattering is now of little importance. 
At very high energies the Coulomb term is a small perturbation. If it is neces- 
sary, however, a comparison by the method used in I can be made for the one- 
particle problem; for the two-body problem it can also be done formally, but the 
coulomb term enters into the calculation of the W, which is undesirable in view of 
its behaviour at infinity. 

It was also pointed out that our equations needed only slight modification to 
become differential equations for a unitary matrix S(r) which was asymptotically 
the Heisenberg S matrix in the representation. This is still true in relativistic 
calculations, the modification being simply to replace (7+ Sm) everywhere by 
(h, — Shy), where h, and h, are asymptotically e~“*” and e’*’, respectively. 

The S matrix form naturally raises the question of calculating the energies of 
bound states, assuming the Dirac equation, by analytical continuation of this 
matrix in k-space. It was shown in I that a necessary condition for a bound state 
atk = —1k, of energy — ky? was lim,.. ,. S(r) exp 2Rgr =0. This agreed, in particular, 
with the wes given by Jost (1947) ey the condition S(co) =0 is sufficient for any 
potential which is cut off at a certain distance, since dS/dr =0 for such a potential 
outside this range. It is now clear than the same condition applies using the Dirac 
equations at all energies, and it is also probable that the conditions recently given 
by Pais and Yost (unpublished) also apply to all states. 

‘The appearance of the functions W(r) in the two-body problem may be due toa 
wrong treatment of the equation. Apart from this it would be satisfactory if in 
fact they were zero when the correct potentials were inserted, but, of course, this 
cannot be shown until the potentials are known. 


The Two-Body Scattering Problem: IT IOI 
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ABSTRACT. The energy spectra of neutrons scattered in the forward direction when 
different target nuclei are bombarded with high energy protons are calculated and 
compared with recent measurements. The calculation is subdivided into three stages : 
(i) as a crude analogue to the nuclear scattering, the scattering of a single neutron ina 
square well is calculated using Born approximation; (ii) single scattering processes are 
estimated classically using the results of (i) to allow for the momentum distribution of the 
target neutron, and allowing for the attenuation of incident and scattered beams; 
(iii) multiple scattering is estimated similarly. Several simplifying assumptions have to 
be made and double collisions only are estimated. 

The caiculated differential cross sections are in fairly good agreement with experiment, 
but the calculated energy spectrum falls off too quickly. The effect of neglecting certain 
processes, such as triple collisions, is discussed, as well as the applicability of the nuclear 
model used. 


Si UN PROD UWET TON 

HE object of this paper is the calculation of the energy spectrum of neutrons 

scattered in the forward direction when different target nuclei are bombarded 

with 171 Mev protons. Such spectra have recently been measured for 
beryllium, carbon and aluminium (Cassels, Randle, Pickavance and ‘Taylor 1951 a). 
A comparison of these results will indicate to what extent the nuclear model used 
is adequate at high energies. In particular, the experimental results show two 
characteristic features which should be explained. Firstly, the differential 
cross section is much smaller than N times the free particle cross section. 
Secondly, the energy spectrum has a maximum at about 150 Mev; for smaller 
energies it decreases only slowly. For example, the intensity of 100 Mev neutrons 
is still 67°% of the maximum intensity. 

The general features of high energy nuclear reactions are well known (Serber 
1947, Fernbach, Serber and Taylor 1949), but in detailed calculations considerable 
difficulties are encountered. There are no satisfactory nuclear models. ‘There 
are indications that nuclear forces are many-body forces, but even if one avoids 
these, there are no satisfactory two-nucleon potentials to account for the known 
free particle scattering data (see Christian and Noyes 1950, Christian and Hart 
1950, where references to experimental papers are given). ‘This difficulty can 
be overcome by expressing all processes approximately in terms of single-particle 
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collisions and then using the experimental data. Finally one is here dealing with 
a many-body problem, and only quite crude approximate calculations are 
practicable. 

Apart from these general considerations, to see the points of importance 
involved let us consider qualitatively the scattering of protons by nuclei. 

When an incident proton hits a nuclear particle many processes may occur, 
but on the Serber picture the most frequent result will be for two nucleons to 
be ejected from the nucleus. There is a much smaller chance at these high 
energies of only one nucleon being ejected and the nucleus remaining in an 
excited state; in any case such transitions will only affect the intensity of neutrons 
of almost the maximum energy (at most a few Mev below). As the target nucleon 
is bound in the nucleus it is no longer possible to assume conservation of energy 
and momentum for the two nucleons colliding. We assume conservation of 
energy, but any spare momentum can always be taken up by the nucleus. 

At 170 Mmev the mean free path of protons in nuclear matter is of the same 
order as the nuclear radius (A=2-8 x 10-8 cm), hence multiple collisions will 
matter. ‘Their effect is to cause an attenuation of incident and scattered waves 
and to produce a broadening of the energy spectrum. ‘These two effects should 
be responsible for the characteristics of the experiments mentioned above. 
The Coulomb scattering, being due to long-range forces, can be estimated for 
the nucleus as a whole and turns out to be negligible at these energies. 

The calculation is subdivided into three sections : 

(i) As a crude analogue to the nuclear scattering process that we want to 
consider, we calculate the scattering of a 170 Mev proton by a single neutron in 
a square well, of such depth that its binding energy is 8 Mev, using Born approxi- 
mation, which should be quite good at these energies. One can avoid using a 
specific (p,n) interaction potential and can express the results in terms of the 
known differential cross section for free (p,n) scattering. This calculation 
corresponds to the scattering from a ‘typical’ nuclear neutron (§ 2). 

(ii) Single scattering processes are estimated classically, using the results 
of (i) to allow for the momentum distribution of the target neutron and allowing 
for the attenuation of incident and scattered beams. ‘The nucleus is considered 
as consisting of individual non-interacting particles distributed uniformly 
throughout the nucleus. The Pauli principle, however, is taken into account 
by forbidding transitions to ‘occupied states’ in the calculation of the mean 
free paths (Goldberger 1948) (§3). 

(iit) Multiple scattering can be estimated similarly but leads to very heavy 
computation, and some simplifying assumptions must be made. We have only 
calculated double collision processes, allowing for the various possible processes 
of direct and exchange scattering (§4). ‘The effects of neglecting triple and higher 
order processes should not be serious. It will be discussed together with the 
results of these calculations and the applicability to high-energy scattering 
problems of the nuclear model used (§ 5). 


§2. PROTON SCATTERING BY A SINGLE NEUTRON 
IN A POTENTIAL WELL 
Let $(r,) be the wave function of a neutron in the Is state of a spherical square 
well potential of radius R and depth Uy, where for R we take the nuclear radius 
and U, is so chosen that the binding energy « of the 1s level is 8mev. Then 
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$(72) = Ary sinary, .72<R; (r2)=Brytexp —Br,, r,>R; A and B are 
determined from the normalization of ¢ and the continuity of its logarithmic 
derivative at r,=R. Describe the incident proton, energy E,, momentum /ik, 
by a plane wave exp(ik,.r,). In the final state both neutron and proton are 
ejected from the well, with momenta /ik,’, /ik,’, into elements of solid angle dQ,’, 
dQ,', with energy Ey’, FE,’ respectively. The differential cross section for the 
incident proton of energy E, to be ejected into the solid angle dQ,’ with energy 
in the range (F,’, E,’+dE,') is then given by 


/ / , M 2 g / / 
o( Ey, By’) dQy’ dB’ = Fe 5" | | Hi, |*8(By’ + By’ — Ey +6) 
1 
1 dk,’ 


x (anys Re dE,’ Glens? Ahyd@, cds Ie 4) evo (1) 
the integration being over k,’._ The 5-function ensures conservation of energy. 
H,, is the matrix element for the transition. If the (p, n) interaction is described 
by a Serber potential }(1+ P)V(r), (r=r,—r,), then 
A, = J$(r2) exp (thy .r,) {3 V(r) exp (—tky’.1r —tk,’. ry) 
+4V(r) exp (—zk,’.r, —7k,’.r,)} dr, . dr, 
and 
| Ais P = | Jo(rs) exp { —2( ky’ + ka — ky )ro} dir, [? 
x |4JV (7) exp {2(k, —k,’)r} dr + 4fV(r) exp {i(k, —k,’)r}dr|?. ...... (2) 

The first term on the right-hand side gives the momentum distribution of the 
neutron in the square well. If we were dealing with a free neutron, ¢(r,) would 
also be a plane wave and this first integral would reduce to the 6-function 
expressing conservation of momentum. Of the two Fourier transforms of the 
potential V(r), the first corresponds to exchange, the second to direct scattering. 

The differential cross section X(E,,26) in the centre-of-mass system, for a 


proton of energy £, to suffer an exchange scattering through an angle 25 when 
striking a free neutron, is given in Born approximation by 


2\ 2 
X(E,, 28) (=) = |4JV(r) exp {a(K —K’)r} dr + 3[V(r) exp {0(K + K’)r} dr [?, 
ae (3) 

where iK, —/K’ are the momenta of the proton before and after collision 


respectively (in the centre-of-mass system). In eqns. (2) and (3) the last terms 
can be neglected as we shall only consider final states in which a neutron is 
scattered in the forward direction. We can then identify these two equations 
if we make | k, —k,’| = | K—K’|, 1.e. define 6 by 


sind =|k, —ky’[/ky=1-(Ey/E)!@ — (0<8<4n). 


Using the Berkeley data for X(£,, 26), the calculation of o(£,, E,’), eqn. (1), 
is quite straightforward, though lengthy. If calculation is limited to the energy 
spectrum of neutrons in the forward direction, all the integrations can be carried 
out analytically except the radial integration in the Fourier transforms of 4(r,). 
o(E,, E,’) is, however, not quite the quantity required to compare our results with 
the measurements of Cassels et al. for the proton beam in the cyclotron is 
not monoenergetic but, due to multiple traversal of the target, has an energy 
spectrum given approximately by J(£,)=J,exp[(171—,)/15]. o(£,, £,') must 
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be integrated over this energy spectrum to obtain the differential cross section 
o,(E,’) in the forward direction per unit solid angle per unit energy for the exchange 
scattering process considered. o,(E,’) is plotted against E,’ in fig. 1, the radius 
of the well R having been taken equal to the radius of the Be nucleus. 

We shall, however, also use the same function o,(£,’) for scattering by Cand 
27Altargets. The error thereby introduced should not be serious, as all we require 
is a typical ‘smearing function’. Furthermore, the nuclear radius R (i.e. the 
well radius) varies only as the cube root of the atomic number 4. 


Au 


Neutron Spectrum (millibarns per MeV per unit solid angle) 


0 | | it 
60 80 100 120 140 (60 


Neutron Energy (Mev) 


Fig. 1. Energy spectrum os () of neutrons ejected in forward direction from a square well as 
a function of the neutron energy EL. 


We have also calculated the quantity analogous to o,(£,’) for (p, p) scattering. 
Though the absolute magnitude differs, the energy dependence is very nearly 
the same in the two cases. We shall need to make use of this result when 
considering double collisions. 


§3. SINGLE SCATTERING PROCESSES 

We next consider those processes where the incident proton collides with 
a nuclear neutron so that the neutron is scattered in the forward direction and 
emitted from the nucleus without further collisions, allowing for the attenuation 
of incident and scattered beams. 

Consider a nucleus with N neutrons and Z protons and assume that these 
are uniformly distributed with density m throughout the nucleus, pictured as 
spherical of radius R. 

If the incident proton, energy Ey, moves in the direction of the positive 
x axis of a coordinate system whose origin is at the centre of the nucleus, and it 
enters the nucleus at the point with polar coordinates (R, 0, $), then the differential 
cross section o,(£,) dk, dQ, for a neutron with energy in the range (E,, E, + dE,) 
to be emitted into dQ, in the forward direction, due to a single collision, is given by 


o,(£y) = { 


7—2Rcosd 


exp (— ¢/ly) 


6=n / C=0 


N ' 
x (Ga errs sin 4 cos 6 dé dé o(E,) | exp {(2Rcos 6+ €)/1,}. eee 


va 


i 
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Here the square bracket gives the probability of a collision of the desired kind 
occurring in the volume element 27R?sin@cos6d0d¢. The two exponentials 
give the attenuation of the incident and scattered beams. The mean free path 
/,=K£,) of a proton of energy E, in the nucleus is calculated from Tone oper TE 
free proton—neutron and proton-proton cross sections: 

3 
where «,,, and ,,, are numbers smaller than unity allowing for the fact that on 
account of the Pauli principle some transitions are forbidden (Goldberger 1948, 
Cassels, Randle, Pickavance and Taylor 1951b). All energies occurring in (4) 
are measured relative to a potential well of the same depth (U, =21 Mev) as the 
one used in $2. /, is calculated at the mean energy EF, of the incident proton 
beam, but ((£) is a slowly varying function of E over energy ranges of 10 to 20 Mev 
when E is of the order of 100 to 200 Mev. Apart from this simplification, we have, 
however, allowed properly for the energy spectrum of the proton beam, by using 
o({E,). Eqn. (+) simplifies to 


TEN Care Aas Arcee pes) ee Uae (5) 


pp 


3N 1 
o,(£;) cs Gs RiilEo» E,)o(E;), seca (6) 
2 1 2Ru 
where F(4, A)=7| wdu df exp, — ees : — Z BF ves Grane (7) 
ly Jo 0 ke ly fy ely 


is the probability of making exactly one collision with the scattered neutron 
moving in the forward direction with energy F,. 


§4. MULTIPLE SCATTERING PROCESSES 

Apart from the scattering processes considered in the last section there 
are two other types of processes to be considered. 

Firstly, there are processes in which only one nucleon is ejected from the 
nucleus which remains in an excited state. On a shell model where the A lowest 
levels are all filled up and the spacing of levels corresponds to a mean binding 
energy of 8 Mev, such processes lead to ejected nucleons with energy within 
a few Mev of the incident protonenergy. These processesare thus not of importance 
at more than 8 Mey below the proton energy. Furthermore it seems likely that 
their contribution to the cross section will not be large. 

Secondly, there are multiple collision processes. As a result of a collision a 
high-energy particle will in general be deflected through only a small angle into 
a cone of not more than 40° half-angle, retaining most of its original energy. 
(In an exchange collision it will of course change from proton to neutron or 
vice versa.) The target particle will gain only a small fraction of this energy. 
As we are only interested in neutrons which are ultimately ejected in the forward 
direction with more than 60 Mey, this being the lowest observed energy, we shall 
throughout assume that only the high-energy component of a collision, scattered 
through not more than 40°, will lead to observable neutrons. 

The calculation of multiple processes, i.e. processes in which the high-energy 
component makes several collisions before being ejected in the forward direction 
as neutron, is much more complicated than for single collision processes. For 
we must now also consider cases where the high-energy particle is deflected. 
We shall treat only double collisions, as adequate calculations for triple and higher 
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order processes are very heavy and do not seem worth while in view of the rather 
crude model used. The effect of neglecting these processes will be discussed in § 5. 

Double collisions in which a high-energy neutron is ejected can be subdivided 
into four groups according to the nature of the particles involved: (i) (p, p) 
collision followed by exchange (p, n) collision, (ii) direct (p, n) collision followed 
by exchange (p,n) collision, (iii) exchange (p,n) collision followed by (n,n) 
collision, (iv) exchange (p,n) collision followed by direct (p,n) collision. For 
' definiteness we shall consider a process of the first class and then generalize our 
result to allow for all four possibilities. 

Using the same coordinate system and notation as in §3, suppose the first 
collision of the incident proton with a nuclear proton occurs at the same point P, 
as the single collision in §3. Suppose o,,(Eo,«; E,)dQ,dE, to be the cross 
section for a proton with energy E, to be scattered through an angle « into a 
solid angle dQ,, with energy in the range (£,, E,+dE,), the velocity vector v, of 
the scattered high-energy proton having polar angles «, 8. Suppose the scattered 
proton, after travelling a distance s, makes a collision with a neutron in the volume 
element s*dssinadadf8 such that the neutron is scattered into a solid angle dQ, 
in the forward direction, with energy in the range (£,, FE, +dE,), the cross section 
for such a process being o,,*(E,,«; E,)dQ,dE,, and let ¢’ be the distance this 
neutron has to go before leaving the nucleus. Then the cross section for processes 
of this group is given, analogously to (4), by 


toa ff I ft 


s=0 


oe nZ exp (—s/L,) 
: {2nR® sin 6 dd dé E Pop Lon; zB) | dQ, aE, | 7a, = 
aed Dee pe a 

x <s*dssina da dB om (4&3 zs) | dQ, dE, > exp (—€'/l,). 


‘The term s° dQ, in the denominator allows for the divergence of the proton beam 
between the two collisions. The integration over s goes from s=0 to s=Sp, 
corresponding to the point where the velocity vector v, intersects the nuclear 
surface. In practice the ranges of integration of , and «are restricted by 
conservation of energy considerations. 

We obtain equations similar to (8) for processes of groups (ii), (iii) 
and (iv) with the appropriate factors replacing [(nZ/A)o,,(Ey,«; E,)] and 
[ (nN/A)oy,* (£,, a; E.)). 

The quantities o,,(E),«; E,) etc. allow for the dependence on the energy Ey 
of the incident particle and the resulting energy spread for the various processes. 
In §2 we calculated this quantity for exchange (p, n) scattering for the case where 
the emitted neutron is ejected in the forward direction. In principle we should 
repeat this calculation for all the processes considered obtaining the angular 
distributions. Such a calculation is clearly not worth doing, and we have hence 
made several approximations and simplifying assumptions. 

Firstly, we take for o,,,(Ey,«; £,) (and similarly for the other cross sections) 


Opp (Lo, %; Ly) =4 cos aX,,,(Eo, 2x)y( Ey cos? x — F;) E, <E, cos? « 
== () EF, >, cos? « 
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where &,,,(Eo, 2«) is the differential cross section in the centre-of-mass system for 
free (p, p) scattering through an angle 2x, and y(x)=N~10,(E,—x) where o, is 
the function defined in §2, and N is a normalizing factor so chosen that 


bau y(x)dx=1. The function y(«) can be represented analytically to a very good 


approximation by an expression of the form y= Ax? exp(—vx), A and v being 
constants. 

We thus assume that a proton with initial energy E,, scattered through an 
angle «, will have a maximum energy F, cos’ «, corresponding to the free particle 
case. This is accurately true for «=0, and at high energies it is approximately 
true for deflected particles. Furthermore we assume the energy spread to be 
a function of E,cos?«—E, only, but otherwise independent of the angle «. 
Without carrying out the complete calculation this seems the best and simplest 
assumption. In practice one expects a larger energy spread for larger angles 
of deflection, but the error made by our assumption should not be too serious. 

Secondly, it is desirable to introduce approximations in (8) so that all the 
integrations except that over « can be carried out analytically (for each value 
of E,). For this purpose we take EF, =(E)—25)cos’& in &,,,*(Fj, 2x), which for 
a=0 reduces to EF, =E,—25, the mean energy of the scattered proton for the 
energy distribution (9), and in J, =l(E,) we take FE, =(E,EF,)"*, corresponding to 
the case of free particle scattering. In (8) we put «=0 in sy) and @’. This 
corresponds to assuming, for the purpose of calculating the fraction of particles 
making exactly two collisions, that they are not deflected. Again, this assumption 
will not introduce a large error as we are dealing with scattering directed strongly 
forward, and furthermore, the effects of inward and outward scattering from 
different regions of the nucleus roughly cancel out. 

Finally, using the symmetry of &,,,, and the fact that &,,, (taken equal to &,,,) 
can be taken as constant over the energy range concerned, all cross sections can 
be replaced by either X,,, or X,,. In this we have assumed, as in §2, that the 
experimental cross section corresponds entirely to ordinary or exchange scattering 
for small or large deflections respectively. 


§5. RESULTS AND DISCUSSION 
We have calculated the differential cross sections in the forward direction 
for neutrons with energy FE of more than 60 Mev for single collision processes 
(o,(E >60)) and for double collision processes (o,(£ >60)) for beryllium, carbon 
and aluminium. These values, as well as o(f > 60) =o,(E > 60) + 0,(H >60), are 
given in the table. 


Theoretical cross sections Experimental cross sections 
(millibarns per unit solid angle) (millibarns per unit solid angle) 
‘Target element o,(E>60) o,ff£ 260) o(£ >60) 24° ye 
Beryllium 55 11 66 eee 99 +19 
‘Carbon 56 Wy 68 Rae 72+15 
Aluminium 84 21 105 +61 
22 4 — 


The last two columns give the experimentally measured differential cross 
sections (Cassels et al. 1951a); the two columns refer to different methods of 
measurement and angles of observation of 24° and 5° respectively. ‘The calculated 
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values correspond to 0°. The measurements at 5° are much more accurate than 
those at 21°. In the latter there was some uncertainty as regards the absolute — 
cyclotron current due to multiple traversal effects. In particular, the cross 
section for carbon should be multiplied by a factor 1-5 approximately, giving a 
differential cross section of 70 millibarns per unit solid angle. 

The energy spectra o,(E), o(E) for single and double scattering in the forward 
direction have also been calculated for Be, C, Al. These energy spectra, together 


Neutron Spectrum (millibarns per Mev per unit solid angle) 


60 80 100 120 «+4140 4160 180 
Neutron Energy (MeV) 


Fig. 2. Energy spectra of neutrons ejected in forward direction from beryllium as a function of 
the neutron energy E:: a, single collision spectrum o,(£); 65, double collision spectrum o,(F); 
c, o(£)=0,(E)+0,(E); d, experimental curve of Cassels et al. 


with o(£) =0,(E)+0,(£), are plotted against E for beryllium in fig. 2. Curve d 
of fig. 2 gives the experimental curve, with ordinates adjusted so that the 
maxima of the two curves agree. Very similar curves were obtained for the 
other two elements, and we shall therefore only discuss this case. 

As is seen from the table, the calculated differential cross sections are in 
fairly good agreement with experiment. The cross section for triple collisions 
one would expect to be of the same order as o, though somewhat smaller. Cross 
sections for fourth order and higher collisions should become less and less 
important. We have also neglected transitions leading to excited nuclei, which 
should only matter for the high-energy end of the spectrum but hardly for the 
cross section. ‘Thus, allowing for these processes will increase the cross section 
o(E > 60) by not more than 20°%, so that the agreement persists. 

As regards the energy spectrum, the intensity falls off too rapidly with 
decreasing energy, compared withexperiment. At100 Mev, where double collisions 
are most effective, the theoretical result is too small by a factor 4. This will be 
somewhat improved, but clearly not enough, by allowing for triple and higher 
collisions which lead to neutrons of lower energy. 

The energy spectrum depends much more delicately on the nuclear model 
used than does the differential cross section. For the latter it suffices to take 
account of the screening effects etc. in an approximate way, otherwise using the 
free particle cross sections. The energy spectrum, however, depends very 
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sensitively on the way the incident particle and the nucleus can exchange energy. 
Thus on our crude model better agreement is hardly to be expected. 

But one can see how the model has to be changed. We considered a nucleus 
consisting of individual non-interacting particles. 'To obtain a broader energy- 
spectrum one requires strong interactions between the nucleons, as in a liquid 
drop, for example. During a collision energy can then be exchanged in a large 
number of ways, so that the nucleus can remain in a much more highly excited 
state than our assumptions would allow. That an independent nucleon model, 
for instance the model we have used, or a shell model, is inadequate for treating 
nuclear reactions, and that for these one must assume strong interactions between 
nucleons, has recently been pointed out by Weisskopf (1950). 
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ABSTRACT. The frequency spectrum of a two-dimensional ionic lattice has been 
evaluated. A method is given for calculating the exact density of the spectrum at special 
points in phase space. ‘T'wo logarithmic infinities, first found by Montroll using a quasi- 
elastic forces model, are confirmed for an ionic lattice. The bearing of these infinities on 
the interpretation of the second-order Raman spectra is discussed. 


Spe LN TROD UCG TLOWN 
KNOWLEDGE of the frequency spectrum is necessary for the calculation 
of those properties of a crystal which depend on the vibration of its © 
constituent particles. A review of the investigations into the form of 
the spectrum has been given by Blackman (1941). 
Montroll (1947) deduced an exact expression for the frequency spectrum of 
a two-dimensional square lattice; quasi-elastic forces were used, adjusted to give 


* This paper forms part of the thesis submitted to the University of London for the degree 
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elastic isotropy. The exact analytic form of the spectrum had previously been 
unknown except for a one-dimensional lattice. Montroll found two infinities — 
in the frequency spectrum which had appeared as finite peaks in earlier work 
because of the use of numerical methods. These infinities were further shown 
to arise generally with any type of quasi-elastic binding. 

In this paper Montroll’s work has been extended to a two-dimensional ionic 
lattice, similar in structure to the (100) face of a rock-salt crystal. ‘The ions 
introduce long-range Coulomb forces. No analytic expression for the whole 
spectrum has been found, but a method for deducing the exact density of vibra- 
tions at special points is given. This method is valid in general for any type of 
force as well as Coulomb-law forces. The two infinite densities found by 
Montroll are confirmed for the spectrum of a square ionic lattice. Some new 
features arising in the phase-space contour maps and the frequency spectrum have 
also been described. 


§2. THE CALCULATION OF THE PHASE-SPACE MAPS AND 
THE FREQUENCY SPECTRUM 

The methods given by Blackman (1935 a, b) and Kellermann (1940) are used 
in this paper. The secular determinant is set up and a suitable set of frequencies 
is evaluated. ‘These frequencies are then used to construct phase-space maps 
from which the frequency spectrum can be deduced. 

These calculations have been performed for an ionic lattice of the type 
described in the introduction. The masses of the particles in the lattice are taken 
to be equal. ‘This model would appear to require two particles per cell, i.e. the 
negative and positive ions; but if we assume the repulsive forces between any 
pair of particles to be independent of the type of particle, then the equations of 
motion are the same for each ion and only one particle per cell need be used. 

We can write the vibrational equations of motion in a form similar to that 
given by Born (1923): 

(w? M+ Agata A yyy =0, 
A ,yu,t+(wM+ A,,)u,=0. 
In the above equation 
Aig= ‘Do OR™ expt =r sh eee (2) 


MMe 

The unit vector s specifies the direction of propagation of a plane wave in the 
lattice; +=27/A, where A is the wavelength; w=2z7v, where v is the frequency ; 
U,, U, are the components of the amplitude of an ion in the reference cell; 
lee a is the distance of a particle from the reference cell. For a square lattice 
yee oper Ying where x,,, = myo and Ying = Mas To is the distance between 
adjacent positive and negative ions and is taken to be 3A in the following 
calculations.-The m,, m, values are integral and range from — oo to +0. 

The potential energy between a particle at the point r,, ,,, and the particle in the 
reference cell is 0"; @7;""" is the appropriate derivate of the potential energy, 
which is assumed to be oh the usual form for alkali halides: O(r) = + e?/r + B/r?, 
where the first term is the Coulomb potential between two ions (the sign of the 
term varies according to the charges). ‘The second term is the repulsive potential 
energy; Bis a constant. Full details of the method for calculating the value of 
A,,, and similar terms by means of Epstein zeta-functions have been given in a 
paper by Smollett and Blackman (1951). 
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The secular determinant for the frequencies follows immediately from 
equation (1) and is 


47 Mv? + A, ae ayy 

vine 4a Mv" + A,, 

If we write 47° Mv? =f? we have at once the quadratic equation 
P= —(Avet Ayy) + (Ase Ay) +4A ay Fl. 0... (3) 


Now the boundary conditions of a cyclic lattice restrict the components of 
the phase vector to a discrete set ranging from 0 to z (see Born 1923). The 
phase components for the wave are 779s, and z7s,, and the conditions become 
$1 =77S,=a,7/N and ¢,=7rys,=a,7/N, where a,, a, are integers ranging 
from —N to + N and N is the total number of particles along the positive axis. 
This discrete variation of the phase components gives the correct number of 
normal modes, i.e. 8? in the secular determinant. 

The A,, now become functions of ¢, and ¢,. Also, as we may regard f as a 
frequency measured in units of 1/27M4?, then by eqn. (3) the frequency 
is also a function of ¢,, ¢2. The variation of the frequency with the phases is 
plotted in figs. 1 (a) and 1(d) in the form of a contour map. The two diagrams 
correspond to the upper and lower frequency branches given by the two roots 
in eqn. (3). The phase has been 4livided into units of 7/10 and the frequency 
is given in units of 10%e/27M??. 


0 es 4 6 8 10° 
a 

Fig. 1(a). The upper-frequency contour map. Fig. 1(b). The lower-frequency contour map- 

The frequency values are marked on the contours in units of 10%e/27M"?; ¢, is in units of 7/10- 


The frequency spectrum has been constructed from the contour maps. 
The spectra of the low- and high-frequency branches are shown in fig. 2 together 
with the complete spectrum. The units of density are arbitrary; the units of 
frequency have been given. The position of the infinities is indicated by lines. 
A description of the features of these diagrams is given later. 


§3. THE DENSITY OF THE SPECTRUM AT SPECIAL POINTS 
AND MONTROLL’S WORK 


Montroll (1947) has discussed the nature of the frequency spectrum of a 
two-dimensional square lattice with the masses of the particle equal and quasi- 
elastic forces between them. Montroll obtained an explicit formula for the 
frequency spectrum in the special example where the binding constants gave 
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elastic isotropy. This formula involved the integration of the expression for the 
density of normal modes. If g(f) - the py of modes for a frequency f then 


(= ga 5y|| rd 


where the double integral is taken over all portions of the phase-space map in which 
the frequency is less than f. The result of the integration was expressed in terms 
involving complete elliptic integrals. 

The density function confirmed the general appearance of the spectrum found 
by Blackman (1935b); but it also showed the important feature that the peaks 
in the spectrum were actually logarithmic-type infinities. Montroll proved as 
well that for quasi-elastic forces in a square lattice two logarithmic infinities would 
arise generally in the spectrum, one infinity belonging to the low-frequency 
spectrum and the other to the high-frequency spectrum. Montroll deduced 
this last result by a method of approximation that would appear to be useful for 
quasi-elastic forces only. 

The main purpose of this paper is to extend Montroll’s work to Coulomb 
forces. The density of the spectrum has been found at a number of special 
points depending on the features of the contour. lines in phase-space. The 
method which is developed here can be used for any type of force provided that the 
first three derivatives of f? with respect to ¢, and ¢, are continuous and finite. 
The two logarithmic infinities found by Montroll have also been found for 
Coulomb forces. 

First let us consider this method at the “end points’ in phase-space, that is 
at the points (7, 7), (7,0) and (0,7). If the derivatives of f? are suitably behaved 
we may expand f? about the point (z, 0) with respect to ¢, and d-as follows: 


rome) o(B) 
“ste, 9 Ph). 07 GB), 


where x and y are measured in the direction of ¢,,¢, from the point (7, 0). 

Now for the Coulomb and repulsive forces which we are using 0f?/0¢,, 6f?/dds 
and 0?f?/0¢,0¢, vanish at the point (7,0). This result follows from symmetry 
relations and also the fact that f? is in general a maximum or a minimum at this 
point. 

If we put fz =A, $(0°f?/0¢,"), 9 =B, F(Of?/0de?), >= C, where B and C 
are assumed to be positive, then, neglecting higher order terms, we find 
for a contour of frequency f close to (7,0) that f?=4A+Bx*+Cy?, giving 
x= Bl f? — A) —Cy?]}/#, Now in the Res for the density, 


(N= Ze 5 | a db 


we may replace | | dd, dh. by {Hh dx dy = ibe F x dy, 


where y, is small and is the value of ¢, where the chosen contourc uts this axis. 
Using the approximation for x given above we have 


Maes Maric Ne ae 
Ae) eta. pial(f?— A) — Cy? PP dy 


] fi f=A ot? 1 f eC? 
ae (BC)? | ‘| C y| dy = Aq? (BOR aE 
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Then with this result we may find the exact value of g(f) at the point (7, 0) 
by letting x and y tend to zero. 
If «=0 and y,—0 then y,C¥?/(f?— A)#?-+1 and we have 


Cater BO en re (4) 


This result is true for a contour that tends to an elliptical shape like that shown 
in the (7,0) corner in fig. 1(a). If the contour is hyperbolic as in fig. 1(6), at 
the (7, 0) corner, then the equation for f? becomes f? = A — Bx? + Cy”. Proceeding 


as before, we have 
A—f? 1/2 
As + ae [ [a4 eS a dy. 


Integrating, we find 
©. ie A—f? 
s)=— Zareoye 9 \—G 


If f?— A then g(f)— o. 

The density in this example rises to a logarithmic infinity. This method is 
justified because the large contribution to the density near (7,0) allows us to 
approximate to the rest of the contour. 

The results in eqns. (4) and (5) are the same as those obtained by Montroll. 
The equations deduced here, however, can be applied to any type of force at a 
point in phase-space where the contours have the necessary elliptical or hyper- 
bolic form. 

Equations (4) and (5) have been applied to an ionic lattice. ‘The second 
derivatives at the corner points have been calculated and the density of modes 
deduced. 

For the points (0, z) and (z, 0) in the high-frequency map, fig. 1 (a), the density 
of modes is found to be finite, as might have been expected from the shape of 
the contours. The agreement between the calculated and measured values of 
the density is found to be close. ‘The calculated value from eqn. (4) is 5-9 and 
the value measured from the contour map is 5:8 (these units being arbitrary). 

For the low-frequency map in fig. 1(5) the contours are hyperbolic at the 
(0,7), (7,0) points, and from egn. (5) a logarithmically infinite density: is found. 

At the (z,7) point the second derivatives become infinite and an expansion 
is not possible. In this example recourse must be had to the measured density 
in fig. 2 which shows that at the frequency 5-45 (in the arbitrary units specified) 
the density is a constant and the spectrum has no irregularities. ‘The derivatives 
were investigated at the (7,7) point for a general potential law 1/r”. It was 
found that the second derivatives approach infinity as nm tends to 2. For all 
values of m less than 2 the derivatives are infinite and no expansion can be made. 
‘The density can be shown, by means of eqn. (5), to approach ylny+D as 
n->2 where y+ 0 as the (7, z) point is approached; Disaconstant. ‘The density 
thus becomes a constant and not an infinity when n—>2. 

The density was also investigated in the region of the high-frequency contour 
map, fig. 1(a), where the contour lines form closed curves. ‘The density where 
the curves converge to apoint is found to be a constant of value 9-6, in agreement 
with the measured value of the density. which is 9-7. No unusual features arise 
here. 

The one other region where the density can be gauged is the region lying 
between the contour lines 5-7 and 5-8 in the high-frequency spectrum. In this 
region the contours can be seen to be approaching a saddle-point whose exact 
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) - a constant terms (a/c (5) 
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position cannot easily be found directly. However, at a saddle-point we know 
that contour lines become hyperbolic. In accordance with our analysis in eqn. (5), 
we know that hyperbolic contours will give rise to a logarithmic infinity in the 
frequency spectrum as long as the second derivatives are finite. Examination 
of the second derivatives shows that these can only become infinite at the point 
(7,7), but not at any other point. So we may conclude that between the fre- 
quencies 5-7 and 5-8 the spectrum has a logarithmic infinity. This conclusion 
is substantiated by the pronounced maximum shown in the measured value of 
the density. 


Fig. 2. The vibrational spectrum of a two-dimensional ionic lattice. 


1—Upper-frequency branch. 2—Lower-frequency branch. 3—Resultant curve giving the 
complete spectrum. v is in units of 10'%e/27M1?; g(v) is in arbitrary units. 


§4. DETAILS OF THE PHASE-SPACE MAPS AND THE 
FREQUENCY SPECTRUM 

The contour maps for a lattice with Coulomb forces differ in appearance from 
those found by Montroll using quasi-elastic forces. ‘The main difference in the 
high-frequency map, fig. 1 (a), is the existence of a peak-like prominence in the 
three-dimensional space of ¢,, ¢, and f, rising to a maximum frequency of 6-45: 
and lying on the diagonal. ‘This frequency is the highest in the upper-frequency 
branch and is also the highest frequency of the lattice. In quasi-elastic force 
models the highest and lowest frequencies lie on a corner point in phase space. 
This peak is a general feature of the contour map, and it has been calculated that 
the position of the peak varies very little with the variation of the repulsive power 
law. ‘The peak is not found for quasi-elastic forces. 

The low-frequency contours also have their highest frequency at a point other 
than the phase-space corners. It occurs at the points (7, 0-67) and (0-67, 7). 
In the region of the (7,7) point the contours split to form saddle-back surfaces 
in the three-dimensional space described. This phenomenon does not occur 
for quasi-elastic forces. 

The complete spectrum in plotted in fig. 2. The positions of the two 
logarithmically infinite peaks are indicated by vertical lines. The chief difference 
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from Montroll’s work is a subsidiary peak at the frequency 5-25 in the low- 
frequency spectrum. ‘This peak is not infinite as there are no ‘saddle-back” 
contours in the region of the contour map giving rise to this maximum. 


$3. CONCLUSION 

The infinite peaks found in the frequency spectrum are important in the 
Raman effect in crystals. We may conclude from Blackman’s work on lattices 
(1935 b, 1937) that the frequency spectrum of three-dimensional crystals will 
be similar to the spectrum of two-dimensional crystals: in fact the spectra of 
ionic crystals calculated by Kellermann (1940) and Iona (1941) show a strong 
resemblance to the spectrum of the two-dimensional ionic lattice shown in fig. 2. 
In the spectrum of a three-dimensional ionic crystal two steep peaks are found 
which we may expect by analogy to correspond to infinite densities, though this 
has not been proved. These peaks will give rise to steep narrow maxima 
associated with a continuous background in the second-order Raman spectrum. 
This spectrum was first observed by Fermi and Rasetti (1931). Born and 
Bradburn (1947) gave a detailed theoretical account of these maxima using 
Kellermann’s data (1940) to construct the theoretical spectrum. 
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ABSTRACT. ‘The band spectrum of manganese deuteride has been photographed at a 
dispersion of 1:25A/mm. The heads of the (0,0) and (0,1) bands of the 7I[I--*2 system 
are at A5703-7 A and A6073°3 A. respectively. It has been possible to make a rotational 
analysis of the two bands. The constants of the final state have been derived from the 
analysis and the fine structure resulting from the interaction of spin and rotation has been 
studied. The 7II state is perturbed and the terms have been evaluated from the term 
expressions for the final state and the wave numbers of the P, Q and R branches. ‘The 
results throughout are in good agreement with the results for the corresponding bands in 
the spectrum of manganese hydride. 
* Miss Conway’s part of the work was completed at Dunsink Observatory. 
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; §1. INTRODUCTION ; 

HE band spectrum of manganese hydride, first investigated by Heimer (1936) 

and Gaydon and Pearse (1937), consists of anumber of very complex bands 

distributed through the visible spectrum. Three of these bands degraded 

to shorter wavelengths, with ‘ heads’ at 5683 A, 6209 A and 5206 A, are respectively 

the (0,0) and (0,1) and (1,0) bands of a single system the rotational structure of 

which has been studied by Nevin (1942, 1945) and Nevin and Doyle (1948), to be 

referred to as I, II and III respectively, who have shown that the bands involve a 

transition from a7II toa’ state*. A diagram showing the principal branches for 

this type of transition is giveninI. Each band contains seven P, seven Q and seven 

R branches, six °P branches and six SR branches, which with other satellite branches 
bring the total of branches observed in the (0,0) band to fifty. 

The spin fine structure of the 7& state is well resolved and is determined mainly 
by the interaction of the resultant spin S with the rotation of the molecule. ‘The 
7II state as shown by the lines missing at the origin of certain branches is normal, and 
the appearance at low values of K, of N, Oand S branches, shows that at low quantum 
numbers the state tends to case a. At high quantum numbers it tends to case b. 
In both the vibrational levels v’=0 and v’=1 the lower rotational levels are 
perturbed in a regular manner from their expected positions. ‘The effect of this 
perturbation on the appearance of the P, branch of the (0,0) band can be seen very 
clearly in the reproduction of this bandinI. The lines of the branch appear to be 
converging to a head at 5683A, but instead of turning back again to shorter 
wavelengths in the usual way, diverge again to longer wavelengths. 

Because of the very small number of bands of multiplicity higher than three 
recorded in the literature it appeared worth while to photograph the spectrum of 
manganese deuteride, especially as a rotational analysis of the *I1—’> transition 
should be possible. 


§2. EXPERIMENTAL 
The discharge tube used to produce the spectrum is a modification of the design 
of tube used by Pearse and Gaydon (1938) in their work on manganese hydride. It 
consists of a Pyrex tube 50 cm long, 3-5 cm internal diameter, with three side 
tubes (see figure). A brass cap A waxed over one end carries two water-cooled 


Pump E Deuterium 
| 
i A 
ag a B 
{ee Same ae D 
oe —  — ———— 


Water 


terminals sheathed along their lengths with close-fitting Pyrex tubes. A spiral B 
of 12-16 turns of tungsten wire of diameter 1 mm is wound round an alumina 
capillary 7-5 cm long, 4 mm internal diameter, with a wall thickness 1-5 mm. The 
ends of the spiral fit into holes drilled in the water-cooled terminals and are held in 


* In The Identificati "3 
is Pecoptived Ges of Molecular Spectra (Pearse and Gaydon 1941 and 1950) the transition 
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place by set screws. One end of the capillary fits snugly into a hole in the centre of 
a steel disc C 6 mm thick which fits fairly closely inside the Pyrex tube. 

The tungsten spiral forms one electrode of the discharge tube. The second 
electrode D is formed by a sheet of nickel 10 cm long bent into a cylinder sprung 
into the Pyrex tube. Connection to the cylinder is made by a tungsten wire E 
sealed into a capillary tube which sheaths it along its length. This capillary is 
waxed into the middle side tube of the discharge tube. The other side tubes. 
serve as inlets for the deuterium and outlet to the Hyvac pump used to evacuate the 
discharge tube. ‘The discharge tube is placed in a wooden box through which 
water flows, the open end projecting through a hole in the side of the box with the 
joint made watertight. This open end of the tube is closed by a window fastened 
on with sealing wax. Before sealing on the window, five or six small pieces of pure 
manganese metal of total mass 0-2 — 0-3 gramme are introduced into the capillary 
tube and uniformly spaced along the part of the tube surrounded by the tungsten 
spiral. 

When a discharge is passed, the steel disc, as in Pearse and Gaydon’s tube, 
constrains the discharge to pass through the capillary, which becomes very hot. 
Additional heat is supplied by passing a current up to 50 amp through the tungsten 
spiral by means of a 15 v transformer joined to the water-cooled terminals. This 
heating current can be varied by a resistance in the primary. The use of auxiliary 
heating forthe capillary is advantageous as it enables the temperature of the capillary 
to be varied within limits independently of the discharge current. The discharge 
current used in the experiments was usually about 1-5 amp, obtained from the 
secondary of a 3000 Vv transformer. 

At the beginning the tube was cleared by a discharge current gradually raised 
to 1-5 amp with the tube connected to the pump. After about twenty minutes’ 
running the pressure became too low for the transformer to maintain the discharge. 
Deuterium gas stored in a two-litre flask was allowed to enter the tube through a 
needle valve. By proper adjustment of the pressure the nickel electrode could be 
made red hot, and this helped further to outgas the tube. By raising the pressure 
with a current of 30-40 amp through the spiral and a discharge current of 1-:2— 
1:-4amp the manganese deuteride spectrum could be made to appear strongly. 
Under the best conditions the discharge converged from the nickel electrode to the 
mouth of the capillary in the form of acone. At too high a pressure the discharge 
tended to pass from the nickel electrode to the tungsten spiral between the edge of 
the steel disc and the wall of the Pyrex tube. During a run deuterium gas was 
allowed to flow slowly through the tube, the rate of outflow to the pump being re- 
duced to conserve gas and the rate of inflow adjusted to maintain the proper 
conditions for the appearance of the spectrum. 

Although a considerable number of photographs were obtained, the tube 
cannot be regarded as entirely satisfactory. Owing to the necessity of allowing 
for expansion due to heating, the steel disc cannot be made too close a fit in the 
Pyrex tube. The small gap between the wall and the disc was increased by 
deviations of the Pyrex tube from cylindrical form, and occasionally, while 
adjustments were being made to obtain the best conditions, the discharge would 
strike through this gap, leading to fracture of the tube in a few seconds. The 
average life of a number of tubes was thirty-five minutes, the individual lives 
ranging from two minutes to two and a half hours. In each case the life was 
terminated by fracture at the steel disc. Another source of trouble was failure of 
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the alumina capillaries by cracking or puncture. Quartz discharge tubes did not 
prove satisfactory. They were not sufficiently cylindrical, and two specimens we 
had constructed tapered outwards from the ends to the centre, so that a steel disc 
which would enter fitted too loosely in the central part. 

In the course of the experiments the region 3800 A to 6200 A was photographed 
inthesecond order of the21-ft. grating inthis laboratory at adispersion of 1-25A/mm. 
Ilford Special Rapid panchromatic plates were used for the range 6200-4600 A and 
Ilford Special Rapid plates for 4600-38004. ‘The panchromatic plates were 
hypersensitized by bathing for two minutes in a 0-5°% solution of triethanolamine 
and drying rapidly with a fan. This procedure increased the red sensitivity by a 
factor of two and the green sensitivity by a factor of onehalf. The 75cm long plate 
holder enabled a range of 900 A to be photographed in one exposure at a particular 
setting of the grating. The exposures varied between fourteen and forty-five 
minutes. 

Each of the two bands with which the present paper is concerned was measured 
on two plates, one of moderate intensity to give good resolution of close strong 
lines and one of high intensity to show weak lines. The plates were measured with 
respect to second and third order J.A.U. standard iron lines and additional iron 
lines, the wavelengths of which were taken from the M.1.T. tables. 

The intensities of the lines are eye estimates checked against a microphotometer 
trace and large-scale prints of the bands. The microphotometer trace has also been 
used to check visual estimates of the resolution of very close groups of lines. Both 
bands are very well resolved in the second order. In the (0,1) band three greatly 
over-exposed Mn lines, 46013-5, A6016-5, A 6021-8, fall on the early part of the Q 
branches and on the middle part of the P branches but, beyond masking a single P 
and single Q line in some components, do not interfere with the analysis. 
The deuterium oxide used to prepare the deuterium gas was of 99°, purity, but 
the intensity of the MnH bands is about 10% that of the MnD bands, presumably 
due to occluded hydrogen. In the circumstances the intensity of the MnH 
bands is satisfactorily low and, though increasing the labour of measuring the 
plates, did not add to the difficulty of the analysis. 


§3. ANALYSIS OF THE BANDS 

Inspection of a twenty-fold enlargement of the plates of the (0,0) band shows 
a head at 5703-74 (17527-5 cm“), heads at 5699-74 (17540cm~) and 5696-14 
(17551cm~')andacomplex group of lines at 5693-1 A (17 560 cm7), all of which were 
considered to be °P heads. In the region 17570cm™! to 17610cm™! the P, 
branch and its satellite *Q,, can be picked out on the prints. From 17600cm7 
to 18500 cm“ the band is very complex and without any apparent regularities. 
Because of the closer rotational structure which corresponds to the higher moment 
of inertia, the occurrence of groups of seven to fourteen lines which had been 
observed in MnH at a distance of some hundreds of wave numbers from the heads 
is not evident in MnD. 

Rough calculations based on the isotope effect showed that at high rotational 
quantum numbers the separation of the lines of the Q and R branches should be 
about two-thirds of the corresponding separation in MnH and that the second 
wave number differences of successive lines in a branch should be about 0:4 cm—. 

Corresponding R and P branches are related by the combination difference 


ALF ,"(K)= F(K+1)— F,"(K-1)=R(K-1)— P(K+1). 
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The spin fine structure of the 7X state causes A, F,"(K) to decrease slightly with 
increasing z from component to component, but the effect is small, so that one can 
write 

A, F"(K)=(4B,+6D,)(K+4)+(8D,—34H,)(K+4)?+12H,(K+4)...... (2) 


The known values of the constants B, D and H for MnH and the well-known 
relations connecting the values of corresponding constants for isotopic molecules 
enabled the values of A, F’(K) to be calculated for MnD. 

Starting about 18 400 cm a search was made by trial and error for series of 
strong lines of the expected separation, decreasing to lower wave number in 
arithmetical progression, and eventually twelve such series or branches were found. 
These twelve branches could be divided into two sets, one containing seven and the 
other five branches, tentatively assumed to be Q and R branches, because at a few 
hundred wave numbers from their beginning the separation of the lines of each of 
the branches in each set was practically identical, while the separation in the R set 
was two wave numbers greater than in the Q set. ‘The two remaining R branches 
were found subsequently. 

The R branches were tentatively numbered R,, Rg, R; etc. inorder ofincreasing 
frequency at high quantum numbers. Successive values of K starting with K =0 
were given to the R, branch, and eqn. (1) was used in conjunction with the 
calculated values of A, F”(K) to calculate the P, branch. The correct numbering 
was obtained when the wave numbers of the calculated branch agreed closely with 
the P, branch already picked out onthe enlargements. Once the correct numbering 
of the R, branch had been found, the numbering of the other R branches was evident, 
and the corresponding P branches were readily calculated and the lines picked 
out on the enlargements and in the wave number tables. 

The °Q, ;,; branches are the high-frequency satellites to the P; branches with 
the numbering of each line the same as that of the adjacent P, line. The °Q,, 
branch is particularly clear. It is first resolved from the P, branch for K=6, 
the separation being 0-18cm™, and can be traced to about K=25, where the 
separation is 0:-44cm™~1. For convenience the separation of adjacent lines in a 
P, and ?Q, ;,, branch is denoted by Af/’;.,(K) = F,"(K) — Fi,(K). To number 
the Q branches and obtain the °P branches the relation 

A, Fi,1(K)=Q({K—1)—°P,(K+1)+4fii(K-1)  ....- (3) 
was used with an approximate allowance for Afj;,,(K—1). The procedure was 
similar to that used to correlate the R and P branches. The °P,, branch forming 


a head at 17526-5 cm“ represents the start of the band and is very clear. ‘The 
SR branches were found by the relation 
A, F;(K) = *Rit,i(K— YO, +1) TAG (Kl)... ssa (4) 

Substantial P,, and ‘P,, branches have been identified. All of the NP,, and a 
considerable part of the NP,, branch lies beyond the °P,; head ina region free from 
overlapping lines of other branches. ‘The lines of the branches can be calculated 
from theexpressions ‘P, ;,(K) ="Q),:41(K —2)— A, Fi',o(K — 1) + Afis, i42(K — 2), 
the identification of the lines in the calculated positions giving a valuable check on 
the correctness of the analysis of the components 7II]_,~’= and “II_,;~>*% at low 


values of K. 
The study of the (0,1) band was undertaken when the analysis of the (0,0) band 
was fairly complete, and the details do not requirecomment. The °P,, branch has 
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a head at 16461 cm™ (6073-34). From about 16500cm~ the lines are closely 
spaced, and this dense region continues for some hundreds of wave numbers. The 
Qand R branches can be traced right back to the start of the (0,0) band. 

When the analysis of the two bands had been carried as far as possible, the 
combination differences were tabulated in the form used for MnH in I. ‘To save 
printing, the tables of the wave numbers of the branches and combination 
differences have not been reproduced in this paper* as they can be reconstructed 
from the wave-number tables. In the final checking of the analysis and in tracing 
the branches to the lowest values of K the relations in eqns. (3) and (4) are very 
useful. Because of the nature of the fine structure, Af’, ,,(K) is nearly independent 
ofi,anditis correct towrite®R,,, (K —1)—Q,,,(K+1)=Q,(K—1)—°P, ,,,(K+1). 

These expressions furnish twelve sets of differences very nearly equal to one 
another in each band and slightly smaller than the true A, F(K) values. Since 


A, Fi,(K) = F,,(K +1) — Fi,(K— 1) =*R,,-1(K—1)— Q({K—-1) + Af_, {K—1) 
=Q(K+1)—°P, (K+ 1) + Mffin(K+D), 


the values of °R,;_,(K—1)—Q,(K-—1) and Q(K+1)—°P, ,,,(K+1) should be 
closely the same for a particular component in both bands as they have a common 
initial state. 

The various tables were systematically examined for disagreements and 
discrepancies. Any discrepancy exceeding 0-15 cm™ was specially studied. In 
some cases one of the lines involved was a blend or was diffuse because of imperfect 
resolution. Insome cases the wrong line had been chosen. ‘The structure of the 
bandsis much closerthanin MnH and thebranches go to higher rotational quantum 
numbers. Moreover, though the MnH bands appear weakly on the plates, a 
number of lines, because of their great intensity in MnH, are moderately strong. 
These factors mean that there is more blending and overlapping of lines, and 
although the inherent accuracy of measurement is the same as on the MnH plates, 
the actual wave numbers are not quite so accurate in parts of the bands, especially 
where the density of the lines is very great. The large number of lines on the 
plates makes it difficult in some cases to determine the exact end-point of 
the branches, a line of low intensity being followed by a strong line which by chance 
coincides with the next line of the branch in question. For the same reason it is 
difficult to be certain of the ®Q branches which lie in a crowded region of the band. 
A number of intensity irregularities may be attributed to the same cause. 

The analyses of both bands are internally consistent and consistent with each 
other and with the MnH bands.. It has not been found possible to trace the 
branches of the inner components to the values of K at which for a normal 7II state 
they should begin. At low values of K the lines of the branches become weak and 
it is difficult, if not impossible, to be sure of the correct assignments. In addition, 
it is possible that because of systematic perturbations the lines deviate more and 
more from their expected positions. It is possible to make assignments of lines 
but not possible to exclude the possibility that agreements in the combination 
differences are due to chance. It does seem, however, that a number of fairly 
strong °P lines belonging to the inner components are present on the plates of lower 
K values than are givenin the tables, but in the absence of combination differences 
correct identifications cannot be made. 


* A limited number of copies of the wave-number tables are available upon application to 
the authors or to the offices of the Physical Society. 
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§4. DISCUSSION OF THE INITIAL AND FINAL STATES 
The spin fine structure of 7% states is discussed in II. In MnH effects due to 
spin-spin interactions are negligible and the observed structure is determined 
almost completely by the interaction of the resultant spin S with K in accordance 
with the quantum-mechanical cosine rule 


T(K, J—- BK) =ty [JU +1)—- K(K + 1)—S(S +1), 


where J =K+3forf,...J=K-—3forf,. The separation of adjacent components 
of a rotational level is given by : 

Oia) =(4—2)y tyYK ae (5) 
and increases linearly with K. Consideration of the experimental data for MnD 
suggests that, because of blending, the most accurate data on the fine structure 
should be obtained from Af, , ,,(K) ="Q, ;.,(K)— P,(K), ignoring the results given 
by the other possible expressions. 

The observations in the (0,0) band are fairly complete and, within the limits of 
observational error, the slopes of the straight lines obtained by plotting Af, 4 
against K are identical and yield a value 0-0141 cm™ for y. As in MnH, the 
expression for the constant term in equation (5) is not completely confirmed by 
experiment. ‘The disagreement is undoubtedly due to the spin-spin interaction. 
The data for the (0,1) band are much less accurate and extensive but yield a 
value 0-0145 cm“! for y. 

It was shown in II that A,F7(K+4) is given accurately by eqn. (2), and 
this relation has been used to calculate By D; and H?, by the method already used 
for MnH. The final values of B; and Dj were calculated by the method of 
least squares. 

The agreement of the term values derived from 


F,!(K) = BY K(K + 1)+ D) K*(K +1)? + A, K(K-+ 1)?+47[2K(4—1) + (¢—8)(¢—1)] 
has been compared with the experimental values derived from A, F’,"(K) for Fi(K), 
Fi(K) and F7(K). While in MnH very good agreement was obtained for F'j(K), 
the difference between the experimental and the calculated values for F{(K) and 
F7(K) increased in proportion to K? to about 0-5 cm™ at K =35, the difference 
being positive for F, and negative for F,. ‘This discrepancy was attributed to a 
dependence of the fine-structure constant y on K and could be removed by assuming 
y=y(1+pK). In MnD satisfactory agreement between experimental and 
calculated term values is obtained throughout for a constant value of y, and if any 
dependence of y on K exists it is too small to be detected. 

Using the ¢alculated term values of the 7X state, the terms of the “II state 
can be calculated from the expressions 

F,,(K) =R{K—-1)+ F,'(K-1)=P(K+1)+ £,"(K+1) 

and F,(K) =Q{K) — Fi"(K). 
These terms are given in table 1. The A doubling for the seven components 
defined by Av,,, = F,(K) — F,,(K) can be obtained at once from the term values. 

For multiple states intermediate between case a and case b 


t=7 
LY AR) =(4B, £6D,\I +4) +8D (J +9) + 12H (J +h. ...... (6) 
=i! 


To reduce the effect of the A doubling A, F,(/) is taken as the mean of A, F,,(/) and 
A,F,,(J). In spite of the perturbation of the ‘II state, the mean of the seven values 
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of A, FJ) from J =14, the lowest value observed for all components, to J =45, 
the highest value observed, agrees remarkably well with eqn. (5). However, in 
view of the perturbations, it is not certain that the constants have their usual 
meaning, and their values have been determined graphically. 

In the absence of theoretical expressions with which to compare the terms it is. 
difficult to discuss adequately the overall distortion of the “II state especially as it 
has not been found possible to trace the branches to such low values of K as in 
MnH. A comparison of table 1 with the corresponding table 2 for MnH in II 
indicates that the general character of the perturbations is similar in the two 
cases. In particular a distortion of the P, branch occurs similar to that in MnH, 
but not of so marked a character. The lines are closest together about K=7 
but, instead of forming a head, diverge again to longer wavelengths. 


Table 2. Molecular Constants derived from the (0,0) and (0,1) Bands 
of the 7I1—7X Transition in MnD 


oa —0 C1 

B’, (cm~) 2°86959 2°81427 

Dj (cm?) —7:9514 x 10-5 —8:2212 x10 

ELS(Cinve:) 1-41 x 10-° 217A Gee 

2y (cm?) 0-0141 0:0145 ; 
Be 89 72 cies Di=— 78165: 10" cms oe =0:05532 cm=}, 
Bo=0°2698 x 10-4 cm=!, He=AO3p-a1 Ome cress Ve=0-70 <10-" chia 
Ig=9-6598 x 10-*° g.cm?, re=1-7234 X 10=2cm: 


Initial state v'=0 : 
By’ =3-244 cm DD, == 9-51 x10 car | a ee 

The constants of the initial and final states derived from the rotational analysis. 
are collected in table 2. It is well known (Herzberg 1950) that for deuterides and 
hydrides the spectroscopically determined value of p”, the ratio of the reduced mass 
of XH to XD, is greater than the value determined from the atomic masses. 

For MnH and MnD, p?=8B,MnD)/B,(MnH) =0-50959, while the value 
determined from the atomic masses is 0:50932. 
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ABSTRACT. An apparent selective excitation to the v’=1 level of the CH 39004 
(x11) system occurs in hydrogen flames containing small amounts of organic compounds 
since the (1, 0) and (1, 1) bands are abnormally strong. This is accompanied by a suppres- 
sion of the usual predissociation in this system. An additional band of this system, the 
(1, 2) at 4494-5 A has been observed for the first time. 

The suppression of the CH predissociation is interpreted as showing that conditions 
in hydrogen flames approach thermal equilibrium. The appearance of the CH predis- 
sociation in normal flames indicates that the excitation is not thermal. 


§1. INTRODUCTION 
LTHOUGH CH, together with C, and OH, bands form the strongest 
AN feature of the spectra of all hydrocarbon—oxygen (and air) flames the 
mechanism of the formation and excitation of the CH and C, radicals 
is unknown (Gaydon 1950). 

Under normal flame conditions where CH is excited, only the (0, 0) bands of 
two different systems appear. The 4315 A (0,0) CH band 1s always the strongest 
and arises from a 2A->*I] transition. The 3872 A (0,0) CH, due to a 22>? 
transition, is also reasonably strong. In very hot flames (1.e. oxyacetylene) the 
3628 A (1,0) and the 4025 4 (1,1) bands of the CH 3900A ?X —?II) system 
(which has (0,0) head at 3872 A) also appear; under these conditions the 3143 A 
CH (?2+-?IT) band may also appear weakly (Gaydon 1948, p. 44). 

During an investigation of the spectra of hydrogen—air flames containing small 
amounts of the fluorocarbons, the C, and CH (and other) bands were observed 
but there appeared to be an apparent selective excitation to the v’ =1 level of 
the 3900 A CH system since the (1,0) and (1,1) bands were unusually strong 
(Plate, (a)). The same effect was also observed when small amounts of carbon 
tetrachloride, methyl iodide or ethylene were added to a hydrogen—air flame. 
The phenomenon thus appeared to be a characteristic of the hydrogen flame 
rather than of the organic additives. 

This effect has been subjected to a. detailed investigation and the results 
together with a possible explanation are given in this paper. 


§2. EXPERIMENTAL 

The change in the spectra of flames of both ethylene and acetylene with both 
air and oxygen has been studied as the hydrocarbon is gradually replaced by 
hydrogen. During this replacement the total mixture strength, defined by A, 
was kept constant. A is the ratio of the amount of oxygen actually used to that 
required for complete combustion of the total fuel (i.e. hydrocarbon + hydrogen) ; 
A=1 for a stoichiometric mixture, A<1 for a fuel-rich mixture and A>1 for a weak 
mixture. ‘These experiments were carried out for a wide range of A. 
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The spectra of hydrogen—oxygen flames containing small amounts of methane, 
hexane, benzene, methyl alcohol, acetone and acetic acid have also been studied. 

All flows were measured using capillary type flowmeters. A quartz Smithell’s 
flame separator was used for the air flames and an oxyacetylene welding burner 
for the oxygen flames. 

The spectra were photographed on a Hilger E518 Raman spectrograph 
(aperture f/4) using both the quartz and glass trains. Large dispersion plates 
were obtained with a quartz Littrow spectrograph. 

To obtain an estimate of the relative intensity of the (0,0) and (1,1) bands 
the plates were calibrated using a tungsten strip lamp and rotating stepped sector. 


§3. RESULTS AND DISCUSSION 
Vibrational Intensity Distribution in the CH 3900 A System 


For the hydrogen-ethylene—air flame the high excitation of the (1,0) and 
(1,1) CH bands only occurs for low ethylene concentrations in rich mixtures 
(i.e. A low). The effect does not set in suddenly, but gradually appears as A 
and the ethylene concentration are reduced, becoming noticeable with 3% 
ethylene (expressed as a percentage of the total fuel, i.e. hydrocarbon plus 
hydrogen) and A=0-4; it is strongest at approximately 1-5°% ethylene at A=0-4. 
It is observed for all ethylene concentrations less than 3% when A is less than 0-4 
until eventually the ethylene concentration is too small to give an inner cone and 
therefore no C, or CH appear at all. 

In the hydrogen—ethylene-oxygen flame the effect is apparent when A is 
less than 0-4 and ethylene concentration not greater than 2°. Similar results 
have been’ obtained in the hydrogen—oxygen flames when acetylene, methane, 
benzene, acetone or acetic acid are used. 

Under conditions where the (1,1) and (1,0) CH bands are strongest an 
additional red-degraded band appears with head at 4495-5 A (see Plate, (d)). 
This band appears to be the R head of the previously unknown (1, 2) band of the 
3900 A CH system since it gives a value of x,”w,” =62:5 and w,” =2844 in good 
agreement with the values already known for the CH ground state (Herzberg 
1950, p.518). The appearance of this band provides further evidence of the 
high excitation to the v’ =1 level. 

An estimate of the relative intensity of the (0,0) and (1, 1) band in a hydrogen— 
ethylene—air flame (2°% C,H,, ’=0-4) gives the intensity of the (1,1) band as 
being approximately one quarter the intensity of the (0,0) band. 


Suppression of the CH Predissoctation 


When flames showing the above effect are studied under higher dispersion 
(i.e. with Raman glass spectrograph) the high excitation to the v’ = 1 level is found 
to be accompanied by an unusual appearance of the rotational structure of the 
CH 3872 & (0,0) band (see Plate, (d)). 

In the usual discharge and arc sources and in normal flames the 3900 A CH 
system exhibits a well-marked predissociation (Gaydon 1948) in the (0,0), (1,0) 
and (1,1) bands. Each branch in the rotational structure remains normal in 
intensity until a certain K’ value beyond which the next line becomes diffuse and 
faint and the branch suddenly stops. Under these conditions the last observed 
line is usually K’=17 for (0,0) and K’=7 for the (1,0) and (1, 1) bands (Shidei 
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1936). However, in the hydrogen flames studied the predissociation does not 
set in and the branches extend to higher K’ values. Lines up to K’ =20 in the 
(0,0) band and K’=13 in the (1,0) and (1, 1) bands appear (Plate, (c)). 


Interpretation 


Predissociation may be observed in emission when the excited state of a 
molecule is being depopulated by radiationless transitions, either to an unstable 
state whose potential energy curve intersects the potential energy curve of the 
stable excited state, or, as appears to be the case for the ?&~ state of CH, by a 
quantum-mechanical passage through a potential barrier (i.e. tunnel effect). 
All lines in a band arising from levels above the point where these radiationless 
transitions begin will therefore tend to be considerably weaker or even entirely 
absent (for full treatment of the theory of predissociation see Herzberg 1950, 
p. 405). 

When thermal equilibrium exists the number of molecules predissociating 
in the upper state will be balanced by the formation of new molecules in this state 
by the inverse process and the intensity of the spectral lines in emission will be 
unaffected although the levels affected by predissociation will be broadened by 
the lack of rigid quantization. Therefore the weakening of lines arising from 
the predissociated levels immediately shows that the excitation of the spectrum 
is not of a thermal nature. 

Weare therefore led to conclude that conditions in very rich hydrogen flames 
(where CH predissociation is suppressed) more nearly approach thermal equili- 
brium and the CH is excited thermally when small amounts of organic compounds 
are introduced. The unusually high intensity of the (1,0), (1,1) bands and the 
appearance of the (1,2) band occur because the v’=1 level is now populated 
according to the Boltzmann distribution law. However, since the vibrational 
transition probabilities for CH are not known, it is not possible to determine the 
relative intensity of the CH bands to be expected at the flame temperatures 
(assuming thermal equilibrium) for comparison with the observed relative 
intensities. 

This conclusion is in complete agreement with the observation that the 
radiation in the (0, 0) OH band is probably thermal in origin in the oxy-hydrogen 
flame (Gaydon 1948, p.30). Further, the rotational temperature of the OH 
(0,0) band in low pressure oxy-hydrogen flames is around 2300°k compared 
with a theoretical flame temperature of 2500°K (Gaydon and Wolfhard 1948). 


The Excitation of CH in other Flames 


The observation of the CH predissociation in normal hydrocarbon—air and 
hydrocarbon—-oxygen flames shows that the CH is not excited thermally. ‘This 
supports the observations of Gaydon and Wolfhard (1949) that in low pressure 
flames the CH has an abnormally high translational temperature; this has led 
these authors to favour excitation by collision with some active species in the 
flame (rather than a true chemiluminescence since the rotational temperature 
is approximately normal). 

The gradual setting in of the suppression of the results of predissociation 
when the hydrocarbon is replaced by hydrogen suggests that at some stage the 
CH excitation is occurring by at least two competing processes, i.e. thermal 
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and the usual hydrocarbon mechanism. The thermal excitation gradually 
becomes predominant at low hydrocarbon concentrations and for low values. 

Gaydon and Wolfhard (1951) have shown that the reaction zones of most 
organic flames exhibit abnormally high electronic excitation temperatures since 
when iron carbonyl is introduced lines in the Fe spectrum requiring at least 173 
kcal/mole for their excitation appear. ‘This abnormal excitation does not extend 
into the interconal gases, or the mantle that is often present around the inner 
cone in rich flames. This caused them to suggest that these mantles are probably 
inthermal equilibrium. However, astudy of the CH bands emitted by the mantles 
in oxyacetylene and oxyethylene flames shows that the predissociation is observed 
as in the reaction zone. This shows that thermal equilibrium has not been fully 
established. Further evidence for this is provided by Lewis and von Elbe 
(1943) who have observed that temperatures measured by the sodium-line 
reversal method exceed the theoretical temperatures for this region of the flame 
for rich mixtures of fuel with oxygen. 


§4. CONCLUSION 


On the basis of this work it seems that observations on the predissociation 
in CH when organic substances are introduced into flames may provide a simple 
method of determining if thermal equilibrium has been established. 
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ABSTRACT. Powder patterns have been obtained on the (100) and (011) surfaces of a 
carefully shaped single crystal of 3° silicon-iron, of dimensions 13 mm x 6 mmx 6 mm, 
and these have permitted a more thorough quantitative test of proposed ferromagnetic 
domain structures than has been made hitherto. The patterns on the (011) surface are very 
beautiful and complicated. Although the measurements of domain spacings on both planes 
under Mode III magnetization conditions, as well as measurements of the angles between the 
domain vectors and main domain boundaries, confirm the essential features of the main 
structure proposed by Néel, there are found considerable numerical discrepancies which are 
attributed to the complicated nature of the closure domains on the (011) surface. 


§1. INTRODUCTION 
N this communication we give an account of a careful examination of the 
] domain structure of a single crystal of silicon—iron by the powder pattern 

technique. In essentials, this technique is a kind of field plotting in which 
the coarse iron filings of one’s student days are replaced by the fine particles in a 
colloidal solution of magnetite. When a drop of this solution is placed upon the 
carefully polished surface of a single crystal of a ferromagnetic metal, in general 
the colloidal particles settle upon the surface to form patterns which are 
sometimes called Bitter figures (Bitter 1931) and which are easily photographed 
using a metallurgical microscope with bright field illumination. In the papers 
by Williams, Bozorth and Shockley (1949) and by Bates and Neale (1950) full 
accounts of the method of preparation of the solution and necessary photographic 
procedure are given. Consequently, we deal here only with certain improvements 
which have been effected in the course of the present work. 

The technique has provided much new information concerning the domain 
structure of ferromagnetic crystals. Very beautiful figures and convincing 
interpretations of certain characteristic patterns, such as ‘fir-tree’ patterns, 
have been given by Bozorth and his colleagues, while attempts at quantitative 
measurements have been made by Williams and Shockley (1949) and by Bates 
and Neale. In particular, the latter workers studied the powder patterns on the 
surface of a silicon-iron single crystal strip which was (approximately) a (100) 
plane, in which the [110] axis was (approximately) parallel to the long axis of 
the strip. This closely approaches the case which was studied theoretically by 
Néel (1944), and for which he suggested the domain structure shown in fig. 1 (a) 
when no field, or only a very weak one, is applied along the [011] direction. 

In fig. 1(a) the main domains are leaf-shaped sections with their planes 
perpendicular to the [011] axis and their magnetization vectors alternately 
directed along the [001] and [010] directions. In addition, there are postulated 
certain closure, or antiparallel, q domains, in order to reduce the very high 

* Now with The Steel Company of Wales Ltd., Lysaght Works, Newport. 
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magnetostatic energy which would otherwise arise from the presence of free 
poles distributed along the sides of the specimen. oy rae 

Whenasmall (effective) field acts upon the crystal parallel to the [011] direction, 
we have the onset of what Néel terms Mode III magnetization (see Bates 1950, 
1951b) and the domain structure suggested is that of fig. 1(6). The 
magnetization vectors of the main domains now turn inwards to make an angle @ 
which is a function of the effective field H with the [011] axis. At the same time 
there appears a new set of closure domains, the so-called parallel, p domains, 
which are magnetized as shown in fig. 1 (8). 

It is the latter case which gives the possibility of a quantitative check of the 
domain theory, since Néel proved that for a single crystal of pure iron the 
domain period represented by the distance D, which is the thickness of a pair of 
successive domains, is given by D=[yL(W,+W,)/W,W,})”, where y is the 
energy per cm® of a domain wall or boundary between successive domains, 


4 
(lod) | (@) No field applied 


(b) Field applied 


Fig. 1. Domain structures proposed by Néel. 


r 


L is the width of the strip specimen, and W,, and W,, are respectively proportional 
to the energies in erg/cm* required to change the directions of magnetization of 
the main domains to those of the p and q domains. W,, and W, are therefore 
functions of 8, which in turn is a function of H. It is therefore possible to plot 
a graph of the theoretical values of D against H. 

Bates and Neale (1949, 1950) made a large number of measurements of 
D and H, and their work undoubtedly provided much evidence for the essential 
correctness of Néel’s views, but it was incomplete in certain respects. For 
example, they had to use thin strips of single crystal specimens of silicon—iron 
made from transformer sheet material, whereas the Néel theory was given for 
the case of a thick crystal of pure iron. Consequently, when it was found that a 
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predicted minimum in Néel’s (D, H) curve was not found in the experiments,. 
it was natural to conclude that the excessive thinness of the specimens was the 
cause of the lack of agreement. ‘The magnetic anisotropy coefficients of pure iron 
and of silicon-iron are not sufficiently different to be regarded as a main cause 
of the lack of agreement, particularly as the anisotropy constant enters into the: 
expression for D as K-"*. But, of course, the excessive thinness of the material 
meant that it was impossible to shape the specimen so that the surface was 
exactly a (100) plane. Moreover, there was no appreciable [110] face on the 
specimens, whose edges were rounded by mechanical and electrolytic polishing, 
thus making the closure domain surfaces semicircular instead of plane, so that 
additional magnetostatic energy could arise from free poles on them. In any 
case, at low fields the width of a q domain is of the same order as its depth, hence 
for free poles to be maintained on its ends much energy would be required. 
Indeed, it almost looks a priori as if the periodicity of the whole primary domain 
structure might be profoundly affected by accidental blemishes on the edges of 
the specimen. It was therefore clear that to obtain more reliable quantitative 
results thick crystals had to be used and their surfaces more ideally shaped than. 
heretofore. 
§2. EXPERIMENTAL TECHNIQUES 


An excellent ingot of 3° silicon—iron some 7 cm long and about 1 cm in 
diameter was kindly presented to us by Dr. D. Shoenberg. It contained about 
six large single crystals and the largest of them was cut from the ingot by means: 
of a rock cutting machine, a rapidly rotating diamond-impregnated disc which 
could be pressed very gently against the metal to give minimum strain, the 
metal being water cooled during cutting. 

After etching to reveal [100] planes, the approximate orientation of the 
crystal with respect to the specimen faces was measured with a two-circle 
goniometer as in the work of Bates and Neale. The surfaces were then ground 
parallel to the desired crystal planes to within some 3°. This was done by 
mounting the crystal in a brass face-plate with sealing wax. he orientation of 
the crystal with respect to the mount was calculated from the goniometer readings: 
and the required crystal plane was then set parallel to the face-plate, by placing 
the face to be ground on another face-plate whose setting with respect to the 
first face-plate could be adjusted with a micrometer. ‘The crystal was then 
ground parallel to the first face-plate. When four sides had been ground in this. 
way, the specimen had as its longest direction the [011] direction, and as sides. 
the (100) and (011) planes, and the following treatment permitted us to cut the 
crystal surfaces parallel to crystal planes to within some half a degree. 

A surface grinder was modified to grind off layers of 0-0005 in. maximum 
thickness per cut, in order to minimize strains in the specimen. ‘The same 
replaceable crystal mount was used in the grinder and in an x-ray goniometer. 
After a flat surface had been ground on the crystal, the deformed surface layer 
was removed by a rapid electropolish, and the two x-ray photographs, necessary 
for computing the angles between the surface and the crystal plane measured 
about two axes of rotation on the surface grinder, were taken. ‘The specimen 
was then replaced in the surface grinder and rotated through the known angles 
to bring the crystal plane parallel to the grinding surface. ‘he new surface was 
then ground parallel to the desired crystal plane, and its orientation checked by 
another x-ray examination. This procedure was also carried out for the 
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(100) face and the (011) face; the other two sides of the crystal were then ground 
parallel to these planes, and the ends of the crystal perpendicular thereto. The 
final dimensions were 13 mm x 6 mm x 6 mm, approximately. 

Following the final grinding, the cold-worked surface metal was removed by 
electrolytic polishing in which great care was taken to avoid rounding the crystal 
edges by mounting the specimen in a polystyrene mould. It was found necessary 
to use a special type of polishing cell with a rotating anode and an electrolyte 
with a slower polishing action than that used by Bates and Neale. (In later forms 
of the apparatus used at Nottingham both anode and cathode are rotated.) 
A special bridge (Bates and Mee 1950) was devised for the automatic measurement 
of the resistance of the cell during polishing. This was particularly important 
as we obtained best surfaces with the non-explosive chromic acid-acetic acid 
bath (800 cm® glacial acetic acid, 38 cm® distilled water and 150g chromium 
trioxide), for which the conditions for maximum anode film thickness are rather 
critical. The usual polishing time was about 20 minutes. After polishing, the 
specimen was washed in commercially pure alcohol and dried in a hot air blast. 
The specimen was then annealed for one hour at 1 100°C im vacuo. 


§3. RESULTS 

Bozorth and his colleagues found that the patterns observed on a surface 
of a demagnetized silicon-iron crystal give conclusive information concerning 
the inclination of the surface to a true (100) crystal plane. We obtained many 
excellent ‘tree’ pattern photographs during the initial grinding of the specimen 
and from them we were able to make valuable deductions about the inclination 
of the surface and to satisfy ourselves of the correctness of Bozorth’s interpretation 
of such patterns. Considerations of space prevent their reproduction here, and 
although some of the patterns allowed us to deduce, for example, that in some 
cases of large inclination the mechanism of the change from Mode I to Mode III 
magnetization consisted of elongated ‘tree’ branches growing deep into the 
crystal by 90° boundary displacements and by the elimination of 180° domains 
by 180° boundary displacements, attention will be confined to a few of the 
patterns on the more perfectly ground surfaces only. 

In fig. 2 (a) (Plate I) are reproduced the patterns found near the edge of the 
(100) plane when an effective field of 8 Oe, sufficient to produce the two-phase 
structure of fig. 1(d), is applied-in the [011] direction. A marked decrease in 
domain width occurred with an increase in field of only 2 Oe; incidentally, the 
effective field was always measured directly with a small magnetic potentiometer. 
The reason for the heavy deposits of colloid was that the crystal surface was 
inclined at an angle of about 2° to one easy direction of magnetization, the 
[010] direction, and so free poles were formed. However, the [001] direction 
was in the surface and therefore closure domains magnetized along the latter 
direction are colloid-free and show up sharply against the black background. 
Alternate domains magnetized in the [001] direction have upon them streaks of 
colloid perpendicular to the direction of magnetization. The patterns are 
similar to the ‘tree’ patterns of the second type observed by Bozorth for small 
inclinations of the surface to the (100) plane. However, nearer the edge the 
surface inclination is greater and branches are formed at angles of 45° to the 
main domain boundaries. These branches, in turn, have small secondary 
branches, and it is possible to give a relatively simple diagram of this structure. 
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Similar patterns with 71° walls have been observed (Yamamoto and Iwata 1951, 
see also Bates and Wilson 1951) near the edge of a (110) plane of nickel where the 
main domain structure has 71° walls corresponding to the 90° walls of iron. 

When the surface approached still more closely, within 0-5°, to a (100) plane, 
the two-phase structure of fig. 2(a) had no closure domains. Fig. 2(b) shows 
the pattern on a specially chosen part of the final surface where the striations of 
colloid—perpendicular to the directions of magnetization—show up well, but 
no colloid has settled on the domain boundaries. We think that this is due to the 
effective field having a vertical component in opposition to the local field of the 
domain boundaries in this region, so preventing the attraction of colloid particles 
to them. 

In figs. 2(¢) and 2(d) are reproduced low magnification photographs of the 
line patterns observed on the (100) and (011) planes with moderate fields. We 
made many quantitative measurements with such photographs and the results 
are described below. Attention is now directed to the complicated structure of 
the individual lines of fig. 2(d). The spacing between successive lines on the 
(011) plane was identified with the period of the main structure; quantitative 
measurements are described later. In higher fields the lines under low 
magnification have the appearance of doublets, but under high magnification 
the doublets are seen to be a complicated secondary structure which grows in 
size as the magnetization is increased. Incidentally, ‘doublets’ were previously 
observed by Kaya (1934) on a (110) surface of iron magnetized along a 
[110] direction, and probably represented a similar structure. 

In fig. 3 (Plate II) is reproduced a set of high magnification pictures of the 
complicated structures for effective fields ranging from 13 to 110 Oe along the 
[011] direction. It was found that as the field was increased the complicated 
domains at first increased in width, and then multiplication of the pattern 
occurred by a new set of these domains suddenly shooting across the surface. 
‘To make room for this new set, the existing structures were reduced in width 
so that the new distribution was uniform. In fig. 3(d) is shown a wide pattern 
about to divide. 

We shall, provisionally, identify the black, heavily-coated deposits in fig. 3 (a) 
with the theoretical p closure domains. In fig. 3(a), Mode III magnetization is 
beginning, the ratio of the width of a p domain to the width of a q domain is very 
small, and the heavy deposits on the p domains manifest the presence of free 
poles. In fig. 3(5) we see that an increase in field from 13 to 17 Oe has caused 
a complete alteration in the appearance of the p domains, for a central, heavy 
deposit of colloid now separates two colloid-free portions within the p domain. 
The horizontal striations of colloid in the q domains show that the latter are 
magnetized along the [100] direction, which is parallel to the long edge of the 
paper. 

In several parts of these pictures one can see no visible boundaries between 
the q domains and the colloid-free portions of the p domains, e.g. in the region 
of fig. 3(c) slightly to the right and slightly below the centre of the picture, an 
observation which indicates that in such places the p and q structures are 
magnetized approximately 1 in the same direction. 

In fig. 3 (d) an increase in field has again caused the p structures to increase 
and again their appearance is greatly changed, for they are now further 
complicated by a fine lace-like structure; no significant change occurred on 
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reversing the field, and up to a field of 100 Oe the p structures multiplied and 
lost clarity, without further significant change in pattern. With higher fields the 
q domain boundaries are no longer parallel to the [100] direction but assume 
zigzag forms as in fig. 3(f), which are somewhat like the colloid-free areas of the 
p domains of fig. 3(d). Above 250 Oe the patterns became less clear and 
rather ‘woolly’, because the colloid particles were very strongly directed parallel 
to the applied field and were less affected by the localized domain fields. 

Figs. 3(e) and (f) were taken when the thickness of the crystal had been 
deliberately reduced to 0-13 cm by systematic cutting in the course of experiments 
described below. Fig. 3(e) is interesting because it shows a local disturbance 
which was probably due to a sub-surface inhomogeneity, and here again the 
disturbed q domain has the zigzag form similar to the colloid-free areas of the 
p domains in fig. 3 (f). 

It was frequently noted that the colloid deposits were much heavier on one 
side of the p domain structure than on the other. This was attributed to a small 
component of field perpendicular to the surface. This explanation was tested 
by deliberately applying a small field, perpendicular to the surface and 
insufficient in magnitude to alter the main domain structure. The results are 
shown in fig. 4 (Plate III), and are consistent with the suggested explanation. 
They are most helpful in determining the underlying structure of the p domains 
since the sign of the surface free poles can be inferred and any predicted domain 
structure must produce surface poles of the same sign. It must also explain the 
beautiful complementary detail of the patterns when either the horizontal or 
the vertical field is reversed. 


$4. DISCUSSION AND INTERPRETATION OF OBSERVED PATTERNS 

As found by Bates and Neale, and often confirmed in the present work, 
colloid deposits occur only on alternate boundaries on the (100) plane. Mee 
(1950) has given an explanation of this fact based on his observations of the 
patterns observed on a single crystal of cobalt. In his view, the direction of 
polarization of the colloid dipoles by the external field, as well as the direction 
of the surface fields attracting them, must be taken into account. Thus, in 
zero applied field, the particles are attracted to the highly localized stray fields 
at the boundary walls and near the surface they will be oriented almost 
perpendicular to it. On applying a field parallel to the surface, the particles are 
aligned parallel thereto and no longer attracted to the boundaries themselves. 
They are only attracted by the less highly localized surface fields due to free 
poles of opposite sign appearing at the intersection of adjacent domains with 
the surface. But for the (100) plane under investigation these surface fields 
are in opposition to the applied field at alternate boundaries and are therefore 
ineffective, particularly as the colloid particles are so polarized that these opposition 
fields repel them. In this instance it is not possible to make the favoured and 
unfavoured boundaries exchange places by reversing the field, since the whole 
domain structure is altered by so doing, whereas the cobalt structure studied 
by Mee (1950) is unaffected by reversal of a small applied field. 

The complex but regular patterns found on the (011) plane prove that the 
simple domain structures proposed by Néel, fig. 1, are not entirely correct. 
The fact that the domain spacing on the (011) plane is D and not 4D proves also 
that the domain scheme suggested by Lawton (1950) is unlikely. However, the 
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observed patterns have much in common with the Néel scheme. Firstly, the 
q domains have colloid striations running in the [011] direction (fig. 3 (b)), and 
are therefore magnetized in the [100] and [100] directions—a point which we 
proved by some compression experiments in which it was found that a 
compression in the [011] direction caused these domains to increase in size, 
because the magnetostriction coefficient of this silicon-iron is positive. 
Secondly, in the very weak fields required for the onset of Mode III 
magnetization, the q domains are very wide and cover the whole surface. 
Thirdly, with increase in field p domains appear between the gq domains, although 
the observed p structure is quite different from the simple form postulated by 
Néel. 

One naturally enquires whether further information on the closure structures 
can be obtained from observations of patterns on the (100) plane near the edge 
of the crystal. Unfortunately, the inevitable slight rounding of the edge caused 
by electrolytic polishing permits the formation of ‘tree’ patterns which mask 
the closure structure. 

In fig. 3(a) the new domains formed at 180° boundaries as the field is 
increased are heavily coated with colloid except for certain lozenge-shaped, 
colloid-free patches. These patches appear to represent a basic unit from which 
many of the more complicated structures are built (for example, they are found 
in the fine structure of figs. 3(d) and 3(e)). As magnetization proceeds, the 
long colloid-free lozenges divide into two sets with a heavy colloid deposit 
between them, as in figs. 3(b) and 3(c), and this new shape of colloid-free areas 
persists with only minor changes until magnetic saturation obtains. 

We think that these areas, while remaining magnetized parallel to the surface, 
have their magnetic vectors turned away somewhat from those of the adjacent 
q domains. For, if the regular pattern is disturbed by a crystal imperfection, 
the q domains are sometimes directly joined to the p structure, as in fig. 3(e), 
which means an approximate alignment of their magnetic vectors. Again, with 
increase in magnetization the q domains themselves assume the same zigzag 
shape as the p structures, as in fig. 3(f), a most convincing point. ‘Thus, the 
approximate magnetizations of the p structures are fairly well established, though 
their detailed sub-surface shape is not known. 

We therefore see that in devising a closure domain scheme for the (011) 
surface the structure must satisfy the conclusions reached above with regard to 
the p structures and the vertical field patterns of fig. 4, and must roughly agree 
with Néel’s predictions particularly in regard to the q domains. We tried several 
theoretical schemes of which the following two appear to us to be the most 
satisfactory. 

The first, shown in the simplified diagram, fig. 5, contains the q domains 
predicted by Néel without the intervening p domains. ‘The flux emerging from 
the q domains is therefore carried partly antiparallel to the applied field by a 
regular but complicated arrangement of domains magnetized approximately in 
the [100] and [100] directions, and partly parallel to the field outside the crystal 
between adjacent free poles of opposite sign. We may thus explain the wide 
black deposits of colloid between colloid-free areas, and the asymmetric 
distribution of colloid found when the colloid particles are polarized 
perpendicular to the surface, as in fig. 4, a close study of which shows that 
asymmetric colloid distribution takes place on the boundaries themselves. In 
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fig. 5 it is seen that the q domain boundaries have the same sign for the surface 
free poles as their adjacent main domains, so that there is a tendency for the 
colloid to deposit on alternate boundaries with a spacing D. Experimentally, the 
p structures also have an asymmetric distribution on their boundaries, but we 
cannot at present interpret it. 

Of course, this simplified structure of fig. 5 does not account for the beautiful 
lace patterns, but an even more serious objection to it is that large magnetostatic 
energy is associated with the free poles in comparison with that associated with 
Néel’s p domains, unless the extent of the surface occupied by free poles is 
restricted. An elementary calculation (Mee 1951) shows that 15% of the surface 
should not be exceeded in this case. 

A second suggested structure is shown in fig. 6, in which the triangular 
closure domains are added to reduce the magnetostatic energy, although some 
magnetostatic energy is required for their formation, and the difference in 
energy between the two systems of figs. 5 and 6 may not be great. In addition 
to giving the same powder patterns, the second theoretical structure suggests 
why the q domains at high magnetizations assume the same zigzag form as the 
colloid-free p structure at small magnetizations. At high magnetizations the 
angle between the main domain vectors and the (011) plane is similar to that 
between the triangular closure domain vectors and this surface at low fields. 
It is therefore possible that closure domains magnetized parallel to the surface— 
q domains and colloid-free p structures respectively—will have the same general 
shape, since they occur on domains whose orientations are similar. 


bil, 


Fig. 5. Closure structure on (011) plane. Fig. 6. Closure structure on (011) plane. 


To sum up, we see that Néel’s p domains are replaced in our scheme by a 
fine structure of domains magnetized approximately in the same direction as the 
adjacent q domains. If we denote the angle between the direction of the 
q magnetization vector and the [011] direction by ¢, calculations (Mee 1951) 
on the same lines as those of Néel show that ¢ may vary between large limits for 
a comparatively small change in the energy of the q domain. It is therefore 
reasonable to consider that the zigzag shape of the q domains at high 
magnetizations is due to rotation of the q vectors from the [100] direction. 


§5. QUANTITATIVE MEASUREMENTS 
Measurements were made of the domain spacing D as a function of H on 
both the (100) and (011) surfaces, both with the original crystal and with the 
crystal when systematically reduced in width. The specimen was exposed to a 
uniform magnetic field by means of an electromagnet with specially shaped 
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pole-pieces. ‘The effective field acting on the specimen was measured by a 
calibrated magnetic potentiometer whose feet rested on the side of the specimen 
opposite the (011) plane surface investigated. The potentiometer had feet 
6 mm wide and 2 mm thick and the distance between its feet centres was 7 mm. 
It consisted of a semi-circular ebonite former uniformly wound with four layers 
of double silk covered No. 46 s.w.c. copper wire. Using a sensitive galvanometer, 
fields of about 4 Oe were measured with a possible error of 3°/, while with fields 
above 20 Oe the error was about 1%. 

The average spacing of the powder lines was found by counting about 
20 lines in the microscope image projected on a ground-glass screen and measuring 
their total width with a steel ruler. The magnification was later determined by 
replacing the specimen by a graticule. The specimen was not removed from the 
pole-pieces for washing between readings under different field conditions, but 
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Fig. 7. Variation of domain spacings on (100) and (110) planes with field. 


the surface was cleaned 7m situ with a water jet; this procedure had many obvious 
advantages. 

Prior to each determination the specimen was first demagnetized. A field 
of say 10 Oe was then applied and reversed many times to obtain cyclic 
conditions. A drop of colloid was now placed on the (100) surface and the line 
spacing D measured. ‘The colloid was then washed away with the water jet. 
Magnetic field measurements were then made by noting the deflection of the 
galvanometer in the potentiometer circuit when the field current was reversed. 
The whole magnet assembly was next turned through 90° to bring the (011) 
plane horizontal and under the microscope. A drop of colloid was placed on 
the (011) plane and the spacing width D’ on this plane was also measured as 
described above. Later it was found possible or preferable to take all readings 
up to saturation on one plane before taking any readings on the other. 

In fig. 7 the encircled points and the filled circles represent the values of the 
domain spacing as a function of H obtained from measurements on the 0-6 cm 
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thick crystal for the (100) plane and (011) planes respectively, while the broken 
curve is that calculated from Néel’s theory based on fig. 1(b). The crystal was 
then systematically reduced in width by careful grinding, and further sets of 
measurements on the two planes were taken for crystal widths of 0-43, 0-37, 
0-24 and 0-13 cm, and curves similar to those of fig. 7 were obtained. In all 
cases there was very good agreement between the two sets of measurements for 
the domain spacings D and D’. 

The values of the spacings D on the (100) plane for chosen fields are plotted 
as a function of \/L, where L is the specimen width, in fig. 8, together with a 
theoretical curve from Néel’s theory. The two points in the top right-hand 
corner of the figure are approximate only, as they were obtained by applying the 
field in the [011] direction, i.e. at 90° to the field direction in the preceding graphs, 
thus giving an effective value of L =1-21 cm, but the values of H could not be as 
accurately determined in this case. 
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Fig. 8. Variation of domain spacing with width of specimen. 


A close examination of the photographs of lines on the (100) plane shows 
that the colloid striations found on the sides of these lines changed direction 
as the applied field was increased, and it was possible to measure their directions 
on the original photographs. Now, it follows from the work of Williams, 
Bozorth and Shockley (1949) that colloid striations always form perpendicular 
to the main domain magnetization vector. Consequently, the angle @ between 
a striation and a colloid line should be equal to the angle between the domain 
vector and the resultant magnetization which is along the [011] direction. We 
can therefore measure @ directly. 

The angles between striations and lines were measured by projecting a 
magnified image of a photographic negative on to a screen and tracing the colloid 
lines and striations. ‘The mean of about 50 values of 6 were taken for each of a 
series of negatives for a crystal in different fields. D was plotted against cos 0 
and compared with the theoretical curve based on Néel’s calculations; the two 
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curves are reproduced in fig. 9. The general agreement in outline of the two 
curves provides rather striking, fundamental support for Néel’s theory, although 
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Fig, 9, Variation of domain spacings with angle between domain magnetization vector and field. 


as in the case of the (D, H) curves the experimental curve lies well above the 
theoretical curve. 


§6. CONCLUSIONS 


The general shape of the experimental (D, H) curves agrees well with that 
predicted by the Néel theory, when allowance is made for the anisotropy of 
silicon-iron, and we have the supporting evidence of the (D’, H) curves. 
Although the agreement between the (D, H) and (D’, H) curves is reasonably 
good, the (D’, H) curves always lie above the (D, H) curves at high fields, but 
here the sensitivity of the powder pattern method is lower. Moreover, line 
splitting occurred more frequently on the (100) plane than on the (011) plane and 
this tends to depress the D values with respect to the D’ ones. The shape of the 
(D’, H) curve conforms better the Néel’s theoretical curve although no minimum 
was ever observed. 

But the most outstanding feature is the fact that the experimental values of D 
are often as much as three times as great as the theoretical values. We had 
expected that the relatively small discrepancy observed by Bates and Neale 
would be removed when a bigger and better crystal and improved technique 
were available. In the measurements of Bates and Neale the discrepancy varied 
with the crystal width and they did not find a linear relation between D and \/L, 
but in the present work we have very good proportionality between D and \/L. 

We think that all these discrepancies must be due, at least in part, to the 
complicated p structure observed, for which the energy must be smaller than W, 
in the Néel expression D=[yL(W,+ W,)/W,W,}}”. A suitable reduction in 
W,, would, of course, give higher theoretical values for D and a greater slope to 
the theoretical (D, \/L) graph. It is equally possible that a modified closure 
structure on the very narrow (011) faces of the specimens used by Bates and Neale 
would account for their results. We hope in a later communication to give a 
mathematical treatment of some of these cases. 

It would be interesting to extend the above measurements to single crystals 
of very pure iron or, indeed, to any similar ferromagnetic crystals with suitable 
anisotropy constants, and to a specimen of iron or silicon-iron with a circular 
instead of a rectangular cross section, since, according to Néel, a two-phase 
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structure exists when such a specimen is magnetized along the axis of the 
cylinder. The closure structure on the circumference of the cylinder might be 
very interesting, since Néel predicts that a continuous change from p to q 
structures will exist. 
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A Study of Surface Closure Domains by the 
Powder Pattern Technique 
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ABSTRACT. ‘The powder patterns in the neighbourhood of cavities and inclusions on 
the surfaces of single crystals of silicon-iron have been studied, and the photographs 


obtained are evidence for the essential correctness of current views on closure domain 
structures. 


$1. 5INTRODUCTION 

N a series of important papers, Williams, Bozorth and Shockley (1949) and 
] Williams and Shockley (1949) have described the ‘dagger’ and ‘tree’ 

patterns which may be formed by the Bitter powder technique on the 
carefully polished surface of a single crystal of silicon-iron, and have given 
convincing explanations of the modes of formation of special patterns in terms 
of ferromagnetic domain theory. We have in the course of recent researches 
made many similar observations and have confirmed the essential correctness of 
these explanations. In addition, we have made a careful study of the closure 
domains found in the neighbourhood of cavities or inclusions in the crystal 
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surface, and of the way in which such closure domains are modified by the 
presence of a main 180° boundary wall in the neighbourhood of the imperfections. 
A full description of the colloid or powder technique and references to other 
work will be found in papers by Bates and Neale (1950), Bates (1950), and 
Bozorth (1951). 


§2. A STUDY OF DOMAIN MOVEMENTS 

In fig. 1 (Plate I) is reproduced an interesting series of pictures of powder 
patterns formed on the (100) plane of a 3° silicon-iron crystal specimen, some 
13 mm long, 6 mm wide and 6 mm deep, which illustrate the domain formations 
around irregularly shaped cavities and the way in which they are modified as a 
main 180° boundary moves from one cavity to another. In fig. 1 (a) we see that 
a main boundary wall appears to be attached somewhere near the lower edge of 
the cavity system. The actual placing of this boundary was achieved by applying 
a small magnetic field parallel to the surface of the specimen and to the long edge 
of the paper. T'wo distinct types of closure structure are visible around the two 
main cavities, which appear as large irregular blobs inthe picture. The structure 
around the cavity which is further to the left in the picture is of the type shown 
diagrammatically in that part of fig. 3 which lies to the left of the main domain 
wall, while that around the cavity to the right is of the type which lies to the right © 
of the main domain wall in fig. 3. In both cases the possible directions of the 
domain magnetization vectors are given by the arrows. 


a | 


Fig. 3. Schematic representation of fig. 1 (a). Fig. 4. Closure domains around an inclusion 
(after Néel). 

On increasing the magnetic field the boundary moved slowly from the 
position shown in fig. 1 (a) to reach the positions successively shown in figs. 1 (4) 
and 1(c). In fig. 1(b) the closure structures to the left are extended and still 
attached to the boundary, but in 1 (c) they have broken away. In fig. 1 (c) the main 
boundary now attempts to pass through the cavity on the right, and the closure 
structure round this cavity changes from the earlier ‘dagger’ formation to a 
single type which was predicted by Néel (1944) and which is shown 
diagrammatically in fig. 4 (a very good example is seen in fig. 2(a)). We note 
that the lower closure domain of the cavity on the left in fig. 1 (c) is still attached 
to the main domain wall, but is complicated by being joined to another closure 
domain system originating from a cavity somewhere on the left of the picture. 
On the upper side of the cavity the domain structure is now that really 
characteristic of a large cavity isolated within a large domain, and it bears a 
striking resemblance to the structure, shown diagrammatically in fig. 5, which 
was predicted by Néel for such a case. The twin-triangle structure is clearly 
pronounced in both figs. 1(c) and 1(d), to which we next refer. 
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Fig. 1(d), shows what happened when a slight reduction in the applied field 
was made; the closure structure is seen to extend up to the main wall from the 
cavity on the right. On again increasing the field the wall moved across the 
cavity with the newly formed structure attached to it, as shown in fig. 1(e). 
Finally, with still greater field, the wall moved away from both cavities as in 
fig. 1(f), where we now see that the structure around the cavity on the left is 
definitely that of fig. 5, while the structure around that on the right is somewhat 
similar thereto. 

Some experiments were also carried out on the (011) surface of the above 
specimen. In the demagnetized state the powder deposits showed that on this 
surface the domains were magnetized along the [100] and [100] directions. Some 
six 180° boundaries were irregularly spaced over the surface, and each domain 
was covered with striations of powder deposits perpendicular to the direction of 
domain magnetization. 

In fig. 2 (6) are shown some dagger-shaped domains of reversed magnetization 
close to three comparatively large cavities. If we picture the magnetization 
vector of the main domain in which they are situated to be directed downwards 
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Fig. 5. Closure domains around an isolated inclusion (after Néel). 


parallel to the edge of the paper, then the vectors in the daggers must be directed 
upwards in order that a minimum of free poles shall occur on the boundaries of 
the shallow daggers. Such dagger domains must decrease the magnetostatic 
energy of the free poles which exist on the sides of the cavities. Hence, it is. 
probable that the complete structure consists of dagger and small triangular 
closure domains similar to the system depicted in fig. 5. The total energy of 
the system is thus reduced, since the extra energy located in the walls of the 
daggers must be less than the change in magnetostatic energy which results from 
their formation. 

Applying a field of a few oersteds parallel to the direction of magnetization 
of a domain caused the crystal to become saturated, and it was found that the 
size of the daggers was independent of the field, thus indicating that the change: 
in the energy associated with them in the presence of the field was small compared. 
with the change in magnetostatic energy caused by their initial formation. 

It is considered that daggers appeared on only one side of a cavity in fig. 2 (6): 
because the surface in this case was not a true (011) plane but was inclined slightly 
thereto. Thus only one of the daggers, which are orientated along the [100] 
direction, intersects the surface. Following upon much more accurate grinding 
of the surface, checked by orientation measurements, it was found that at an 
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undulation in the surface dagger-shaped domains of reversed magnetization 
were obtained as in fig. 2(c). The vertical line, about one-third of the width 
of the picture from the left, is a 180° boundary. To the left of this boundary a 
main domain may be considered magnetized in the direction from bottom to 
top of the printed plate; to the right of the boundary a main domain may be 
considered as magnetized from top to bottom. It is then reasonable to consider 
the lozenge-shaped daggers on the left of the boundary as magnetized in the 
top to bottom direction, and those on the right to be magnetized in the bottom 
to top direction. Again, on magnetizing the crystal along a [100] direction, the 
180° boundary in fig. 2(c) was displaced laterally across the picture, but no 
change in position or appearance of the daggers was observed. 

Similar experiments were made on a 2:6%, silicon-iron single crystal sheet 
specimen some 0-5 mm thick, with the (110) plane inclined at about 2° to the 
surface, and the edge of the surface parallel to a [110] direction. Very regular 
deposits were observed on the demagnetized crystal as shown in fig. 2(d). The 
black lines running parallel to the long edge of the plate are 180° boundaries 
between domains magnetized alternately in opposite directions parallel to that 
edge, for the powder striations perpendicular to such directions are clearly seen. 
In addition, there is an interesting arrow-head formation on the 180° walls. 


Fig. 6. Schematic explanation of arrow-head formations. 


On applying a small field perpendicular to the surface the arrow-heads were 
more clearly manifested, since certain (alternate) domains collected colloid, 
whereas the arrow-heads situated in them did not, while in the adjacent, 
comparatively colloid-free, domains the arrow-heads collected colloid. This is 
clearly shown by figs. 2(e) and 2(f). ‘The asymmetric distribution of colloid on 
the main domains is explained by the fact that the [001] direction does not lie 
exactly in the crystal surface, since that surface is inclined to the (110) plane 
about the [110] direction. Hence, with no external applied field, the neighbouring 
domains have an even distribution of positive and negative free poles and the 
distribution of colloid is symmetrical. When a vertical field is applied, however, 
the domain vectors rotate slightly towards the field direction and the surface 
density of free poles is no longer uniform, successive domains now having greater 
or lesser free pole density. Moreover, the colloid particles are now polarized 
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perpendicularly to the surface in such a way that they are repelled from domains 
of lesser free pole density and attracted to those of greater free pole density. 
This explanation also covers the colloid distributions on arrow-head formations. 

Because the main domain structure consists of 180° domains, and because 
the specimen surface is slightly inclined to the easy direction along which the 
domains are magnetized, it is natural to think that the arrow-head formation 
represents some kind of closure structure similar to the ‘trees’ observed on the 
(100) plane of a silicon-iron crystal. Very similar arrow-head patterns have 
recently been observed on a (110) plane of a nickel crystal by Yamamoto and 
Iwata (1951) for which they give a ‘tree’ pattern explanation. But in the case 
of silicon—iron the [001] direction is the only direction of easy magnetization 
lying in the surface, and therefore the arrow-heads must be magnetized along 
a more difficult direction of magnetization. This direction can be estimated 
roughly by measuring the angle, «, between the sides of the arrow-heads and the 
[901] direction, since « must be half the angle between the magnetization vectors 
of the arrow-heads and the [001] direction if the condition of continuous magnetiza- 
tion over the boundary is obeyed. The mean value of «, as measured from the 
photographs, was 32°, so that the domain vectors in the arrow-heads must make an 
angle of about 64° with the [001] direction. Consequently, the domain structure 
for the arrow-heads may be depicted as in fig. 6. This is, in effect, a statement 
that a kind of ‘tree’ structure has been observed on the (110) plane. 

In our view, the above series of pictures provide conclusive evidence of the 
essential correctness of the closure domain concept. 
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LETTERS TO THE EDITOR 


On the Radiative Correction to the Coulomb Scattering of 
Electrons 


The result obtained by Schwinger (1949) for the radiative correction to electron 
scattering by nuclei has been checked by the use of the method of Feynman (1949), and 
agreement with Schwinger’s formula has been obtained. However, there is some doubt 
about the method of Schwinger’s proof. 

Using the Feynman approach we have to calculate the matrix elements corresponding 
to the first and second Born approximations and to the first radiative correction. The 
second Born approximation must be included as it is of the same order in e?/fic as the 
first radiative correction. The inclusion of the latter leads to the well-known infra-red 
divergence which is exactly cancelled by the infra-red divergence arising from inelastic 
scattering. To obtain the correct cancellation it is important to use the same limiting 
process in both cases. This is easily achieved, if a vanishingly small photon rest-mass 
(the Feynman A,jn) is used (Dalitz 1951b), which is relativistically invariant. Schwinger 
uses a minimum photon momentum &,jn, but he applies the kj, cut-off in the radiative 
correction after first changing the origin in momentum space (Schwinger 1949, eqn. (1.88)), 
while no such change of origin is made in the bremsstrahlung calculation. In fact, the 
divergent integral in the radiative correction is, in the Feynman notation, 


fi d‘k al al ( d‘k 
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where p=xyp,+x(1—y)p2, and the change of origin (Schwinger’s eqn. (1.88)) is equivalent 
to the transformation (Feynman 1949, after eqn. (11 a)) 


kas kp. apr Ke ee RAP ae. oc kedcree (2) 
which transforms (1) into 
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where the vector K is (k,, ks, kj). Feynman makes this integral finite by giving the photon 
an infinitesimal mass Ain; Schwinger makes it finite by cutting off |K| at an infinitesimal 
momentum kyin- However, because of the transformation (2), which is clearly a function 
of the integration variables x and y, this ky,j, has no direct physical meaning and is certainly 
not identical with the kj, employed in the corresponding bremsstrahlung. integral, 
which is a cut-off on the true momentum. ‘Thus different regions are cut out of the 
momentum space in the two cases. The infinite parts of the two cross sections still cancel, 
but the remaining finite part comes out quite incorrectly. The reason that Schwinger’s 
final result is nevertheless correct is due to the fact that a second error exactly cancelling 
the first is introduced in his eqn. (2.89), the right-hand side of which should be zero 
(Ward, private communication). 
Further, it can be shown that the Rpjy used in the radiative correction is, most 
unexpectedly, relativistically invariant. On comparing the result of integrating (3) which 
is, using Schwinger’s method and Feynman’s notation, 
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with the corresponding expression, using the Amjn process, which is (Feynman’s 


eqn. (23 a)) 
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we see that log Amin=log 2kmin—1 at all energies. In contrast with this, the relationship 
between Amin and the kmin of the bremsstrahlung calculation is (Dalitz 1951 b) 


1 
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which reduces in the non-relativistic limit to log Apin=log 2kmin—% so that this Rmin 1 is 
not relativistically invariant. Here sinh d=f(1—f?)-1/? sin36, B=v/e and @ is the 
scattering angle. 

The cross section for essentially elastic scattering, i.e. scattering where the 
bremsstrahlung photon has an energy up to AE where AE < W the kinetic energy of the 
electron, is then 

do Ze" : 
— (0, AE)= {| ———_—_——~. } (1—8?)(1 —? sin? $0)(1 +0,—9y), 
qo. (9, AE) (oa sin? a) (1—f?)(1 —B? sin? 38)(1 +6, — dp) 
where (McKinley and Feshbach 1948, Dalitz 1951 a) 
§,(0) =7Ze?B sin $0 (1 —sin $0)(1 —f? sin? $0) 


and 


5,(0, AE <W)=(e2/n) {20 —2¢ coth 24)(1-+log 2AE/m) +¢ tanh¢ 


4 
+ 2(1—¢ coth ¢)(1—4 coth? pms Z + glee : a 2¢ coth 2¢ log ie 
262 sin? 40 
Seater i) cosech 


—(1—B?) cosh 2¢ cE dl :(*E 3(1—BC) os log 3(1 a) } 7 
cos 9/2 B sin $0(£?—cos* 38)"? \  1+BC 1—p¢ 
Here natural units have been used. 

Lastly we consider the question of experimental verification. "The second Born 
approximation gives correct results (MacKinley and Feshbach 1948) for Z less than 
approximately 15, so that we must restrict ourselves to light nuclei. For these, any deviation 
from pure Coulomb scattering due to the finite size of the nucleus is negligible at energies 
less than 10 mev (Elton 1950, Acheson 1951). A further restriction is that 5, must be small, 
as otherwise the next term in the perturbation series becomes important. The table of the 
radiative correction dz (as a percentage) for various energies W of the incident electron and 
AE for the bremsstrahlung photon has been worked out bearing these considerations in mind. 


W (Mev) 25 4-0 : D5 
6 (deg.) 45 90 135 45 90 ASS 45 90 135 
AE (kev) : 
10 © AS Aare, 0:95 99-9 Wilss} IW re alte) UGS 
ZS Se) OO) il Oi Des} 10°5 135 14:8 
50 BI, SAO Sy) 4-7 6:8 US EO = aliloy/ 12°8 
100 Des AyD oH Sitey Lys)os) 6-4 7°6 8) 10-8 


Our thanks are due to Professor H. S. W. Massey and Mr. L. Castillejo for many 
helpful discussions and to Dr. J. C. Ward for informing us of the error in Schwinger’s 
eqn. (2.89), which started us on this investigation. 


Wheatstone Physics Department, L. R. B. ELTon. 
King’s College, London. 
University College, London. 


H. H. Ropertson. 
9th November 1951. 
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Properties of Films of Non-Metallic Antimony 


A recent investigation of photoconductivity in the elements (Moss 1951) has 
shown the elements B, C; Si, P, S; Ge, As, Se; Te and I to be photoconductive. 
It is now known that gray tin is a semiconductor (Busch, Weiland and Zoller 1951) 
and thus, presumably, should be photoconductive under appropriate conditions. 
Thus in the fourth row of the periodic table there is a gap between Te and Sn, and 
to retain an ordered progress in the table requires the filling of this gap by a 
semiconducting form of antimony. 

From the behaviour of the other elements we are led to expect this semiconducting 
form to have an activation energy between those of Te (0:37 ev) and Sn (0:1 ev), 
and thus we might expect it to exhibit photoconductivity in the 8 region of the 
spectrum. It was to check this interesting speculation that the present work was 
carried out. 

From experience with arsenic, which gives semiconducting layers if evaporated 
on to substrates not much hotter than room temperature, attempts were made to 
produce semiconducting antimony layers by evaporation (im vacuo) on to substrates 
at 195°x or 90°K. Semiconducting layers were obtained in this manner which 
changed irreversibly to the metallic form at temperatures near 0° c. 

Measurements of the temperature dependence of resistance were repeatable 
provided the layer was not allowed to reach too high a temperature. Typical results 
are shown on the usual (log (resistance), reciprocal temperature) plot by curve A of 
the figure. Although the points do not all fall on one straight line, representative 
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A. Resistance-temperature curve. B. Spectral sensitivity curve. 


lines may be drawn through parts of the curve to obtain estimates of the activation 
energies. The values found are: low temperatures, 0:07 ev; high temperatures, 
- 0-11 ev. Two other layers gave very similar values for the same temperature regions 
as the above. 

The resistance of the layers decreased on irradiation. Using a double 
monochromator with radiation interrupted at 33c/s, it was possible to measure 
spectral sensitivity curves for some of the layers. With one cell fitted with a 
KRS5 window it was possible to extend the measurements to wavelengths greater 
than 15. The results for a layer temperature of 195°K are shown in curve B of 
the figure. The curve is characterized by a marked rise in sensitivity (~36: 1) 
in the 4-9 region. Similar behaviour was found in all cells measured, increases of 28:1, 
33 : 1 and 37 : 1 being recorded for three other cells fitted with periclase, KRS5, and 
silver chloride windows. 

10 2 
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For wavelengths greater than 9y the sensitivity is relatively constant, with the 
possibility of a significant fall near the limit of measurements. As there is no 
well-defined drop of sensitivity at long wavelengths it is not possible to conclude 
from the spectral sensitivity curves whether the layer is acting as a photoconductor 
or bolometer. 

Measurements of response times of layers gave values near 10~* sec, which might 


equally apply to photoconductors or thin, backed bolometers. In one set of | 


measurements the response time increased by 30% on cooling from 195°k to 
90°, in another set it decreased by 40% for the same cooling. 

It is thus established that layers of antimony can be prepared which show 
semiconducting properties, thermal activation energies of 0:06 to O-1lev being 
indicated. A paper by Suhrmann and Berndt (1940), which has recently come to 
the notice of the author, gives values of 0:055-0:18 ev. Layers have been made 
which were sensitive to radiation of wavelengths as great as 15p. 

It is probable that with this investigation of non-metallic antimony (together 
with the elements mentioned above) all the possible photoconducting elements have 
now been investigated. 


Acknowledgment is due to the Chief Scientist, Ministry of Supply, and to the 
Controller, H.M. Stationery Office, for permission to publish this letter. 


Telecommunications Research Establishment, TS. Miess 
Gt. Malvern, Worcs. 
11th October 1951. 


Buscu, G., WEILAND, J., and ZoLLer, H., 1951, Semi-conducting Materials, ed. H. K. Henisch 
~ (London : Butterworths Scientific Publications). 
Kress, H., 1951, Semi-conducting Materials, ed. H. K. Henisch (London: Butterworths Scientific 
Publications). 
Moss, T. S., 1951, Proc. Phys. Soc. A, 64, 590. > 
SUHRMANN, R., and BERNDT, W., 1940, Z. Phys., 115, 17. See also Krebs, H., 1951. 


The Dissociation Energy of SnO 


Dissociation limits have been obtained by short and apparently reliable extrapolations 
of vibrational levels in excited electronic states of several of the diatomic molecules of the 
group IV (b)+VI (b) elements (for example : SiS, Vago and Barrow (1946 a); GeS, GeSe, 
Drummond and Barrow (1952); SnSe, Vago and Barrow (1946 b)). These extrapolations 
are based on the observation of long v’=0 progressions in the absorption systems 
E<X'x. The outstanding problem in the determination of the dissociation energies of 
these molecules is the identification of the products of dissociation of the states E. It 
seems most probable that MX(E)— M(°P)+X@P) (Vago and Barrow 1946 a), but the 
multiplet separations in the *P atomic states are by no means negligible in all cases, and 
it is desirable to attempt to identify the products still more closely. The object of this 
note is to point out that, at least for SnO, a comparison of the thermo-chemical dissociation 
energy with the spectroscopic limit leads to a tolerably certain identification of the atomic 
products of dissociation of the state E. 

The dissociation limit obtained from the E-X system of SnO corresponds to 
Do=131+, kcal (298°) (Eisler and Barrow 1949). The thermo-chemical value is 
obtained from the heat of the reaction 


Sn,+ SnO,,=Sn0, se cneeri (1) 


for which Brewer and Mastick (1951), from an analysis of the data of Vesselovsky (1943) 
derive AH.5g=136-5+0-7 kcal/mole, combined with the following values : 


AH gy 9, kcal/mole Reference 
Sn,+ Og,=SnO0,, — 138-8 N.B.S. 1947 
Sn,=Sng 62 Baughan 1951 (private communication) 


Oo,=20, 118°3 N.B.S. 1947 


a 
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We obtain Do’(298)=122-, kcal/mole. The value derived for the excitation energy of the 
products of dissociation of the state E is thus F,—131-,—122-,=9-, kcal. The energies 
of the relevant atomic states are given below (values in kcal). 


Sn O 
8% =P lB =P eD: Bs aE Rs 
24-6 9-8 4:8 0:0 45-4 0:65 0-45 0-0 


‘The assumption that the products are Sn(#P)-} O(?P) is strongly supported by these results. 
Moreover, the calculated value of Fy, 9°, kcal, fits very well with the excitation energy of 
Sn(*P.), 9-8 kcal. (The multiplet separations of O(?P) are of course far smaller than the 
experimental uncertainty in E,.) The present evidence is therefore consistent with the 
process SnO(E) > Sn(?P,)+O(@P). 

The validity of this argument depends primarily upon the accuracy of the values of the heat 
of sublimation of tin and of the vapour pressures of SnO used in the determination of the heat 
of reaction (1). The uncertainty in Lg, is unlikely to be greater than + 2 kcal (Baughan, 
private communication); only the one set of observations on the vapour pressure of SnO is 
known to us ; the measurements are self-consistent, but there is no means of assessing the 
probable error in the value derived for AH of (1). We must therefore conclude that although 
the agreement between the thermo-chemical and spectroscopic values, on the basis that 
SnO(E)> Sn(?P,)+O(?P) is at the least highly suggestive, we cannot regard the present 
identification of the products of dissociation of the state E as more than tolerably certain. 

The best value available for D)’(0° K) of SnO is then 121 kcal or 5-2; ev. 


Note added in proof : 

' Brewer and Mastick (1951) derive Dy”(0°K) =131., kcal/mole. The discrepancy between 
our values arises because they have used Lg,=70 kcal, and because in correcting from D)” 
from 298° to 0°xK they have added 1:4 kcal instead of subtracting 1 kcal. 


Physical Chemistry Laboratory, G. DRUMMOND. 
Oxford University. R. F. Barrow. 

27th October 1951. 
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Line Widths in Ferromagnetic Resonance 


In recent accounts of the theory of ferromagnetism (Bloembergen 1950, Kittel 1951, 
Van Vleck 1950, 1951) attention has been drawn to the difficulty of explaining the line widths. 
Various possible mechanisms have been proposed, and while some of these do give appreciable 
contributions, the different authors have come to the conclusion that the present known 
mechanisms are insufficient. It is the purpose of this letter to describe another possible 
mechanism. 

It is probably best described by means of an example. Consider the case of nickel. ‘The 
free atom has a ground state which, apart from inner shells, is built from 3d’ 482. In the 
metal the wave functions of electrons on different lattice sites overlap and two different 
approaches, the Heisenberg and the collective electron models, have been used in attempts to 
describe the actual situation. Ferromagnetic resonance gives a g-value of 2:23 at 24° c, and 
it is usual to explain this as arising in some way from an electron which has nearly a free 
spin. On the other hand, measurement of the saturation magnetic moment gives an 
average magnetic moment of 0-606 Bohr magneton per nickel atom. ‘This latter result is 
easily explained on the collective electron theory of ferromagnetism, while the former result 
has been shown to be a consequence of the Heisenberg approach. Suppose now that, instead 
of using the conventional collective electron model to explain the saturation magnetic 
moment, we postulate that 0-606 of the total lattice sites contain unpaired holes (i.e. are d°) 
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and that the remainder are d!°. (The remaining electrons are paired off in a conduction band.) 
Then the saturation moment is given correctly and the resonance would arise from the d® 
sites, giving g~2. The arrangement of d® and d”° sites will presumably fluctuate: an electron 
might move from a d” site to a d? or, alternatively, a d® site might trap a conduction electron 
and become d!°. Thus one has a model in which magnetic and non-magnetic sites are 
continually interchanging, but in such a way that on the average 0-606 sites are magnetic. If 
the lifetime of a given magnetic site is long compared with a Larmor precession time, this 
model will give a width to the resonance line which will persist down to T=0, and will be of 
the order of the field set up at one lattice by a magnetic ‘ hole’ at a neighbouring site. On 
the other hand, if the lifetime is short a very broad resonance might be expected. 

The theoretical determination of the lifetime of a given magnetic site is not easy, but it is 
possible to devise an experiment to test the validity of the model. Suppose resonance 
measurements are made on an isotropic substance which has an integral number of Bohr 
magnetons per lattice site. There will not then be a fluctuation and the model would predict 
very narrow absorption lines. The choice of a suitable substance is very limited; gadolinium 
is the most obvious metal to use as it is known to have an integral number of electrons per 
lattice site, but so far no experiments have been made on it. If such a fluctuation does exist 
in the common ferromagnetic metals it is possible that a lower limit to the lifetime might be 
set by experiments on the scattering of neutrons by the magnetic sites. 


Clarendon Laboratory, K. W. H. STEVENS. 
Oxford. 
29th November 1951. 


BLOEMBERGEN, N., 1950, Phys. Rev., 78, 572. 
KiTTEL, C., 1951, 7. Phys. Radium, 12, 291. 
VAN VLECK, J. H., 1950, Phys. Rev., 78, 266 ; 1951, Physica, 17, 234. 


Some Regularities in Nuclear Energy Levels 


On several occasions nuclear energy levels have been found whose excitation energies 
bear integral ratios to each other. Guggenheimer (1943) carried out a general survey of the 
data then available and attempted, with some success, to explain low-lying nuclear levels 
in terms of a rotator model of the nucleus. 

More recently the accuracy of the experimental results has very greatly increased, and 
anew survey of the data has been made to see whether or not there is now more trustworthy 
evidence for the general existence of simple numerical relations between nuclear excitation 
energies. Results obtained by B- and y-ray absorption methods were rejected as being of 
too low accuracy to be useful; so also were ‘ continuum’ levels of light nuclei since, although 
their spacing is generally known with considerable accuracy, their distance above the 
ground state is not. 

The results of the survey may be summarized as follows: 

(i) Out of 45 nuclei where the positions of two or more energy levels are known 
accurately (involving a total of 175 levels), 38 (involving 140 levels) obeyed the integer 
‘rule’ to within 1 or 2%. The ‘rule’ was also obeyed in a further 16 nuclei out of 23 
where y-ray energies were known but not the positions of the levels. 

After making allowance for the experimental errors, the probability of the random 
occurrence of all the observed integral sequences can be calculated approximately. The 
value obtained is 10~*°, within 2 or 3 orders of magnitude. Nuclear energy levels are 
probably not randomly distributed, so that this calculation is hardly a fair one; nevertheless 
the evidence for the physical significance of the integral ratios seems extremely strong. 

The table gives as an example the experimental levels of 11B (van Patter, Beuchner and 
Sperduto 1951) together with the calculated levels, multiples of 0-1784 Mev, for 


comparison. 
Experimental 2:138 4-459 5-034 6-758 6-808 7:298 8-568 8-926 9-190 9:276 mev 
Calculated 27141 4-460 4-985 6:779 7:314 8:563 8-920 9-098 9:277 Mev 


(ii) if the highest common factor of the nuclear excitation energies (e.g. 0-1784 Mev 
for '}B) is plotted against the nuclear charge Z, then it is found that the majority of the 
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points lie on or near two smooth curves, one for even—even nuclei and the other for nuclei 
having either odd Z or odd N. 

This is especially striking in the case of even—even nuclei where, in addition, the nuclear 
shell structure (Mayer 1950) appears to influence the shape of the curve, since there are 
peaks in the region of nuclei having N=50 and N=82. This is shown in the figure. From 


06 


O5- ® 


Wee Even-Even Nuclei 


MeV 


this result it appears that the highest common factor of the excitation energies of the levels 
in a nucleus may have some physical significance as a nuclear constant. 


Department of Natural Philosophy, P. J. GRANT. 
The University, Glasgow. 
20th November 1951. 
GUGGENHEIMER, K. M., 1943, Proc. Roy. Soc. A, 181, 169. 
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The Spectrum of Lutecium Monoxide 


An intense spectrogram of lutecitum oxide was first obtained by Eder and Valenta (1911) 
from the sulphate burning in a carbon arc source. ‘They recorded the approximate wave- 
lengths of about 30 band heads. The prominent system of bands at 14660 was later obtained 
at the Bureau of Standards on spectrograms of 5-amp d.c. arcs at 220 volts between silver 
electrodes tipped with LuO, and a vibrational analysis of that system was reported by Watson 
and Meggers (1938). Band heads at 45170 and 44093 were also measured by Meggers, but 
as a sufficient number of band heads did not develop in these regions (presumably due to 
the high excitation conditions of their 5-amp arcs) they were ignored. The three systems at 
5170, 44660, A4093 were obtained very clearly on 3-amp arc spectrograms taken in this 


B y 5 og 


(0, 0) (0, 1) (0, 0) (0, 0) 
| || a ene | . if PG gh pau 


4004 4326 4660 5170 5900 A 


laboratory using carbon electrodes soaked in a solution of lutecitum oxide in HCl. The 
complete spectrogram as obtained later from the flame above a carbon rod steeped with 
the same solution and excited by a current of 600 amp at 15 volts is shown in the photograph. 

In addition to the above three systems named «, f, y, the spectrum shows two continua 
8 and ¢ in the red region. On the original negative a hazy sequence of bands could be 
recognized in the feature at A6000, but reliable measurements could not be made. The 
entire group of band systems and continua is attributable to lutecitum oxide, as none of 
the features could be found on similar spectrograms taken with other rare earth oxides 
{Gatterer and Junkes 1945) free from lutectum impurity. 
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Watson and Meggers analysed system f. The vibrational analysis of the other two 
systems is presented in the table. 


New Bands in Lutecium Oxide 


Int. A v v’ v" Int. r v v v" 
2 3978 37 25128:8 com(i2eed)) 10 5171-43 19331-6 y (0, 0) 
a 4094-00 24419-2 i ra , 0) 7 5185-42 19279°5 del (eai)) 
4 4105°81 24348 -9 ay a 5199-65 19226-7 Byl C424) 
6 4241-80 23568°3 1) 8 5217-40 19161:3 (33'S) 
4 4252-49 23509-1 i ah yD) 5 5233-69 19101-7 », (4, 4) 
2 4983-83 20060:9 y (2, 1) 4 5250:65 19040-0 (So) 
5 5019-28 19917-6 % G 3) 4 5418-64 18449-7 Bel 2) 
6 5034-64 19856:9 4) 4 5432-95 18401:1 Sa P48). 
5 5050:52 19794:5 ie ie 5) 8 5448 -63 18348 -2 » (3, 4) 

7 5463-87 18297:0 «~~ ,, (4,5) 


A short sequence of bands, at vy 23 417-6, 23 349°3, 23273-7 and 23 197-0, follows the (0, 1) 
sequence of the a-system, but does not continue it unless a high perturbation is supposed at 
the disappearing (2, 3) band of that sequence. 

The well-known characteristics of the band systems of the rare earth oxides, viz. (i) the 
value 840 cm~! for w¢ of a lower electronic state, (ii) the small difference in the values of the 
vibrational quanta OE the lower and upper states, and (iii) the delimitation of the band 
systems to a few sequences round the (0, 0), are all observable in these systems too. 

The spectrum of the oxide of lutecium with its 5d 6s? configuration may be expected to 
be similar to that of the oxide of lanthanum of the same electronic configuration if the effect of 
the completed 4f orbit in lutecium is negligible. Lanthanum oxide has a ?Z lowest electronic 
state, and transitions to this level from a?II, 87%, c?II, and a D level are known. The 
A°IIT—x?% transition is in the red region and c?J[—x’?& in the region A4360—4415. The 
intense double system 8 and y in the present spectrum appears on this analogy to be the 
c*II—x°X of lutecium oxide. The interval 2114 cm of the c?II term of the molecular 
spectrum of the oxide is of the same order as the 7D ground term interval 2000 cm! (approx.): 
of Lui. The CN band at 4216 overlaps the a-system in the arc spectrograms and it 
is very brief in the flame spectrum shown in the photograph. The dispersion is too small, 
but from what could be judged from the use of the microscope, this system is of the simplest 
structure and may be a?X—*X transition. It is believed that the continua too would develop: 
recognizable band structure under lower excitation conditions. 


Laboratorio Astrophysico, A. GATTERER. 


Vatican City, Rome. S. G. KRISHNAMURTY.- 
30th November 1951. 


Ener, J. M., and Vatenra, E., 1911, Atlas typischer Spektren (Vienna: Hof- und Staatsdruckerei). 
GATTERER, ae and JUNKES, J., 1945, Spektren der seltenen Erden (Rome: Vatican Observatory). 
Watson, W. W., and Meccers, W. F., 1938, Bur. Stand. F. Res. Wash., 20, 125. 


Meinel’s Infra-Red Auroral Bands 


Meinel (1951) has described a band system in the infra-red auroral spectrum which 
has been ascribed by him to the hitherto unknown a*II—x23 transition of N,+. Bands 
which appear to be identical with the most prominent of Meinel’s have been obtained in 
this laboratory from two sources (see also Herman 1951). The first was an electron beam 
tube, similar to one used by Hanle and Schaffernicht (1930) for studying excitation 
functions. It contained nitrogen at a pressure of about 0:01 mm Hg. Accelerating 
potentials of about 200 volts were used. ‘The spectrum was photographed with an 
instrument with dispersion 185 A/mm at 8000 A. The second was a hollow cathode tube 
operated at 440 volts in nitrogen at a pressure of 0-2 mm Hg. The spectrum was 
photographed with a dispersion of 50 A/mm at 8000 A. In both cases, hypersensitized 
Kodak II N plates were used. With the first source, the first positive bands were relatively 
weak, and isolated bands near the positions of Meinel’s (1, 0) and (2, 1) bands were recorded. 
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There was also a very faint band near Meinel’s (3, 2) band at c. 8 298 A. With the second 
source, the new (1, 0) and (2, 1) bands were feebler relative to the first Positive bands, 
and largely obscured by the (7, 6) and (6, 5) heads of the latter. The wavelengths found 
for the maxima of intensity are given in the table, together with Meinel’s data. The 


Wavelengths of Intensity Maxima of Infra-red Bands 


Hollow cathode 220-volt beam Aurora (Meinel) 
78256 7826°5 [8] 7828-2 (1, 0) 
7859-4 7859-1 [11] 7861:9 
7878-0 7877-8 [13] 7879-6 

(7915)* (not recorded) 7912-2 
8053-0 8054-0 [5] 8056-8 (2, 1) 
8090-8 8089-4 [9] 8094-0 
8106-9 8109-8 [10] 8107-0 


* Very feeble. 


numbers in square brackets are rough estimates of intensity from densitometer records, 
after allowing for the background due to scattered light from the filament. The 
wavelengths found here are mostly a little shorter than Meinel’s, but it is thought the 
differences may be due to different excitation conditions, and to a certain arbitrariness in 
locating intensity maxima. The figure is a photograph of the spectrum of a 220 volt beam, 
with an argon+ neon reference spectrum. ‘The spectrum of the beam shows a few O1 lines 
in addition to nitrogen bands and Ni lines; no spectra not due to nitrogen were found on 
the higher dispersion plate. ‘The wavelengths and appearance of the bands are different 
from those of CN bands (Herzberg and Phillips 1948). 


FIRST POSITIVE 


(2, 1) N#? 


tac Gh 
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The writer wishes to thank Professor D. R. Bates and Professor K. G. Emeléus for 
suggesting the use of the two sources, and Dr. O. White and Mr. L. W. Kerr for help in 


building the tubes. 
Physics Department, N. D. SAYERs. 
Queen’s University, Belfast. 

11th December 1951. 
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REVIEWS OF BOOKS 


Die Mathematischen Hilfsmittel des Physikers, by E. MapELUNG. Pp. xx +531. 
4th Edition. (Berlin, Gottingen, Heidelberg: Springer Verlag, 1950.) 
DM. 49.70. 


Madelung’s book was very popular with German students of mathematical physics 
tight from its first publication in 1922, and was doubtless also widely appreciated among 
physicists outside Germany. That the present fourth edition is about twice the siz: of the 
first corresponds with the great increase in the 1equired mathematical equipment of the 
physicist, which has resulted mainly from the development of quantum mechanics. 

One of the many attractive features of the present edition is that it combines a 
systematic compilation of mathematical results and formulae with the connected treatment 
of a textbook. Within the compass of 531 pages the author covers a wide range of 
mathematical topics, and the incorporation of adequate introductory and explanatory remarks 
much facilitates the use of the book. Throughout, the needs of the practising physicist 
are placed first. A resulting sacrifice of some mathematical rigour, which will be regretted 
only by the purist, is obviously necessary to avoid unwieldiness. 

The book is divided into two main parts. The first entitled Mathematics contains 
thirteen sections dealing with almost every aspect of mathematics which is of use to 
present-day physicists. The main subjects discussed are : higher algebra, transformations 
and invariants, vector and tensor analysis, group theory, differential and integral equations, 
variational calculus and statistics. The second part is entitled Physics and constitutes a 
highly condensed treatment of modern theoretical physics. It is extremely well presented 
and will form a ready reference source. The chapter on quantum mechanics deserves 
special mention as it contains a somewhat new and elegant derivation of Schrédinger’s 
equation : statistical concepts are introduced right from the start through the equation of 
continuity for the probability density and its current. To account for the characteristic 
interference effects the equation is split into two linear differential equations in the 
probability amplitude and its conjugate. The derivation is finally based on the exact 
validity of the classical equations of motion for the expectation value of any dynamical 
variable in its time dependence. 

Madelung’s book in its new edition will continue to be of excellent service to all who 
are engaged in physical research and teaching and it should make a valuable addition to 
any physicist’s library. The publishers deserve congratulation for their share in making 
the book so attractive. J. AHARONI. 


Thermodynamics of Irreversible Processes, by S. R. DE Groot. (Selected Topics 
in Modern Physics III.) Pp. xvi+242. (Amsterdam: North Holland 
Publishing Co., 1951.) 35s. 


This highly specialized monograph by the Professor of Theoretical Physics at Utrecht 
is devoted to the systematization of the theory of irreversible processes, and to the 
formulation of the general theory on a logical basis. The treatment is made to depend 
to a very great extent on a theorem due to Onsager which in its turn depends on a 
hypothesis, of a plausible character, about the nature of the fluctuations which can occur 
in a statistical-mechanical ensemble. 

The book is clearly written and contains a number of very useful results in its chapters 


on the application of the theory to heat conduction, diffusion, thermal diffusion and 
thermo-electricity. Jo Howe 


Universal Mechanics and Hamilton’s Quaternions, by O. F. Fiscuer. Pp. vi+ 
356. (Stockholm: Axion Institute, 1951) $10. 


This book describes itself in its preface as ‘“‘a book written by a civil engineer on 
universal mechanics with an attempt to introduce a certain order in its mathematical 
structure by means of the calculus of Hamilton’s quaternions ’’. By universal mechanics 
s meant ordinary mechanics, elasticity, hydrodynamics, aerodynamics, electromagnetism, 
elativity and quantum theory. The description is a just one. 
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After summarizing vector analysis, quaternions are introduced and applied to mechanics 
and electromagnetism. The algebra is treated from a naive point of view which resembles 
the original one of Hamilton, and this is probably an advantage. In the course of these 
applications certain electro-mechanical analogies are discussed. The treatment is by no 
means confined to quaternion methods; for instance, theories of earth magnetism and 
Blackett’s empirical relation are described. Special and general relativity are summarized, 
and Einstein’s 1950 theory is described in 34 pages. Clifford’s numbers, or E-numbers, 
are considered as ‘ double quaternions’, and some of Eddington’s ‘ fundamental theory ’ 
is discussed. This discussion is uncritical of the physical aspects, but wisely remarks that, 
instead of rejection, ‘ It would be more useful to science to give careful thought . . . And 
then attempt to clear up what is dubious ’. Some description is also given of H. G. Kiissner’s 
Principia Physica (Gottingen 1948). This includes a short account of his combinatorial 
theory of arithmetical invariants of mechanical systems, and its applications to fundamental 
physics. 

The notations employed are not always such as to show the power and simplicity of 
quaternions; in the opinion of the reviewer it is a mistake to depart from Hamilton’s use of 
simple juxtaposition to denote the quaternion product. On the other hand, the book 
contains a great deal of information drawn from a wide field ; the bibliography lists references 
to 124 papers and books. This makes the work of some value to anyone using quaternions 
in physical applications. There has been no previ us attempt to fill such a need. It is 
printed by a lithographic method in a rather small print. There is a very satisfactory index. 

C. W. KiLMISTER. 


Fluoreszenz organischer Verbindungen, by TH. FOrsTER. Pp. 312. (Gottingen : 
Vandenhoeck and Ruprecht, 1951.) In paper, DM. 29.50; in cloth, 
DM. 32.50. 


A book specifically concerned with the luminescence of organic compounds is welcome 
at the present time. Recent developments in this field have already been reviewed by 
Professor Férster (Naturwissenschaften, 1947, 6, 166), but his book provides a much more 
comprehensive and up-to-date appraisal. It includes all aspects of the subject, from 
experimental methods to the modern theories of the various phenomena involved. The 
book begins with a short introductory chapter containing necessary definitions of terms and a 
brief historical survey. Chapter II reviews the experimental techniques involved in lumines- 
cence studies, such as measurements of fluorescence intensities, degree of polarization of 
fluorescence, decay of emission and fluorescence spectra. Chapter III introduces the reader 
to the general properties of organic compounds, with particular reference to the aromatic 
hydrocarbons and dyestuffs. Their structural chemistry and the physical conditions 
necessary for fluorescence are discussed. 

After a chapter on the theoretical aspects of light absorption and emission, the transfer of 
energy between molecules, and the conditions for molecular dissociation, the fifth chapter 
again takes up the question of relations between chemical constitution and the probability 
of fluorescence emission, dealing in detail with the aromatic hydrocarbons, with and without 
substituted groups, heterocyclic compounds and dyestuffs. After describing the older 
phenomenological theory which involved the assumption of chromophoric and fluorphoric 
groups, the later physical theory of fluorescence based on wave mechanical ideas is dis- 
cussed. The emission and absorption spectra of organic compounds in vapour phase and in 
solution are next considered with respect to their electronic and vibrational energy states. 

In chapter VII we find a discussion of the absolute yield of fluorescence and its dependence 
on the wavelength of exciting radiation, temperature and nature of the solvent, while the 
next chapter, a very brief one, deals with the decay of fluorescence. It is natural that a 
large space should be given to the important studies of polarization and quenching of fluores- 
cence in solution, with which chapter IX and chapters X and XI are concerned respectively. 
Measurements of fluorescence polarization and its dependence on solution viscosity, tempera- 
ture and concentration of the fluorescent compound furnish important information on 
molecular structure. A knowledge of the dependence of quenching of fluorescence by 
solvent impurities on similar factors to those affecting polarization indicates the nature of 
energy transfer in molecular systems and the relation of fluorescence to other phenomena 
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such as photochemical dissociation and the state of molecular aggregation. With respect 
to the latter it is gratifying to see a section of chapter XI devoted to a discussion of the 
Scheibe polymers (pseudoisocyanin polymer threads), which show such remarkable fluores- 
cence characteristics. 

The text concludes with a chapter on phosphorescence of long duration in organic 
molecules, being particularly concerned with studies by Lewis and his collaborators of 
fluorescein in boric acid and the various theoretical interpretations of phosphorescence in 
condensed systems. 

There are 28 pages of references, including a special list of recent literature on fluorescence 
which, however, is not quite complete. A short, but apparently quite adequate, subject 
index is provided. 

It is apparent from the references in the text that the author is conversant with almost all 
of the recent and past studies in the field by workers in all countries. He has given adequate 
indications of the important contributions made during the last decade by Wawilow and 
other Russian workers, which are, perhaps, not so well known because of their relative 
inaccessibility. . 

Although the author lays stress on the theoretical aspects of his subject, he preserves a 
reasonable balance between them and experimental studies. It is valuable to have the recent 
and more detailed theoretical treatment of such phenomena as quenching and _ polarization 
of fluorescence since these have not been given collectively elsewhere. G. F. J. -GARLICKS™ 


Ferromagnétisme et Antiferromagnétisme (Colloques Internationaux du Centre 
National de la Recherche Scientifique, XXVIII). Pp. 360. (Paris: 
GCNES, 1951) 


The present volume of the series of reports issued by the C.N.R.S. contains the 
49 papers, in French or English, read at the international conference at Grenoble ii 
July 1950. Its contents reflect the marked advances made in the field of ferromagnetism 
and antiferromagcetism since the last conference was held, under the chairmanship of 
the ate Prof. Pierre Weiss, at Strasbourg in 1939. At the same time it is most noticeable 
what wide gaps still exist between well-established experimental facts relating to the 
magnetic properties of even the simplest solids, and a fully satisfying theoretical 
coordination of these facts. 

On the credit side must be placed the enormously increased understanding of the 
domain structure and hysteretic properties of ferromagnetic metals and alloys. This 
advance is due in no small measure to the recent powder pattern experiments, yielding 
such fascinating results, described here-by Dr. R. M. Bozorth (p. 160) and Prof. L. F. Bates 
(pp. 174, 356). Theoretical work on coercivity and related properties is described in an 
admirable review article by Dr. L. Weil (p. 289). The corresponding enormous recent 
advances in the performance of permanent magnets are described by Dr. K. Hoselitz 
(p. 300). An almost completely new, and already flourishing field of research, ferro- 
magnetic resonance, is covered by Prof. C. Kittel (p. 143). The novelty of this field may 
be judged from the fact that spectroscopic splitting factors g which are considerably larger 
than 2 are here the rule. The long way theory has yet to advance, on the other hand, is 
exemplified by the results of precision measurements of the gyromagnetic ratio g’ 
described by Dr. A. J. P. Meyer (p. 156), who finds g’=2-00 for iron and nickel, n 
disagreement with Kittel’s simple rule g—2=2~—g’. Further, as stressed in discussion, 
present ideas as to the interpretation of ferromagnetic anisotropy seem to be invalidated 
by Meyer’s results. The realization of the fundamental importance of anisotropy makes one 
regret that the theoretical aspects of this subject were not discussed more fully at the 
conference. ‘The two main reviews of the present status of the basic atomic theory of 
ferromagnetism in metals and alloys, by Prof. E. C. Stoner (p. 224) and Prof. R. Furrer 
(p. 254), indicate that differences of viewpoint still exist between the various schools of 
thought, roughly described as the ‘ collective electron’ school and the ‘ discrete levels’ 
school. An international conference of this sort may do much to bring these viewpoints 
closer together and to help in overcomir g old-established prejudices. 
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Perhaps even more spectacular advances, fully described in this volume, have recently 
been made in the field of antiferromagnetism and ‘ ferrimagnetism’; the latter term is 
one coined by Prof. L. Néel, the organizer and chairman of the conference and the present 
distinguished doyen of French magnetism. Prof. J. H. Van Vleck (p. 114), in a 
comprehensive review, states the theoretical positi n regarding antiferromagnetism, but 
ends on a note of caution by saying ‘‘ the remarkable thing is that . . . idealized and 
over-simplified descriptions work as well as they do”’. In this connection it is to be 
regretted that the recently developed neutron diffraction technique was not described more 
fully in a separate article. Related to Van Vleck’s review are the various accounts of the 

"present state of experiment, by Prof. H. Bizette (p. 13), Prof. F. Trombe (p. 22) and others. 
Ferrimagnetism is best exhibited by the non-metallic solids grouped under the name 
ferrites. Their properties are described by the late Dr. J. L. Snoek (p. 80), 
Dr. C. Guillaud (p. 91), Dr. E. F. Bertaut (p. 104) and others. The fact that the properties 
of these substances may be understood by regarding them as arising from incomplete 
antiferromagnetism was first stated by Néel; bis already classical papers on this subject 
underlie the discussion to a marked extent. 

It is not possible to mention here more than a few of the papers included in the volume 
under review. Many of the accounts of which no mention is made deal with particular 
properties of particular solids in some detail, and frequently fail to consider these properties 
against the wider background afforded by the excellent more general review articles. Such 
specialized articles could have been more appropriately published as research pap rs in 
the scientific journals. Indeed, some of them are no more than slightly amended reprints 
of papers which have, in fact, already been published elsewhere. The space saved by 
their omission from the present volume might have been more profitably filled by further 
detailed reviews. On the other hand, a formal presentation of papers read at a scientific 
conference can never reproduce the more informal atmosphere usually prevailing on such 
occasions, and, in particular, from all accounts, on the occasion of the Grenoble conference 
during the hot French summer of 1950. E. P. WOHLFARTH. 


Semi-Conducting Materials—Proceedings of the Reading Conference, 10—15th 
July 1950. Ed. H. K. Heniscn. (London: Butterworths Scientific Publi- 
cations, 1951.) 40s. 


The sub-title and the foreword of 14 lines are the only indication that this collection 
of twenty-eight papers on semiconducting materials is, in fact, the Proceedings of an 
International Conference on the subject, and, of the foreword, only four lines specifically 
refer to the Proceedings here presented : 

“The Reading Conference 10—15th July 1950, of which this volume constitutes 
the Proceedings, provided research workers from this country and abroad with an 
opportunity for discussion of the most recent progress. Over 200 physicists and 
engineers took part, among whom it was a pleasure to welcome 80 visitors from 
other countries.’ 

As no discussions of the papers are reproduced, the book appears at first sight as an 
issue of a scientific periodical in which, by some chance, all papers deal with semiconductors. 
Its real significance however, consists in its representing the nearest thing to what used to 
be called a treatise on a scientific subject, i.e. a complete, up-to-date, concise but under- 
standable review. It is to the credit of the organisers that by the selection of the speakers 
this book provides, in fact, a reasonably complete cross section of contemporary semi- 
conductor physics, and to the credit of the authors that they have given well-readable 
accounts with a satisfactory selection of references. 

It is obviously impossible to refer here to each of the 28 papers, and it would be an 
invidious task to refer to a few in particular. ‘The reader will get a good idea of the 
theoretical problems involved, of the classification of materials, of measuring techniques, 
absorption spectra, of some recent ideas and experiments concerning carrier injection, 
transistors and surface phenomena, of the effect of nucleon bombardment, of motion of 
electrons in single crystals and in polycrystalline materials, of the influence of lattice defects 
and, of course, of the results obtained in measurements on a large number of particular 
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materials. Further, there are papers on a particular class of semiconductors which have 
recently been called oxidic semiconductors and on some properties of thin metallic films. 

The most serious omission (not the organisers’ fault) is the lack of reference to recent 
Russian attempts to revise the theoretical basis of conduction in solids, for instance, by the 
introduction of polarons, electron pairs, etc. 

The fact remains that this book will be studied and consulted again and again by every-- 
body working on semiconductors or on allied fields, and it should be on the shelves of any 
physics library keeping a representative collection of periodicals. W. EL 


Calcul Conceptuel ou Qualitatif, by RENE CorpeBas. Pp. 38. (Paris: René 
Cordebas, 1951.) 


Graphical methods of research of ideas and laws in the physical sciences. Tables of 
‘ Idéogrammes ’. 
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(a) Comparison of 1, an ethylene—air flame and 2, a hydrogen—air flame containing a small amount 
of perfluoro-methyl cyclo-hexane (C;F,,)—note strong CH (1, 0) and (1, 1) bands in the 
latter. (Raman quartz spectrograph, Ilford Zenith plate, exposure (1) 30 sec., (2) 10 min.) 

(b) Comparison of 1, a hydrogen—air flame containing ethylene (1:2% ethylene and A=0-36) and 
2, an ethylene-air flame. Note the suppression of the predissociation in the Q branch 
in (1). (Raman glass spectrograph, Zenith plate, exposure (hour, @))lSksecs) 

(c) 1. Large dispersion plate of CH (0, 0) and (1, 1) bands in a hydrogen-oxygen flame containing 

acetylene (2°5% acetylene, A=0-43). 
2. CH (1, 0) band in the same flame. 

Note the suppression of the predissociation in all three bands—Pr marks the line where predissociation 

usually sets in. (Quartz Littrow spectrograph, Zenith plate, exposure 4 hours.) 
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ABSTRACTS FOR SECTION B 


A Coil System for an Air-cored Betatron, by R. LatHam, M. J. Pentz and 
M. BLACKMAN. 


ABSTRACT. The paper describes the design of an air-cored coil system for a small 
betatron and measurements of the field configuration produced by this system. 


Magnetic Electron Microscope Projector Lenses, by G. LIEBMANN. 


ABSTRACT. ‘The results of an earlier paper are extended to the special case of the 
electron microscope projector and similar type lenses. It is shown that the focal properties 
of symmetrical projector lenses with different ratios of gap width S to lens bore D can be 
represented by a single curve. ‘This is then applied to the discussion of various cases of 
practical interest. The calculated results are found to be in excellent agreement with 
measurements by Ruska. ‘The question of image distortion as a function of design 
parameters and lens excitation is then discussed; it is shown that all the types of lenses 
investigated can be operated at a point where they are free from ‘ radial distortion’. 
A consideration of scaling operations shows that there is little difference between the 
optimum performance of lenses of different geometry, but equal focal length. Certain 
results of use in practical lens design are combined in a design chart from which the 
relevant design or operational parameters can be read directly. In conclusion, design 
rules are given for the optimum design of projector lenses covering a wide range of 
magnifications. 


The Electrical and Luminous Characteristics of Short Air Sparks suitable for 
High Speed Photography, by W. G. STANDRING and J. S. T. Looms. 


ABSTRACT. The light emitted by various forms of spark gap is examined, and the factors. 
affecting maximum value and variation with time are discussed. Measurements of the 
voltage drop across the spark help to explain the effects on light output of spark duration and 
form of spark gap. A rough spectral analysis is given for an open gap. 


A Note on the Calculation of Principal Ray Aberration, by R. S. LoNGHuRsT. 


ABSTRACT. An expression for the spherical aberration of the principal rays is derived. 
If the primary aberrations are calculated the principal ray aberration may be found with 
very little additional effort. 


Experiments on the Mechanics of Rubber: I—Eversion of a Tube, by A. N. GENT 
and R. S, RIVLIN. 


ABSTRACT. Measurements are made of the changes which take place in the dimensions. 
of rubber tubes when they are turned inside out. These are compared with the values 
predicted from the theory of large elastic deformations of incompressible isotropic 
materials. 


An Interpretation of the Magnetic Properties of Some Iron-Oxide Powders, 
by W. P. OsMonp. 


ABSTRACT. ‘The observed magnetic properties of fine dispersed powders of the 
ferromagnetic iron oxides are discussed in the light of modern theories of the magnetization 
of ferrites and of ferromagnetics containing non-magnetic inclusions. ‘The individual 
powder particles are shown to be probably of single-domain size, and the relative 
coercivities of y-ferric oxide powders of nearly spherical and of markedly acicular grains 
are qualitatively in agreement with theory, on the grounds of shape anisotropy. After 
allowance is made for volume concentration and small departures from stoichiometric 
composition, the ratio of the coercivities of apparently identical acicular powders of 
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y-ferric oxide and of magnetite is found to be approximately that of the saturation — 


magnetization per ferric ion rather than that of the value per unit volume, and both 


coercivities are lower than theory would predict. It is concluded that magneto-crystalline 3 
anisotropy is also involved. In contrast, relative and actual values of intrinsic induction 


at the beginning of the approach to saturation are in satisfactory agreement with theory, 
as are the relative values at remanence although the actual values here are rather low. 
It is hoped to extend the study to larger departures from stoichiometric composition 
involving partial transformation to the antiferromagnetic iron oxides %-Fe,O, and FeO. 


The Emission from Oxide-coated Cathodes in an Accelerating Field, by D. A. 
WriGuHT and J. Woops. 


ABSTRACT. It is shown that the steep rise of current with accelerating field in a 
diode with an oxide-coated cathode cannot be fully explained either by surface roughness 
or by patch effect. This leads to a theoretical treatment, based on the existence of a 
space-charge zone immediately inside the coating. The charge in this zone is shown to 
vary with applied field and with current density, and with certain coating parameters. 
The variation in the charge leads to a variation in work function, and thereby to a 
dependence of emission on field strength, which is to be combined with the normal 
‘Schottky effect. The effects of partial space charge on the relationship between field and 
anode voltage are discussed, and the diode characteristics of oxide cathodes can then be 
calculated for accelerating fields as a function of the coating parameters. ‘The results 
are shown to be in agreement with experiment, and to confirm values of coating 
conductivity, electron density and mobility as determined directly by conductivity and 
Hall effect measurements. 


The Effect of Mercury on Selenium, by H. K. HENiscu and E. W. Saker. 


ABSTRACT. It is shown that the interaction between crystalline selenium and liquid 
mercury or mercury vapour results in the formation of mercuric selenide. When produced 
in this form, this material is an excess semiconductor of high conductivity. The associated 
volume and surface diffusion have been investigated using radioactive °°Hg. Mercury 
added to selenium before the crystallization process results in material of high resistivity. 
This is ascribed to a reduction of the positive hole mobility. 


Impurity Scattering in Oxide Semiconductors, by E. W. J. MITCHELL. 


ABSTRACT. The oxide semiconductors are characterized by relatively high concentra- 
tions of impurity centres and low electron mobility. It is suggested that electron scattering 
by neutral impurity centres makes an important contribution to the resistivity of these 
materials. ‘The familiar relation log a7 =«-+ fe together with the theory given by Busch are 
discussed and an alternative explanation involving impurity scattering and a dependence of 
« on N is given. ‘This is compared with experimental results for several oxides, including 
measurements by the author on the Fe,(Ti)O, system. 
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